Strain-Tunable Harmonic Responses in Valley-Polarized Bilayer Graphene

Technical Summary: Strain-Tunable Harmonic Responses in Valley-Polarized Bilayer Graphene

Problem Statement
Bilayer graphene (BLG) offers a versatile platform for investigating symmetry-controlled phenomena in two-dimensional quantum materials. However, pristine BLG possesses inversion symmetry, which causes the electric-dipole second-order optical response, such as second-harmonic generation (SHG), to vanish. While previous studies have demonstrated that breaking inversion symmetry via perpendicular electric fields, substrate interactions, or twist engineering can activate SHG, the specific role of uniaxial strain in conjunction with valley polarization remains insufficiently understood. In valley-balanced systems, contributions from the two inequivalent valleys ($K$ and $K'$) typically cancel out due to time-reversal symmetry. The authors investigate how the combined effects of valley polarization (an imbalance in carrier populations between valleys) and uniaxial strain can lift this cancellation and modulate nonlinear optical responses.

Methodology
The authors employ a theoretical framework based on a low-energy effective two-band Hamiltonian that incorporates trigonal warping and strain-induced anisotropy. The system is modeled under uniaxial strain, which introduces a symmetry-breaking term ($w\sigma_x$) into the Hamiltonian, where $w$ is proportional to the difference in strain components ($\epsilon_{xx} - \epsilon_{yy}$).

To calculate the optical responses, the study utilizes the quantum kinetic formalism within the density matrix approach. The time-dependent density matrix is expanded perturbatively to solve for linear ($\rho^{(1)}$) and second-order ($\rho^{(2)}$) interband and intraband coherences. The macroscopic current density is derived from these elements to determine the linear susceptibility ($\chi^{(1)}$) and the second-order nonlinear susceptibility ($\chi^{(2)}$).

Key parameters include:

• Valley Polarization: Modeled by assigning a valley-dependent Fermi energy, $E_f(\nu) = E_f + \nu \Delta E_f/2$, creating a population imbalance.
• Strain: Simulated by varying the strain energy parameter $w$ across values of $\pm 5, \pm 10,$ and $\pm 15$ meV, corresponding to realistic uniaxial deformations (0.1%–2%).
• Conditions: Calculations are performed at $T = 5$ K with a scattering rate $\gamma = 10$ meV, focusing on the degenerate semiconductor regime near the Fermi surface.

Key Contributions and Results

1. Activation of SHG via Valley Polarization:
   The study confirms that in valley-balanced BLG, the second-order response vanishes due to the exact cancellation of contributions from opposite valleys. However, a finite valley polarization ($\Delta E_f \neq 0$) lifts this cancellation, enabling a net, measurable SHG signal. This mechanism relies on the cooperative effect of trigonal warping and valley population imbalance rather than an inversion-symmetry-breaking gap.

2. Strain-Induced Anisotropy and Tunability:
   Uniaxial strain significantly modifies the nonlinear response by distorting the low-energy electronic structure and altering the pseudospin texture.

   • Linear Response: The linear optical susceptibility exhibits a resonant peak at $\hbar\omega \approx 2E_f$. The magnitude of this response is modulated by the strain parameter $w$, with the peak intensity decreasing as $w$ increases from negative to positive values along the zigzag direction.
   • Nonlinear Response (SHG): The SHG spectrum displays a dual-peak structure with resonant enhancements at $\hbar\omega \approx E_f$ (two-photon resonance) and $\hbar\omega \approx 2E_f$ (one-photon resonance). Unlike the linear response, the SHG signal is highly sensitive to both the magnitude and sign of the strain.

3. Mechanical Switching Mechanism:
   A primary finding is that changing the sign of the strain parameter ($w$) reverses the direction of the induced second-harmonic current. For instance, a positive strain ($w > 0$) induces a sign inversion in the nonlinear susceptibility components ($\chi^{(2)}_{xxx}$ and $\chi^{(2)}_{xyy}$) compared to negative strain ($w < 0$). This provides a mechanically controlled switching mechanism for nonlinear optical transport without altering the resonance frequencies.

4. Magnitude and Spectral Characteristics:
   The calculated peak nonlinear susceptibility reaches values on the order of $10^{-2}$ nm$^2$/V. While smaller than SHG responses in inversion-symmetry-broken BLG (e.g., via electric fields), the strain-controlled mechanism offers a distinct advantage: the ability to reversibly switch the current direction. The spectral positions of the resonances remain largely invariant under strain, while the amplitudes are strongly modulated.

Significance
The paper establishes strain engineering as an effective route for manipulating valley-dependent nonlinear optical phenomena in bilayer graphene. By demonstrating that uniaxial strain can deterministically enhance, suppress, or reverse the direction of SHG currents, the authors highlight a pathway for tunable mid-infrared photonic and valley-optoelectronic applications. The results suggest that mechanical control provides an independent tuning parameter that can be combined with electrostatic and optical methods to manipulate nonlinear optical responses in 2D materials. The study does not propose specific experimental setups but notes that the predicted effects should be detectable using modern mid-infrared nonlinear optical spectroscopy techniques.