The formalism of the reduced density matrix is pursued in both length and velocity gauges of the perturbation to the crystal Hamiltonian. The covariant derivative is introduced as a convenient representation of the position operator. This allow us to write compact expressions for the reduced density matrix in any order of the perturbation which simplifies the calculations of nonlinear optical responses; as an example, we compute the first and third order contributions of the monolayer graphene. Expressions obtained in both gauges share the same formal structure, allowing a comparison of the effects of truncation to a finite set of bands. This truncation breaks the equivalence between the two approaches: its proper implementation can be done directly in the expressions derived in the length gauge, but require a revision of the equations of motion of the reduced density matrix in the velocity gauge.

In this work, the difficulties inherent to perturbative calculations in the velocity gauge are addressed. In particular, it is shown how calculations of nonlinear optical responses in the independent particle approximation can be done to any order and for any finite band model. The procedure and advantages of the velocity gauge in such calculations are described. The addition of a phenomenological relaxation parameter is also discussed. As an illustration, the nonlinear optical response of monolayer graphene is numerically calculated using the velocity gauge. arXiv:1712.04924v2 [cond-mat.other]

In this paper we theoretically describe the absorption of hexagonal boron nitride (hBN) single layer. We develop the necessary formalism and present an efficient method for solving the Wannier equation for excitons. We give predictions for the absorption of hBN on quartz and on graphite. We compare our predictions with recently published results [Elias et al., Nat. Comm. 10, 2639] for a monolayer of hBN on graphite. The spontaneous radiative lifetime of excitons in hBN is also computed. We argue that the optical properties of hBN in the ultraviolet are very useful for the study of peptides and other biomolecules.

In this work, we present numerical results for the second and third order conductivities of the plain graphene and gapped graphene monolayers associated with the second and third harmonic generation, the optical rectification and the optical Kerr effect. The frequencies considered here range from the microwave to the ultraviolet portion of the spectrum, the latter end of which had not yet been studied. These calculations are performed in the velocity gauge and directly address the components of the conductivity tensor. In the velocity gauge, the radiation field is represented by a power series in the vector potential, and we discuss a very efficient way of calculating its coefficients in the context of tight-binding models.

The structure of the electronic nonlinear optical conductivity is elucidated in a detailed study of the time-reversal symmetric two-band model. The nonlinear conductivity is decomposed as a sum of contributions related with different regions of the first Brillouin zone, defined by single or multiphoton resonances. All contributions are written in terms of the same integrals, which contain all information specific to the particular model under study. In this way, ready-to-use formulas are provided that reduce the often tedious calculations of the second and third order optical conductivity to the evaluation of a small set of similar integrals. In the scenario where charge carriers are present prior to optical excitation, Fermi surface contributions must also be considered and are shown to have an universal frequency dependence, tunable by doping. General characteristics are made evident in this type of resonance-based analysis: the existence of step functions that determine the chemical potential dependence of electron-hole symmetric insulators; the determination of the imaginary part by Hilbert transforms, simpler than those of the nonlinear Kramers-Krönig relations; the absence of Drude peaks in the diagonal elements of the second order conductivity, among others. As examples, analytical expressions are derived for the nonlinear conductivities of some simple systems: a very basic model of direct gap semiconductors and the Dirac fermions of monolayer graphene.