Nonlinear nanophotonics has many potential applications, such as in mode locking, frequency-comb generation, and all-optical switching. The development of materials with large nonlinear susceptibility is key to realizing nonlinear nanophotonics. Nanostructured graphene systems, such as graphene nanoribbons and nanoislands, have been predicted to have a strong plasmon-enhanced nonlinear optical behavior in the nonretarded regime. Plasmons concentrate the light field down to subwavelength scales and can enhance the nonlinear optical effects; however, plasmon resonances are narrowband and sensitive to the nanostructure geometry. Here we show that graphene nanoribbons, particularly armchair graphene nanoribbons, have a remarkably strong nonlinear optical response in the long-wavelength regime and over a broad frequency range, from terahertz to the near infrared. We use a quantum-mechanical master equation with a detailed treatment of scattering and show that, in the retarded regime, electron scattering has a critical effect on the optical nonlinearity of graphene nanoribbons, which cannot be captured via the commonly used relaxation-time approximation. At terahertz frequencies, where intraband optical transitions dominate, the strong nonlinearity (in particular, third-order Kerr nonlinearity) stems from the jagged shape of the electron energy distribution, caused by the interband electron scattering mechanisms along with the intraband inelastic scattering mechanisms. We show that the relaxation-time approximation fails to capture this quantum-mechanical phenomenon and results in a significant underestimation of the intraband nonlinearity. At the midinfrared to near infrared frequencies, where interband optical transitions dominate, the Kerr nonlinearity is significantly overestimated within the relaxation-time approximation. These findings unveil the critical effect of electron scattering on the optical nonlinearity of nanostructured graphene, and also underscore the capability of this class of materials for nonlinear nanophotonic applications.

We study entanglement dynamics in a diagonal dephasing model in which the strength of interaction decays exponentially with distance—the so-called l-bit model of many-body localization. We calculate the exact expression for entanglement growth with time, finding in addition to a logarithmic growth, a sublogarithmic correction. Provided the l-bit picture correctly describes the many-body localized phase this implies that the entanglement in such systems does not grow (just) as a logarithm of time, as believed so far.

Trivial Andreev bound states arising from chemical-potential variations could lead to zero-bias tunneling conductance peaks at finite magnetic field in class-D nanowires, precisely mimicking the predicted zero-bias conductance peaks arising from the topological Majorana bound states. This finding raises a serious question on the efficacy of using zero-bias tunneling conductance peaks, by themselves, as evidence supporting the existence of topological Majorana bound states in nanowires. In the current work, we provide specific experimental protocols for tunneling spectroscopy measurements to distinguish between Andreev and Majorana bound states without invoking more demanding nonlocal measurements which have not yet been successfully performed in nanowire systems. In particular, we discuss three distinct experimental schemes involving the response of the zero-bias peak to local perturbations of the tunnel barrier, the overlap of bound states from the wire ends, and, most compellingly, introducing a sharp localized potential in the wire itself to perturb the zero-bias tunneling peaks. We provide extensive numerical simulations clarifying and supporting our theoretical predictions.

We study the attractive Hubbard model with spin imbalance on two lattices featuring a flat band: the Lieb and kagome lattices. We present mean-field phase diagrams featuring exotic superfluid phases, similar to the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, whose stability is confirmed by dynamical mean-field theory. The nature of the pairing is found to be richer than just the Fermi surface shift responsible for the usual FFLO state. The presence of a flat band allows for changes in the particle momentum distributions at null energy cost. This facilitates formation of nontrivial superfluid phases via multiband Cooper pair formation: the momentum distribution of the spin component in the flat band deforms to mimic the Fermi surface of the other spin component residing in a dispersive band. The Fermi surface of the unpaired particles that are typical for gapless superfluids becomes deformed as well. The results highlight the profound effect of flat dispersions on Fermi surface instabilities, and provide a potential route for observing spin-imbalanced superfluidity and superconductivity.

We present a finite-size scaling for both interaction and disorder strengths in the critical regime of the many-body localization (MBL) transition for a spin-1/$2XXZ$ spin chain with a random field by studying level statistics. We show how the dynamical transition from the thermal to MBL phase depends on interaction together with disorder by evaluating the ratio of adjacent level spacings, and thus, extend previous studies in which interaction coupling is fixed. We introduce an extra critical exponent in order to describe the nontrivial interaction dependence of the MBL transition. It is characterized by the ratio of the disorder strength to the power of the interaction coupling with respect to the extra critical exponent and not by the simple ratio between them.

Spin waves in chiral magnetic materials are strongly influenced by the Dzyaloshinskii-Moriya interaction, resulting in intriguing phenomena like nonreciprocal magnon propagation and magnetochiral dichroism. Here, we study the nonreciprocal magnon spectrum of the archetypical chiral magnet MnSi and its evolution as a function of magnetic field covering the field-polarized and conical helix phase. Using inelastic neutron scattering, the magnon energies and their spectral weights are determined quantitatively after deconvolution with the instrumental resolution. In the field-polarized phase the imaginary part of the dynamical susceptibility ${\chi }^{\prime \prime }(\varepsilon ,\mathbf{q})$ is shown to be asymmetric with respect to wave vectors $\mathbf{q}$ longitudinal to the applied magnetic field $\mathbf{H}$, which is a hallmark of chiral magnetism. In the helimagnetic phase, ${\chi }^{\prime \prime }(\varepsilon ,\mathbf{q})$ becomes increasingly symmetric with decreasing $\mathbf{H}$ due to the formation of helimagnon bands and the activation of additional spin-flip and non-spin-flip scattering channels. The neutron spectra are in excellent quantitative agreement with the low-energy theory of cubic chiral magnets with a single fitting parameter being the damping rate of spin waves.

We experimentally and theoretically study the geometry, as well as the electronic and vibrational properties, of the heterotetranuclear magnetic cluster ${\left[{\mathrm{Co}}_3\mathrm{Ni}\right(\mathrm{EtOH}\left)\right]}^{+}$, which is prepared in the gas phase with molecular beam expansion. We characterize the cluster and identify possible isomers through the comparison of experimentally observed infrared spectra with state-of-the-art quantum chemistry calculations, more specifically by focusing on the OH stretching frequency. Furthermore, we suggest ultrafast, laser-induced, local spin-flip scenarios on every Co atom, and report a cooperative effect, in which the spin density is localized on one Co atom, gets transiently transferred to another, and then bounces back pointing in the opposite direction. Finally, we predict a tolerance of the suggested scenarios with respect to the laser detuning of about 20 meV, which lies within an experimentally applicable range. Our joint investigation is an additional step toward the implementation of laser-controlled nanospintronic devices.

Monolayer superconducting transition-metal dichalcogenide ${\mathrm{NbSe}}_2$ is a candidate for a nodal topological superconductor by magnetic field. Because of the so-called Ising spin-orbit coupling that strongly pins the electron spins to the out-of-plane direction, Cooper pairs in monolayer superconductor ${\mathrm{NbSe}}_2$ are protected against an applied in-plane magnetic field much larger than the Pauli limit. In monolayer ${\mathrm{NbSe}}_2$, in addition to the Fermi pockets at the corners of Brillouin zone with opposite crystal momentum similar to other semiconducting transition-metal dichalcogenids, there is an extra Fermi pocket around the $\mathrm{\Gamma }$ point with much smaller spin splitting, which could lead to an alternative strategy for pairing possibilities that are manipulable by a smaller magnetic field. By considering a monolayer ${\mathrm{NbSe}}_2$-ferromagnet substrate junction, we explore the modified pairing correlations on the pocket at $\mathrm{\Gamma }$ point in hole-doped monolayer ${\mathrm{NbSe}}_2$. The underlying physics is fascinating as there is a delicate interplay of the induced exchange field and the Ising spin-orbit coupling. We realize a mixed singlet-triplet superconductivity, $s+f$, due to the Ising spin-orbit coupling. Moreover, our results reveal the admixture state including both odd- and even-frequency components, associated with the ferromagnetic proximity effect. Different frequency symmetries of the induced pairing correlations can be realized by manipulating the magnitude and direction of the induced magnetization.

The dynamical spin structure factor is computed within a variational framework to study the one-dimensional $J_1-J_2$ Heisenberg model. Starting from Gutzwiller-projected fermionic wave functions, the low-energy spectrum is constructed from two-spinon excitations. The direct comparison with Lanczos calculations on small clusters demonstrates the excellent description of both gapless and gapped (dimerized) phases, including incommensurate structures for $J_2/J_1>0.5$. Calculations on large clusters show how the intensity evolves when increasing the frustrating ratio and give an unprecedented accurate characterization of the dynamical properties of (nonintegrable) frustrated spin models.

We study the photoinduced breakdown of a two-orbital Mott insulator and resulting metallic state. Using time-dependent density matrix renormalization group, we scrutinize the real-time dynamics of the half-filled two-orbital Hubbard model interacting with a resonant radiation field pulse. The breakdown, caused by production of doublon-holon pairs, is enhanced by Hund's exchange, which dynamically activates large orbital fluctuations. The melting of the Mott insulator is accompanied by a high to low spin transition with a concomitant reduction of antiferromagnetic spin fluctuations. Most notably, the overall time response is driven by the photogeneration of excitons with orbital character that are stabilized by Hund's coupling. These unconventional “Hund excitons” correspond to bound spin-singlet orbital-triplet doublon-holon pairs. We study exciton properties such as bandwidth, binding potential, and size within a semiclassical approach. The photometallic state results from a coexistence of Hund excitons and doublon-holon plasma.

We investigate ballistic transport of helical electrons in a three-dimensional topological insulator traversing a nanostep junction. We find that a magnetic field perpendicular to its side surface shrinks the phase space for transmission, leading to magnetoresistance for the Fermi energy close to the Dirac point of the top surface. We also find transmission resonances and suppression of the Fano factor due to Landau-level-related quasibound states. The transmission blockade in the off-resonance case can result in a huge magnetoresistance for Fermi energy higher than the Dirac point of the side surface.

The effects of lattice distortion and chemical disorder on charge transport properties of two-terminal zigzag phosphorene nanoribbons (zPNRs), which shows resonant tunneling behavior under an electrical applied bias, are studied. Our comprehensive study is based on ab initio quantum transport calculations on the basis of the Landauer theory. We use nitrogen and silicon substitutional dopant atoms, and employ different physical quantities such as the $I\text{−}V$ curve, voltage drop behavior, transmission spectrum, transmission pathway, and atomic current to explore the transport mechanism of zPNR devices under a bias voltage. The calculated transmission pathways show the transition from a ballistic transport regime to a diffusive and in some particular cases to localized transport regimes. Current flowing via the chemical bonds and hopping are monitored; however, the conductance originates mainly from the charge traveling through the chemical bonds in the vicinity of the zigzag edges. Our results show that in the doped systems, the device conductance decreases and the negative differential resistance characteristic becomes weak or is eliminated. Besides, the conductance in a pure zPNR system is almost independent of the ribbon width.

We study the phase diagram of the one-dimensional boson gas trapped inside an optical lattice with contact and dipolar interaction, taking into account next-nearest terms for both tunneling and interaction. Using the density-matrix renormalization group, we calculate how the locations of phase transitions change with increasing dipolar interaction strength for average density $\rho =1$. Furthermore, we show the emergence of pair-correlated phases for a large dipolar interaction strength and $\rho \ge 2$, including a supersolid phase with an incommensurate density wave ordering manifesting the corresponding spontaneous breaking of the translational symmetry.

The newly discovered topological Dirac semimetals host the possibilities of various topological phase transitions through the control of spin-orbit coupling as well as symmetries and dimensionalities. Here, we report a magnetotransport study of high-mobility (${\mathrm{Cd}}_{1-x}{\mathrm{Zn}}_x{)}_3{\mathrm{As}}_2$ films, where the topological Dirac semimetal phase can be turned into a trivial insulator via chemical substitution. By high-field measurements with a Hall-bar geometry, magnetoresistance components ascribed to the chiral charge pumping have been distinguished from other extrinsic effects. The negative magnetoresistance exhibits a clear suppression upon Zn doping, reflecting decreasing Berry curvature of the band structure as the topological phase transition is induced by reducing the spin-orbit

A search for new topological quantum systems is challenging due to the requirement of nontrivial band connectivity that leads to protected surface states of electrons. Progress in this field was primarily due to a realization of a band inversion mechanism between even and odd parity states that was proven to be very useful in both predicting many such systems and our understanding of their topological properties. Despite many proposed materials that assume the band inversion between $s$ and $p$ (or $p/d$, $d$/$f$) electrons, here, we explore a different mechanism where the occupied $d$ states subjected to a tetrahedral crystal field produce an active $t_{2g}$ manifold behaving as a state with an effective orbital momentum equal to $-1$, and pushing $j_{\mathrm{eff}}=1/2$ doublet at a higher energy. Via hybridization with nearest-neighbor orbitals realizable, e.g., in a zinc-blende structural environment, this allows a formation of odd parity state whose subsequent band inversion with an unoccupied $s$ band becomes possible, prompting us to look for the compounds with ${\mathrm{Cu}}^{+1}$ ionic state. Chemical valence arguments coupled to a search in the materials database of zinc-blende-like lattice space groups $T_d^2$ ($F\overline{4}3m$) lead us to systematically investigate electronic structures and topological properties of CuY ($Y=\text{F}$, Cl, Br, I) and $\mathrm{Cu}X\mathrm{O}$ ($X=\text{Li}$, Na, K, Rb) families of compounds. Our theoretical results show that CuF displays a behavior characteristic of an ideal Weyl semimetal with 24 Weyl nodes at the bulk Brillouin zone. We also find that other compounds, CuNaO and CuLiO, are the $s\text{−}d$ inversion-type topological insulators. Results for their electronic structures and corresponding surfaces states are presented and discussed in the context of their topological properties.

Positronium formation at Si(111) and Si(001) surfaces has been investigated by changing the doping level systematically over the range 300–1000 K. The temperature dependence of the positronium fraction varied with the doping condition, and there were practically no differences between the two surface orientations. In heavily doped $n$-type Si $(n\gtrsim {10}^{18}{\mathrm{cm}}^{-3})$, the positronium fraction ($I_{\mathrm{Ps}}$) increased above 700 K and reached more than 95% at 1000 K. In undoped and lightly doped Si ($n$, $p\lesssim {10}^{15}{\mathrm{cm}}^{-3}$), $I_{\mathrm{Ps}}$ decreased from 300 to 500 K and increased above 700 K. In heavily doped $p$-type Si ($p\gtrsim {10}^{18}{\mathrm{cm}}^{-3}$), $I_{\mathrm{Ps}}$ increased in two steps: one at 500–600 K and one above 700 K. Overall, the positronium fraction increased with the amount of $n$-type doping. These phenomena were found to be dominated by two kinds of positronium with energies of 0.6–1.5 eV and 0.1–0.2 eV, which were attributed to the work-function mechanism and the surface-positron-mediated process, respectively, with contributions from conduction electrons. The positron work function was estimated to be positive. This agrees with first-principles calculation. The positive positron work function implies that the formation of excitonic electron-positron bound states begins in the bulk subsurface region and transits to the final positronium state in the vacuum.

Simulation of interacting electron systems is one of the great challenges of modern quantum chemistry and solid state physics. Controllable quantum systems offer the opportunity to create artificial structures which mimic the system of interest. An interesting quantity to extract from these quantum simulations is the spectral function. We map a noisy quantum simulator onto a fermionic system and investigate the influence of decoherence on the simulation of the spectral density using a diagrammatic approach on the Keldysh contour. We show that features stronger than the single-qubit decoherence rate can be resolved, while weaker features wash out. For small systems, we compare our Keldysh approach to master-equation calculations.

Synthetic antiferromagnets (SAF) have been proposed to replace ferromagnets in magnetic memory devices to reduce the stray field, increase the storage density, and improve the thermal stability. Here, we investigate the spin-orbit torque in a perpendicularly magnetized Pt/[Co/Pd]/Ru/[Co/Pd] SAF structure, which exhibits completely compensated magnetization and an exchange coupling field up to 2100 Oe. The magnetizations of two Co/Pd layers can be switched between two antiparallel states simultaneously by spin-orbit torque. The magnetization switching can be read out due to much stronger spin-orbit coupling at the bottom Pt/[Co/Pd] interface compared to its upper counterpart without Pt. Both experimental and theoretical analyses unravel that the torque efficiency of antiferromagnetically coupled stacks is significantly higher than the ferromagnetic counterpart, making the critical switching current of SAF comparable to the conventional single ferromagnet. Besides adding an important dimension to spin-orbit torque, the efficient switching of completely compensated SAF might advance magnetic memory devices with high density, high speed, and low power consumption.