Date: Thu, 9 May 2024 | Total: 61
Phonon polaritons are hybrid states of light and matter that are typically realised when optically active phonons couple strongly to photons. We suggest a new approach to realising phonon polaritons, by employing a transverse-pumping Raman scheme, as used in experiments on cold atoms in optical cavities. This approach allows hybridisation between an optical cavity mode and any Raman-active phonon mode. Moreover, this approach enables one to tune the effective phonon-photon coupling by changing the strength of the transverse pumping light. We show that such a system may realise a phonon-polariton condensate. To do this, we find the stationary states and use Floquet theory to determine their stability. We thus identify distinct superradiant and lasing states in which the polariton modes are macroscopically populated. We map out the phase diagram of these states as a function of pump frequencies and strengths. Using parameters for transition metal dichalcogenides, we show that realisation of these phases may be practicably obtainable. The ability to manipulate phonon mode frequencies and attain steady-state populations of selected phonon modes provides a new tool for engineering correlated states of electrons.
Chains and arrays of phosphorus donors in silicon have recently been used to demonstrate dopant-based quantum simulators. The dopant disorder present in fabricated devices must be accounted for. Here, we theoretically study transport through disordered donor-based $3\times 3$ arrays that model recent experimental results. We employ a theory that combines the exact diagonalization of an extended Hubbard model of the array with a non-equilibrium Green's function formalism to model transport in interacting systems. We show that current flow through the array and features of measured stability diagrams are highly resilient to disorder. We interpret this as an emergence of uncomplicated behavior in the multi-electron system dominated by strong correlations, regardless of array filling, where the current follows the shortest paths between source and drain sites that avoid possible obstacles. The reference $3\times 3$ array has transport properties very similar to three parallel 3-site chains coupled only by interchain Coulomb interaction, which indicates a challenge in characterizing such devices.
We investigate a disordered multi-dimensional linear system in which the interaction parameters vary stochastically in time with defined temporal correlations. We refer to this type of disorder as "annealed", in contrast to quenched disorder in which couplings are fixed in time. We extend Dynamical Mean-Field Theory to accommodate annealed disorder and employ it to find the exact solution of the linear model in the limit of a large number of degrees of freedom. Our analysis yields analytical results for the non-stationary auto-correlation, the stationary variance, the power spectral density, and the phase diagram of the model. Interestingly, some unexpected features emerge upon changing the correlation time of the interactions. The stationary variance of the system and the critical variance of the disorder are generally found to be a non-monotonic function of the correlation time of the interactions. We also find that in some cases a re-entrant phase transition takes place when this correlation time is varied.
Quantum technologies would benefit from the development of high performance quantum defects acting as single-photon emitters or spin-photon interface. Finding such a quantum defect in silicon is especially appealing in view of its favorable spin bath and high processability. While some color centers in silicon have been emerging in quantum applications, there is still a need to search and develop new high performance quantum emitters. Searching a high-throughput computational database of more than 22,000 charged complex defects in silicon, we identify a series of defects formed by a group III element combined with carbon ((A-C)$\rm _{Si}$ with A=B,Al,Ga,In,Tl) and substituting on a silicon site. These defects are analogous structurally, electronically and chemically to the well-known T center in silicon ((C-C-H)$\rm_{Si}$) and their optical properties are mainly driven by an unpaired electron in a carbon $p$ orbital. They all emit in the telecom and some of these color centers show improved properties compared to the T center in terms of computed radiative lifetime or emission efficiency. We also show that the synthesis of hydrogenated T center-like defects followed by a dehydrogenation annealing step could be an efficient way of synthesis. All the T center-like defects show a higher symmetry than the T center making them easier to align with magnetic fields. Our work motivates further studies on the synthesis and control of this new family of quantum defects, and also demonstrates the use of high-throughput computational screening to detect new complex quantum defects.
Efficient, reliable and easy-to-use structure recognition of atomic environments is essential for the analysis of atomic scale computer simulations. In this work, we train two neuronal network (NN) architectures, namely PointNet and dynamic graph convolutional NN (DG-CNN) using different hyperparameters and training regimes to assess their performance in structure identification tasks of atomistic structure data. We show benchmarks on simple crystal structures, where we can compare against established methods. The approach is subsequently extended to structurally more complex SiO$_2$ phases. By making use of this structure recognition tool, we are able to achieve a deeper understanding of the crystallization process in amorphous SiO$_2$ under shock compression. Lastly, we show how the NN based structure identification workflows can be integrated into OVITO using its python interface.
Motivated by the recent experimental study on a quantum Ising magnet $\text{K}_2\text{Co}(\text{SeO}_3)_2$ where spectroscopic evidence of zero-field supersolidity is presented [arXiv: 2402.15869], we simulate the excitation spectrum of the corresponding microscopic $XXZ$ model for the compound, using the recently developed excitation ansatz of infinite projected entangled-pair states (iPEPS). We map out the ground state phase diagram and compute the dynamical spin structure factors across a range of magnetic field strengths, focusing especially on the two supersolid phases found near zero and saturation fields. Our simulated excitation spectra for the zero-field supersolid "Y" phase are in excellent agreement with the experimental data -- recovering the low-energy branches and integer quantized excited energy levels $\omega_n=nJ_{zz}$. Furthermore, we demonstrate the nonlocal multi-spin-flip features for modes at $\omega_2$, indicative of their multi-magnon nature. Additionally, we identify characteristics of the high-field supersolid "V" phase in the simulated spectra, to be compared with future experimental results.
We consider a flocking system confined transversally between two infinite reflecting parallel walls separated by a distance $L_\perp$. Infinite or periodic boundary conditions are assumed longitudinally to the direction of collective motion, defining a ring geometry typical of experimental realizations with flocking active colloids. Such a confinement selects a flocking state with its mean direction aligned parallel to the wall, thus breaking explicitly the rotational symmetry locally by a boundary effect. Finite size scaling analysis and numerical simulations show that confinement induces an effective mass term ${M_c} \sim L_\perp^{-\zeta}$ (with positive $\zeta$ being the dynamical scaling exponent of the free theory) suppressing scale free correlations at small wave-numbers. However, due to the finite system size in the transversal direction, this effect can only be detected for large enough longitudinal system sizes (i.e. narrow ring geometries). Furthermore, in the longitudinal direction, density correlations are characterized by an anomalous effective mass term. The effective mass term also enhances the global scalar order parameter and suppresses fluctuations of the mean flocking direction. These results suggest an equivalence between transversal confinement and driving by an homogeneous external field, which breaks the rotational symmetry at the global level.
The integer quantized Hall effect reveals a state of scattering-free carrier transport and quantized resistance in a two-dimensional conductor exposed to a perpendicular magnetic field. The quantized resistance is observed for the Hall bar geometry, i.e., if the current carrying contacts are connected by sample edges. A widely accepted model is the Landauer-Buttiker picture, which assumes an incompressible, i.e., electrically insulating bulk state surrounded by one-dimensional edge channels giving rise to quantized resistance. This model is challenged by the screening theory, which takes into account electron-electron interaction and predicts for the quantized Hall plateaus current flow in incompressible strips, which gradually shifts from the sample edges into the bulk with increasing magnetic field. We present direct proof of the predicted scattering-free bulk transport by exploring a Hall bar augmented with an additional contact placed in its center away from the Hall bar edges. Our result supports the screening theory.
A finite cycle time information engine based on a two-level system in contact with a thermal reservoir is studied analytically. The model for the engine incorporates an error in measuring the system's state and time delay between the measurement and the feedback process. The efficiency and power of the engine in steady state are derived as a function of level spacing, feedback delay time, engine cycle time, and measurement error. For a fixed value of level spacing and feedback delay, there is an upper bound on measurement error such that the engine can extract positive work. This threshold value of error is found to be independent of the cycle time. For a range of values of level spacing and feedback delay time, efficiency has a non-monotonic dependence on the measurement error, implying that there is an optimal measurement error for the information engine to operate efficiently. At high temperatures and with precise measurement, the engine's ability to extract positive work is extended over a larger range of feedback delay time.
We provide a method for calculating anharmonic lattice dynamics, by building a surrogate model based on Gaussian Processes (GPs). Due to the underlying Gaussian form of a GP, the model is infinitely differentiable. This allows us to train the model trained directly on forces (the derivative of PESs) reducing the evaluations required for a given accuracy. We can extend this differentiation to directly calculate second and third order force-constants using automatic differentiation (AD). For the five model materials we study, we find that the force-constants are in close agreement with a standard finite-displacement approach. Our method appears to be linear scaling in the number of atoms at predicting both second and third-order (anharmonic) force-constants.
We report low-temperature powder X-ray and neutron diffraction studies on breathing pyrochlore magnets Cu$M$Cr$_{4}$S$_{8}$ ($M$ = Al, Ga), which undergo a magnetic transition at $T_{\rm N} \approx$ 21 and 31 K for {$M$ = Al and Ga, respectively. X-ray diffraction reveals that the magnetic transition accompanies a structural transition from cubic $F{\overline 4}3m$ to polar orthorhombic $Imm2$ symmetry for both the compounds, with larger distortion observed for $M$ = Ga at low temperatures. Neutron diffraction reveals incommensurate magnetic modulation ${\mathbf Q} = (q_{\rm IC}, 0.5, 0)$ in the orthorhombic setting, where $q_{\rm IC} \approx$ 0.39 and 0.31 for $M$ = Al and Ga, respectively. Our magnetic-structure analysis suggests cycloid-type magnetic order but not proper-screw type for both the compounds. We find strong correlation between the local spin configuration and Cr-Cr bond lengths, indicating that the spin-lattice coupling as well as the magnetic frustration play an important role in determining the ground state. Cu$M$Cr$_{4}$S$_{8}$ potentially offers a platform to explore magnetoelectric effects arising from the helimagnet driven electric polarity.
We analyze the pinning of vortices for a stationary rotating dipolar supersolid along the low-density paths between droplets as a function of the rotation frequency. We restrict ourselves to the stationary configurations of vortices with the same symmetry as that of the array of droplets. Our approach exploits the fact that the wave function of each droplet acquires a linear phase on the coordinates, and hence the relative phases between neighboring droplets allows us to predict the position of the vortices. For a confined system, the estimate accurately reproduces the Gross-Pitaevskii results in the spatial regions where the neighboring droplets are well defined.
We show that magnetic insulators with a collinear and compensated order can exhibit a thermal Hall effect even at zero magnetic field if they have altermagnetic symmetries. We predict a finite thermal Hall conductivity vector $\boldsymbol{\kappa}_\text{H}$ for a rutile-inspired effective spin model with Dzyaloshinskii-Moriya interaction. Within the linear spin-wave theory, we identify two magnon branches that carry identical Berry curvature and give rise to a finite $\boldsymbol{\kappa}_\text{H}$, which can be controlled by the N\'eel vector orientation and by strain. The thermal Hall response is further complemented with a spin Nernst response to contrast spin and heat transport in altermagnetic insulators with those in ferromagnets and antiferromagnets. Our results establish the crystal thermal Hall effect of magnons and we discuss material candidates for experimental realization, such as MnF$_2$, CoF$_2$, and NiF$_2$.
In recent years, simulation methods based on the scaling of atomic potential functions, such as quasi-coarse-grained dynamics and coarse-grained dynamics, have shown promising results for modeling crystalline systems at multiple scales. However, this letter presents evidence suggesting that the spatiotemporal trajectories of coarse-grained systems generated by such simulation methods exhibit a complete correspondence with those of specific molecular dynamics systems. In essence, current coarse-grained simulation methods involve a direct amplification of the results obtained from molecular dynamics simulations across spatial and temporal scales, yet they may lack the capability to adequately capture authentic scale effects. Consequently, the findings of related studies warrant careful re-evaluation. Furthermore, this study underscores the importance of not only verifying the consistency of mesoscale simulation methods with microscopic simulations but also meticulously assessing their capability to accurately forecast mesoscale physical phenomena.
We find an unconventional giant negative exchange bias (EB) of $H_{\rm EB}$ = --14.1~kOe at 2~K (cooling field of 50~kOe) in the cluster spin-glass (CSG) La$_{0.5}$Sr$_{0.5}$Co$_{0.85}$Nb$_{0.15}$O$_3$ perovsikte cobaltites. The magnetic memory effect, aging measurements, and nonlineraity in specific heat capacity reveal the glassy magnetic state at low temperatures. Further, the detailed analysis of {\it ac-}magnetic susceptibility confirms the glassy state below $\sim$58~K and the obtained characteristic spin-relaxation time-scale of $\tau_0$ = 8.4$\times$10$^{-10}$~s indicates the presence of CSG. Moreover, the analysis of magnetic training effect using the classical EB relaxation model reveals that the frozen spins relax slowly as compared to the rotatable spins at the interface of antiferromagnetic/ferromagnetic (AFM/FM) regions in CSG. Interestingly, the dependence of EB parameters is found to be unconventional for cooling field $>$50~kOe as the $H\rm_{EB}$ and $M\rm_{EB}$ show decreasing trend instead of expected saturation at higher fields. This unusual nature emerges due to large negative values of intrinsic interface exchange coupling ($J_i$), i.e., --10.24$\pm$0.22~meV and --12.55$\pm$0.49~meV for the measuring fields of $\pm$50~kOe and $\pm$90~kOe, respectively, whereas the number of spins in the FM cluster ($N_{\rm FM}$) are found to be small in the range of 2.4--3.1. These obtained values of $J_i $ and $N_{\rm FM}$ indicate the dominant AFM interactions and the presence of FM clusters in the AFM matrix, respectively, which correlate well with the observed unconventional behavior of giant negative exchange bias in the present sample.
We theoretically investigate spin pumping into a quantum easy-plane ferromagnetic spin chain system. This quantum spin chain is effectively described by the Tomonaga-Luttinger (TL) liquid despite the ferromagnetic exchange interaction because of the easy-plane magnetic anisotropy. This TL liquid state has an extremely strong interaction that is hardly realized in other quantum antiferromagnetic chain systems or weakly interacting electron systems. We show how the strongly interacting TL liquid affects the ferromagnetic resonance that occurs in the ferromagnetic insulator. In particular, we discuss the dependence of the Gilbert damping on the temperature and the junction length. The Gilbert damping allows us to extract information about the above-mentioned strong interaction within the quantum ferromagnetic spin chain. We also point out that a well-known compound CsCuCl$_3$ will be suitable for the realization of our setup.
We consider a mixture of bosons and spin-polarized fermions in two dimensions at zero temperature with a tunable Bose-Fermi attraction. By adopting a diagrammatic T-matrix approach, we analyze the behavior of several thermodynamic quantities for the two species as a function of the density ratio and coupling strength, including the chemical potentials, the momentum distribution functions, the boson condensate density, and the Tan's contact parameter. By increasing the Bose-Fermi attraction, we find that the condensate is progressively depleted and Bose-Fermi pairs form, with a small fraction of condensed bosons surviving even for strong Bose-Fermi attraction. This small condensate proves sufficient to hybridize molecular and atomic states, producing quasi-particles with unusual Fermi liquid features. A nearly universal behavior of the condensate fraction, the bosonic momentum distribution, and Tan's contact parameter with respect to the density ratio is also found.
We perform first-principles density functional calculations to elucidate structural, electronic and magnetic properties at the interface of \ce{(0001)Cr2O3-(111)Pt} bilayers. This investigation is motivated by the fact that, despite the promise of \ce{Cr2O3-Pt} heterostructures in a variety of antiferromagnetic spintronic applications, many key structural, electronic, and magnetic properties at the \ce{Cr2O3-Pt} interface are poorly understood. We first analyze all inequivalent lateral interface alignments to determine the lowest energy interfacial structure. For all lateral alignments including the lowest-energy one, we observe an accumulation of electrons at the interface between \ce{Cr2O3} and Pt. We find an unexpected reversal of the magnetic moments of the interface Cr ions in the presence of Pt compared to surface Cr moments in vacuum-terminated \ce{(0001)Cr2O3}. We also find that the heterostructure exhibits a magnetic proximity effect in the first three Pt layers at the interface with \ce{Cr2O3}, providing a mechanism by which the anomalous Hall effect can occur in \ce{(0001)Cr2O3-(111)Pt} bilayers. Our results provide the basis for a more nuanced interpretation of magnetotransport experiments on \ce{(0001)Cr2O3-(111)Pt} bilayers and should inform future development of improved antiferromagnetic spintronic devices based on the \ce{Cr2O3-Pt} material system.
Stochastic, spatially extended models for predator-prey interaction display spatio-temporal structures that are not captured by the Lotka-Volterra mean-field rate equations. These spreading activity fronts reflect persistent correlations between predators and prey that can be analyzed through field-theoretic methods. Introducing local restrictions on the prey population induces a predator extinction threshold, with the critical dynamics at this continuous active-to-absorbing state transition governed by the scaling exponents of directed percolation. Novel features in biologically motivated model variants include the stabilizing effect of a periodically varying carrying capacity that describes seasonally oscillating resource availability; enhanced mean species densities and local fluctuations caused by spatially varying reaction rates; and intriguing evolutionary dynamics emerging when variable interaction rates are affixed to individuals combined with trait inheritance to their offspring. The basic susceptible-infected-susceptible and susceptible-infected-recovered models for infectious disease spreading near their epidemic thresholds are respectively captured by the directed and dynamic isotropic percolation universality classes. Systems with three cyclically competing species akin to spatial rock-paper-scissors games may display striking spiral patterns, yet conservation laws can prevent such noise-induced structure formation. In diffusively coupled inhomogeneous settings, one may observe the stabilization of vulnerable ecologies prone to finite-size extinction or fixation due to immigration waves emanating from the interfaces.
Magnonic black holes-analogue event horizons for the spin-wave collective excitations of ordered magnets-can be used for fundamental research, for example for investigating Hawking radiation, but also for technological applications of spin waves. Here we show how to engineer a horizon for spin waves in antiferromagnets, which have the attractive feature of fast magnetization dynamics and linear dispersion relation. We propose a set-up with spatially varying exchange interaction with spin transfer torque to implement the horizon and a second set-up for the amplification of spin waves consisting of an antiferromagnet subject to a spatially varying external magnetic field that is driven by spin orbit torque. We compute the values of parameters needed to implement the horizon and to have amplification of spin waves. We develop the corresponding Klein-Gordon equation and quantify the amplification. Our work paves the way for investigation of Hawking radiation of spin waves and for antiferromagnet-based spin-waves amplifiers.
Elastoplastic properties of nanocrystalline metals are non-uniform on the scale of the grain size, and this non-uniformity affects macroscopic quantities as, in these systems, a significant part of the material is at or adjacent to a grain boundary. We use molecular dynamics simulations to study the spatial distributions of local elastic moduli in nano-grained pure metals and analyze their dependence on grain size. Calculations are performed for copper and tantalum with grain sizes ranging from 5-20nm. Shear modulus distributions for grain and grain-boundary atoms were calculated. It is shown that the non-crystalline grain boundary has a wide shear-modulus distribution, which is grain-size independent, while grains have a peaked distribution, which becomes sharper with increasing grain size. Average elastic moduli of the bulk, grains, and grain boundary are calculated as a function of grain size. The atomistic simulations show that the reduction of total elastic moduli with decreasing grain size is mainly due to a resulting larger grain-boundary atoms fraction, and that the total elastic moduli can be approximated by a simple weighted average of larger grains elastic moduli and a lower grain-boundary elastic moduli.
Accurate and fast prediction of materials properties is central to the digital transformation of materials design. However, the vast design space and diverse operating conditions pose significant challenges for accurately modeling arbitrary material candidates and forecasting their properties. We present MatterSim, a deep learning model actively learned from large-scale first-principles computations, for efficient atomistic simulations at first-principles level and accurate prediction of broad material properties across the periodic table, spanning temperatures from 0 to 5000 K and pressures up to 1000 GPa. Out-of-the-box, the model serves as a machine learning force field, and shows remarkable capabilities not only in predicting ground-state material structures and energetics, but also in simulating their behavior under realistic temperatures and pressures, signifying an up to ten-fold enhancement in precision compared to the prior best-in-class. This enables MatterSim to compute materials' lattice dynamics, mechanical and thermodynamic properties, and beyond, to an accuracy comparable with first-principles methods. Specifically, MatterSim predicts Gibbs free energies for a wide range of inorganic solids with near-first-principles accuracy and achieves a 15 meV/atom resolution for temperatures up to 1000K compared with experiments. This opens an opportunity to predict experimental phase diagrams of materials at minimal computational cost. Moreover, MatterSim also serves as a platform for continuous learning and customization by integrating domain-specific data. The model can be fine-tuned for atomistic simulations at a desired level of theory or for direct structure-to-property predictions, achieving high data efficiency with a reduction in data requirements by up to 97%.
We propose a modification to the study of site-wise dynamically disordered totally asymmetric simple exclusion process (TASEP). Motivated by the process of gene transcription, a study in ref. [39] introduced an extension of TASEP, where the defects (or obstacles) bind/un-bind dynamically to the sites of the lattice and the hopping of the particles on lattice faces a hindrance if the arrival site is occupied by an obstacle. In addition, the particle is only allowed to enter the lattice provided the first site is defect-free. In our study, we propose that the particle movement at the entry of the lattice must face an equal hindrance that is provided by the obstacles to the rest of the particles on the lattice. For open boundaries, the continuum mean-field equations are derived and solved numerically to obtain steady-state phase diagrams and density profiles. The presence of obstacles produces a shift in the phase boundaries obtained but the same three phases as obtained for the standard TASEP. Contrary to the model introduced in ref. \cite{waclaw2019totally}, the idea to introduce the modification at the entrance shows that the limiting case $p_d \rightarrow 1$ converges to the standard TASEP, where $p_d$ refers to the affected hopping rate due to presence of obstacle. The mean-field solutions are validated using extensive Monte Carlo simulations.
Accurate and large-scale atomic simulation is crucial for understanding growth mechanism and precise synthesis control of crystals. However, it is beyond the capacities of both ab initio and classical molecular dynamics (MD). This study takes a feasible way by developing ab initio data to accurate machine learning interatomic potentials (MLIPs) to explore the complex growth processes. It successfully simulated the growth of monolayer, bilayer MoS$_2$/WS$_2$ and MoS$_2$/WS$_2$ heterostructures under variable conditions. Importantly, a SMMS (M is Mo/W atom) is disclosed as a new 2D structure with high stability during the growth. The two-step vapor-deposition can be regarded as the alternating between MS$_2$ and SMMS structures. Besides, the intermediate SMMS structure easy forms alloys can result in production of Mo$_x$W$_{1-x}$S$_2$ alloy. Furthermore, metallic SMMS is an ideal self-intercalating electrode for transition metal dichalcogenides based FET. This research provides a mode that ab initio data driving MLIPs to efficient and accurate atomic growth simulation and design of materials.
Li-ion battery degradation processes are multi-scale, heterogeneous, dynamic and involve multiple cell components through cross talk mechanisms. Correlated operando characterization capable of measuring several key parameters are needed to accelerate understanding on these complex degradation processes. In particular, degradation mechanisms during overcharge of LiNiO2/Graphite-Silicon is well known at the material level featuring O2 gas release and concomitant surface reconstruction of LiNiO2. However, there are still debates regarding the role of high voltage O1 phase formation on gas production and no information on the effect of produced gases on the cell components (anode or sensors), or effect of overcharge on electrode level behavior. In this work, we simultaneously measured the gas produced using operando mass spectrometry while spatially resolving nanostructure and lattice changes using operando micro SAXS/WAXS mapping during the formation and over charge of a LiNiO2/Gr-Si pouch cell. This new correlated operando characterization experiment allowed to (1) confirm the absence of O1 phase even with substantial gas produced at end of charge, (2) unveil the effect of gases on reference and negative electrodes, (3) show that overcharge increases in-plane reaction heterogeneities by creating local degraded spots lagging behind the ensemble electrochemistry. These findings will be important to optimize ageing of devices based on similar chemistries, in particular Ni-rich NMC, while showing the strength of correlated characterization leading to more efficient and robust information on complex mechanisms.