The paper theoretically studies the possibility of using the effects of phonon interference between paths through different interatomic bonds for the control of phonon heat transfer through internal crystal interfaces and for the design of phonon metamirrors and meta-absorbers. These metamirrors and meta-absorbers are considered to be defect nanolayers of atomic-scale thicknesses embedded in a crystal. Several analytically solvable three-dimensional lattice-dynamics models of the phonon metamirrors and meta-absorbers at the internal crystal planes are described. It is shown that due to destructive interference in the two or more phonon paths, the internal crystal planes, fully or partially filled with weakly bound or heavy-isotope defect atoms, can completely reflect or completely absorb phonons at the transmission antiresonances, whose wavelengths are larger than the effective thickness of the metamirror or meta-absorber. Due to cooperative superradiant effect, the spectral widths of the two-path interference antiresonances for the plane waves are given by the square of partial filling fraction in the defect crystal plane. Our analysis reveals that the presence of two or more phonon paths plays the dominant role in the emergence of the transmission antiresonances in phonon scattering at the defect crystal planes and in reduction of the thermal interface conductance in comparison with the Fano-resonance concept. We study analytically phonon transmission through internal crystal plane in a model cubic lattice of Si-like atoms, partially filled with Ge-like defect atoms. Such a plane can serve as interference phonon metamirror with the transmission antiresonances in the vicinities of eigenmode frequencies of Ge-like defect atoms in the terahertz frequency range. We predict the extraordinary phonon transmission induced by the two-path constructive interference of the lattice waves in resonance with the vibrations of rare host atoms, periodically distributed in the crystal plane almost completely filled with heavy-isotope defects. We show that the phonon-interference-induced transparency can be produced by the defect nanolayer with the non-nearest-neighbor interactions, filled with two types of isotopes with relatively small difference in masses or binding force constants. In this case, relatively broad transmission antiresonance is accompanied by the narrow transmission peak close to the antiresonance frequency. We describe the softening of the flexural surface acoustic wave, localized at the embedded defect nanolayer, caused by negative surface stress in the layer. The surface wave softening results in spatially periodic static bending deformation of the embedded nanolayer with the definite wave number. The latter effect is estimated for graphene monolayer embedded in a strained matrix of polyethylene. We analyze the effect of nonlinearity in the dynamics of defect atoms on the one- and two-path phonon interference and show that the interference transmission resonances and antiresonances are shifted in frequencies but not completely suppressed by rather strong anharmonicity of interatomic bonds. The reduction of the Kapitza thermal interface conductance caused by the destructive phonon interference in a defect monolayer is described. We show that the additional relatively weak non-nearest-neighbor interactions through the defect crystal plane filled with heavy isotopes substantially reduces the interface thermal conductance, and this effect is stronger in the three-dimensional system than in the quasi-one-dimensional systems studied previously.

In this paper, we demonstrate that phonons and photons of different momenta can be confined and interact with each other within the same nanostructure. The interaction between confined phonons and confined photons in silicon resonator arrays is observed by means of Raman scattering. The Raman spectra from large arrays of dielectric silicon resonators exhibited Raman enhancement accompanied with a downshift and broadening. The analysis of the Raman intensity and line shape using finite-difference time-domain simulations and a spatial correlation model demonstrated an interaction between photons confined in the resonators and phonons confined in highly defective regions prompted by the structuring process. It was shown that the Raman enhancement is due to collective lattice resonance inducing field confinement in the resonators, while the spectra downshift and broadening are signatures of the relaxation of the phonon wave vector due to phonon confinement in defective regions located in the surface layer of the Si resonators. We found that as the resonators increase in height and their shape becomes cylindrical, the amplitude of their coherent oscillation increases and hence their ability to confine the incoming electric field increases.

When phonon transport is in the hydrodynamic regime, the thermal conductivity exhibits peculiar dependences on temperatures ($T$) and sample widths ($W$). These features were used in the past to experimentally confirm the hydrodynamic phonon transport in three-dimensional bulk materials. Suspended graphene was recently predicted to exhibit strong hydrodynamic features in thermal transport at much higher temperature than the three-dimensional bulk materials, but its experimental confirmation requires quantitative guidance by theory and simulation. Here we quantitatively predict those peculiar dependences using the Monte Carlo solution of the Peierls-Boltzmann equation with an ab initio full three-phonon scattering matrix. Thermal conductivity is found to increase as $T^{\alpha }$ where $\alpha$ ranges from 1.89 to 2.49 depending on a sample width at low temperatures, much larger than 1.68 of the ballistic case. The thermal conductivity has a width dependence of $W^{1.17}$ at 100 K, clearly distinguished from the sublinear dependence of the ballistic-diffusive regime. These peculiar features are explained with a phonon viscous damping effect of the hydrodynamic regime. We derive an expression for the phonon hydrodynamic viscosity from the Peierls-Boltzmann equation, and discuss the fact that the phonon viscous damping explains well those peculiar dependences of thermal conductivity at 100 K. The phonon viscous damping still causes significant thermal resistance when a temperature is 300 K and a sample width is around 1 µm, even though the hydrodynamic regime is not dominant over other regimes at this condition.

We study the effective spin-orbital model for honeycomb-layered transition metal compounds, applying the second-order perturbation theory to the three-orbital Hubbard model with the anisotropic hoppings. This model is reduced to the Kitaev model in the strong spin-orbit coupling limit. Combining the cluster mean-field approximations with the exact diagonalization, we treat the Kugel-Khomskii type superexchange interaction and spin-orbit coupling on an equal footing to discuss ground-state properties. We find that a zigzag ordered state is realized in the model within nearest-neighbor interactions. We clarify how the ordered state competes with the nonmagnetic state, which is adiabatically connected to the quantum spin liquid state realized in a strong spin-orbit coupling limit. Thermodynamic properties are also addressed. The present paper should provide another route to account for the Kitaev-based magnetic properties in candidate materials.

In multicomponent superfluids and superconductors, co- and counterflows of components have, in general, different properties. A. F. Andreev and E. P. Bashkin [Sov. Phys. JETP 42, 164 (1975)] discussed, in the context of ${\mathrm{He}}^3/{\mathrm{He}}^4$ superfluid mixtures, that interparticle interactions produce a dissipationless drag. The drag can be understood as a superflow of one component induced by phase gradients of the other component. Importantly, the drag can be both positive (entrainment) and negative (counterflow). The effect is known to have crucial importance for many properties of diverse physical systems ranging from the dynamics of neutron stars and rotational responses of Bose mixtures of ultracold atoms to magnetic responses of multicomponent superconductors. Although substantial literature exists that includes the drag interaction phenomenologically, only a few regimes are covered by quantitative studies of the microscopic origin of the drag and its dependence on microscopic parameters. Here we study the microscopic origin and strength of the drag interaction in a quantum system of two-component bosons on a lattice with short-range interaction. By performing quantum Monte Carlo simulations of a two-component Bose-Hubbard model we obtain dependencies of the drag strength on the boson-boson interactions and properties of the optical lattice. Of particular interest are the strongly correlated regimes where the ratio of coflow and counterflow superfluid stiffnesses can diverge, corresponding to the case of saturated drag.

Topological metamaterials exhibit unusual behaviors at their boundaries, such as unidirectional chiral waves, that are protected by a topological feature of their band structures. The ability to tune such a material through a topological phase transition in real time could enable the use of protected waves for information storage and readout. Here we dynamically tune through a topological phase transition by breaking inversion symmetry in a metamaterial composed of interacting gyroscopes. Through the transition, we track the divergence of the edge modes' localization length and the change in Chern number characterizing the topology of the material's band structure. This Rapid Communication provides a new axis with which to tune the response of mechanical topological metamaterials.

An atomistic effective Hamiltonian is used to investigate electrocaloric (EC) effects of $\mathrm{Pb}\left({\mathrm{Mg}}_{1/3}{\mathrm{Nb}}_{2/3}\right){\mathrm{O}}_3$ relaxor ferroelectrics in its ergodic regime, and subject to electric fields applied along the pseudocubic [111] direction. Such a Hamiltonian qualitatively reproduces (i) the electric field-versus-temperature phase diagram, including the existence of a critical point where first-order and second-order transitions meet each other; and (ii) a giant EC response near such a critical point. It also reveals that such giant response around this critical point is microscopically induced by field-induced percolation of polar nanoregions. Moreover, it is also found that, for any temperature above the critical point, the EC coefficient-versus-electric-field curve adopts a maximum (and thus larger electrocaloric response too), that can be well described by the general Landau-like model proposed by Jiang et al., [Phys. Rev. B 96, 014114 (2017)], and that is further correlated with specific microscopic features related to dipoles lying along different rhombohedral directions. Furthermore, for temperatures being at least 40 K higher than the critical temperature, the (electric field, temperature) line associated with this maximal EC coefficient is below both the Widom line and the line representing percolation of polar nanoregions.

The perovskite ${\mathrm{BaSnO}}_3$ provides a promising platform for the realization of an earth-abundant $n$-type transparent conductor. Its optical properties are dominated by a dispersive conduction band of Sn $5s$ states and by a flatter valence band of O $2p$ states, with an overall indirect gap of about $2.9\mathrm{eV}$. Using first-principles methods, we study the optical properties of ${\mathrm{BaSnO}}_3$ and show that both electron-phonon interactions and exact exchange, included using a hybrid functional, are necessary to obtain a qualitatively correct description of optical absorption in this material. In particular, the electron-phonon interaction drives phonon-assisted optical absorption across the minimum indirect gap and therefore determines the absorption onset, and it also leads to the temperature dependence of the absorption spectrum. Electronic correlations beyond semilocal density functional theory are key to determine the dynamical stability of the cubic perovskite structure, as well as the correct energies of the conduction bands that dominate absorption. Our work demonstrates that phonon-mediated absorption processes should be included in the design of novel transparent conductor materials.

We report the magnetic excitation spectrum as measured by inelastic neutron scattering for a polycrystalline sample of ${\mathrm{Sr}}_3{\mathrm{CuPtO}}_6$. Modeling the data by the 2+4 spinon contributions to the dynamical susceptibility within the chains, and with interchain coupling treated in the random phase approximation, accounts for the major features of the powder-averaged structure factor. The magnetic excitations broaden considerably as temperature is raised, persisting up to above 100 K and displaying a broad transition as previously seen in the susceptibility data. No spin gap is observed in the dispersive spin excitations at low momentum transfer, which is consistent with the gapless spinon continuum expected from the coordinate Bethe ansatz. However, the temperature dependence of the excitation spectrum gives evidence of some very weak interchain coupling.

We developed an efficient method for Monte Carlo simulations with long-range interacting systems. The present method is based on the stochastic potential switching algorithm as in the case of the stochastic cutoff (SCO) method. Compared with the SCO method, the present method is efficient for systems with complicated lattice structures, e.g., amorphous systems and translationally invariant systems with large unit cells. Using the present method, we investigated magnetic structures of ${\mathrm{Nd}}_2{\mathrm{Fe}}_{14}\mathrm{B}$ thin-film magnets based on the atomistic spin model, which contains 68 continuous spins in the unit cell. We found that the anisotropy of Nd ions contributes to that of the whole system only at relatively low temperatures. If the anisotropy of the Fe ions is not considered, a planer magnetic structure arises owing to the dipole-dipole interactions at high temperatures, while it does not with the anisotropy. This difference indicates that the anisotropy of the Fe ions plays an important role in holding the magnetization parallel to the easy axis at the high temperatures.

Stripe domains are narrow, elongated, reversed regions that exist in magnetic materials with perpendicular magnetic anisotropy. They appear as a pair of domain walls that can exhibit topology with a nonzero chirality. Recent experimental and numerical investigations identify an instability of stripe domains along the long direction as a means of nucleating isolated magnetic skyrmions. Here, the onset and nonlinear evolution of transverse instabilities for a dynamic stripe domain known as the bion stripe are investigated. Both nontopological and topological variants of the bion stripe are shown to exhibit a long-wavelength transverse instability with different characteristic features. In the former, small transverse variations in the stripe's width lead to a neck instability that eventually pinches the nontopological stripe into a chain of two-dimensional breathers composed of droplet soliton pairs. In the latter case, small variations in the stripe's center result in a snake instability whose topological structure leads to the nucleation of dynamic magnetic skyrmions and antiskyrmions as well as perimeter-modulated droplets. Quantitative, analytical predictions for both the early, linear evolution and the long-time, nonlinear evolution are achieved using an averaged Lagrangian approach that incorporates both exchange (dispersion) and anisotropy (nonlinearity). The method of analysis is general and can be applied to other filamentary structures.

We investigate topological magnon band crossings of pyrochlore antiferromagnets with all-in-all-out (AIAO) magnetic order. By general symmetry analysis and spin-wave theory, we show that pyrochlore materials with AIAO orders can host Weyl magnons under external magnetic fields or uniaxial strains. Under a small magnetic field, the magnon bands of the pyrochlore with AIAO background can feature two opposite-charged Weyl points, which is the minimal number of Weyl points realizable in quantum materials, and has not been experimentally observed so far. We further show that breathing pyrochlores with AIAO orders can exhibit Weyl magnons upon uniaxial strains. These findings apply to any pyrochlore material supporting AIAO orders, irrespective of the forms of interactions. Specifically, we show that the Weyl magnons are robust against direct (positive) Dzyaloshinskii-Moriya interactions. Because of the ubiquitous AIAO orders in pyrochlore magnets including ${\mathrm{R}}_2{\mathrm{Ir}}_2{\mathrm{O}}_7$, and experimentally achievable external strain and magnetic field, our predictions provide a promising arena to witness the Weyl magnons in quantum magnets.

$\beta -\mathrm{G}{\mathrm{a}}_2{\mathrm{O}}_3$, a candidate material for power electronics and UV optoelectronics, shows strong room-temperature photoluminescence (PL). In addition to the three well-known bands of as-grown samples in the UV, blue, and green, also red PL was observed upon nitrogen doping. This raises the possibility of applying $\beta -\mathrm{G}{\mathrm{a}}_2{\mathrm{O}}_3$ nanostructures as white phosphors. Using an optimized, Koopmans-compliant hybrid functional, we show that most intrinsic point defects, as well as substitutional nitrogen, act as deep acceptors, and each of the observed PL bands can be explained by electron recombination with a hole trapped in one of them. We suggest this mechanism to be general in wide-band-gap semiconductors which can only be doped $n$-type. Calculations on the nitrogen acceptor reproduce the observed red luminescence accurately. Earlier we have shown that not only the energy, but the polarization properties of the UV band can be explained by self-trapped hole states. Here we find that the blue band has its origin mainly in singly negative Ga-O divacancies, and the green band is caused dominantly by interstitial O atoms (with minor contribution of Ga vacancies to both). These assignments can explain the experimentally observed dependence of the PL bands on free-electron concentration and stoichiometry. The information provided here paves the way for the conscious tuning of light emission from $\beta -\mathrm{G}{\mathrm{a}}_2{\mathrm{O}}_3$.

We present results for many-body perturbation theory for the one-body Green's function at finite temperatures using the Matsubara formalism. Our method relies on the accurate representation of the single-particle states in standard Gaussian basis sets, allowing to efficiently compute, among other observables, quasiparticle energies and Dyson orbitals of atoms and molecules. In particular, we challenge the second-order treatment of the Coulomb interaction by benchmarking its accuracy for a well-established test set of small molecules, which includes also systems where the usual Hartree-Fock treatment encounters difficulties. We discuss different schemes how to extract quasiparticle properties and assess their range of applicability. With an accurate solution and compact representation, our method is an ideal starting point to study electron dynamics in time-resolved experiments by the propagation of the Kadanoff-Baym equations.

We report a detailed theoretical study of the electronic structure, spectral properties, and lattice parameters of bulk FeSe under pressure using a fully charge self-consistent implementation of the density functional theory plus dynamical mean-field theory method (DFT+DMFT). In particular, we perform a structural optimization and compute the evolution of the lattice parameters (volume, $c/a$ ratio, and the internal $z$ position of Se) and the electronic structure of the tetragonal (space group $P4/nmm$) unit cell of paramagnetic FeSe. Our results for the lattice parameters obtained by structural optimization using DFT+DMFT are in good quantitative agreement with experiment, implying a crucial importance of electron correlations in determining the correct lattice properties of FeSe. Most importantly, upon compression to 10 GPa our results reveal a topological change in the Fermi surface (Lifshitz transition) which is accompanied by a two- to three-dimensional crossover and a small reduction of the quasiparticle mass renormalization compared to ambient pressure. The behavior of the momentum-resolved magnetic susceptibility $\chi \left(\mathbf{q}\right)$ shows no topological changes of magnetic correlations under pressure but demonstrates a reduction of the degree of the in-plane $(\pi ,\pi )$ stripe-type nesting. Our results for the electronic structure and lattice parameters of FeSe are in good qualitative agreement with recent experiments on its isoelectronic counterpart ${\mathrm{FeSe}}_{1-x}{\mathrm{S}}_x$.

We have synthesized single crystals of Dirac semimetal candidates $A{\mathrm{ZnBi}}_2$ with $A=\text{Ba}$ and Sr. In contrast to $A=\text{Sr}$, the Ba material displays a local Zn vacancy ordering, which makes the observation of quantum oscillations in out-of-plane magnetic fields possible. As a Dirac semimetal candidate, ${\mathrm{BaZnBi}}_2$ exhibits a small cyclotron electron mass, high quantum mobility, and nontrivial Berry phases. Three Dirac dispersions are observed by angle-resolved photoemission spectroscopy and identified by first-principles band-structure calculations. Compared to $A\text{Mn}{(\text{Bi}/\mathrm{Sb})}_2$ systems which host Mn magnetic moments, ${\mathrm{BaZnBi}}_2$ acts as a nonmagnetic analog to investigate the intrinsic properties of Dirac fermions in this structure family.

Spectra of magnetoplasma excitations have been investigated in two-dimensional electron systems in AlAs quantum wells (QWs) of different widths. The magnetoplasma spectrum has been found to change profoundly when the quantum well width becomes thinner than 5.5 nm, indicating a drastic change in the conduction electron energy spectrum. The transformation can be interpreted in terms of transition from the in-plane strongly anisotropic $X_x-X_y$ valley occupation to the out-of-plane isotropic $X_z$ valley in the QW plane. Strong enhancement of the cyclotron effective mass over the band value in narrow AlAs QWs is reported.

Generation of high harmonics in a monolayer graphene initiated by a strong coherent radiation field, taking into account electron-electron Coulomb interaction, is investigated. A microscopic theory describing the nonlinear optical response of graphene is developed. The Coulomb interaction of electrons is treated in the scope of dynamic Hartree-Fock approximation. The closed set of integrodifferential equations for the single-particle density matrix of a graphene quantum structure is solved numerically. The obtained solutions show the significance of many-body Coulomb interaction on the high harmonic generation process in graphene.

By means of a nonlinear two-component hydrodynamic model, we study the valley-polarized collective motion of electrons in a strained graphene sheet. The self-consistent numerical solution in real space indicates the existence of valley-polarized edge plasmons due to a strain-induced pseudomagnetic field. The valley polarization of the edge pseudomagnetoplasmon can occur in a specific valley, depending on the pseudomagnetic field and the electron density in equilibrium. A full valley polarization is achieved at the edge of the graphene sheet for a pseudomagnetic field of tens of Tesla, which is a realistic value in current experimental technologies.

A repulsive Coulomb interaction between electrons in different orbitals in correlated materials can give rise to bound quasiparticle states. We study the nonhybridized two-orbital Hubbard model with intra- (inter)orbital interaction $U\left(U_{12}\right)$ and different bandwidths using an improved dynamical mean-field theory numerical technique which leads to reliable spectra on the real energy axis directly at zero temperature. We find that a finite density of states at the Fermi energy in one band is correlated with the emergence of well-defined quasiparticle states at excited energies $\mathrm{\Delta }=U-U_{12}$ in the other band. These excitations are interband holon-doublon bound states. At the symmetric point $U=U_{12}$, the quasiparticle peaks are located at the Fermi energy, leading to a simultaneous and continuous Mott transition settling a long-standing controversy.