Many-body Interaction Effects Research Articles

A successful 3D core-collapse supernova explosion model

In this paper, we present the results of our three-dimensional, multigroup, multineutrino-species radiation/hydrodynamic simulation using the state-of-the-art code fornax of the terminal dynamics of the core of a non-rotating 16 M⊙ stellar progenitor. The calculation incorporates redistribution by inelastic scattering, a correction for the effect of many-body interactions on the neutrino–nucleon scattering rates, approximate general relativity (including the effects of gravitational redshifts), velocity-dependent frequency advection, and an implementation of initial perturbations in the progenitor core. The model explodes within ∼100 ms of bounce (near when the silicon–oxygen interface is accreted through the temporarily stalled shock) and by the end of the simulation (here, ∼677 ms after bounce) is accumulating explosion energy at a rate of ∼2.5 × 1050 erg s−1. The supernova explodes with an asymmetrical multiplume structure, with one hemisphere predominating. The gravitational mass of the residual proto-neutron star at ∼677 ms is ∼1.42 M⊙. Even at the end of the simulation, explosion in most of the solid angle is accompanied by some accretion in an annular region at the wasp-like waist of the debris field. The ejecta electron fraction (Ye) is distributed between ∼0.48 and ∼0.56, with most of the ejecta mass proton-rich. This may have implications for supernova nucleosynthesis, and could have a bearing on the p- and νp-processes and on the site of the first peak of the r-process. The ejecta spatial distributions of both Ye and mass density are predominantly in wide-angle plumes and large-scale structures, but are nevertheless quite patchy.

Effect of Magnetization on the Tunneling Anomaly in Compressible Quantum Hall States.

Tunneling of electrons into a two-dimensional electron system is known to exhibit an anomaly at low bias, in which the tunneling conductance vanishes due to a many-body interaction effect. Recent experiments have measured this anomaly between two copies of the half-filled Landau level as a function of in-plane magnetic field, and they suggest that increasing spin polarization drives a deeper suppression of tunneling. Here, we present a theory of the tunneling anomaly between two copies of the partially spin-polarized Halperin-Lee-Read state, and we show that the conventional description of the tunneling anomaly, based on the Coulomb self-energy of the injected charge packet, is inconsistent with the experimental observation. We propose that the experiment is operating in a different regime, not previously considered, in which the charge-spreading action is determined by the compressibility of the composite fermions.

Many-body depletion forces of colloids in a polydisperse polymer dispersant in the long-chain limit.

We study a system of spherical non-adsorbing particles immersed in a polydisperse polymer fluid. We derive an analytic expression for the many-body depletion interactions between the colloidal particles in the limit of very long chains. We argue that this expression is essentially exact for long chains and justify this using explicit simulations. In this way we are able to elucidate the profound effect of many-body interactions on the particle thermodynamics. We show that using truncated 2-body depletion interactions leads to strong particle segregation, while the complete many-body description predicts that the total depletion force becomes weak so that the system approximates one which interacts via a so-called Kac potential. This implies that the depletion interactions can be described using mean-field theory. We show that many-body effects cause a significant contraction of the 2-phase region of the particle dispersion. We also investigate the system approaching the (tricritical) θ point, which terminates the line of first-order critical points of the polymer dispersion in a poor solvent and show that many-body effects suppress particle phase transitions.

Path integral Monte Carlo simulations of dense carbon-hydrogen plasmas.

Carbon-hydrogen plasmas and hydrocarbon materials are of broad interest to laser shock experimentalists, high energy density physicists, and astrophysicists. Accurate equations of state (EOSs) of hydrocarbons are valuable for various studies from inertial confinement fusion to planetary science. By combining path integral Monte Carlo (PIMC) results at high temperatures and density functional theory molecular dynamics results at lower temperatures, we compute the EOSs for hydrocarbons from simulations performed at 1473 separate (ρ, T)-points distributed over a range of compositions. These methods accurately treat electronic excitation effects with neither adjustable parameter nor experimental input. PIMC is also an accurate simulation method that is capable of treating many-body interaction and nuclear quantum effects at finite temperatures. These methods therefore provide a benchmark-quality EOS that surpasses that of semi-empirical and Thomas-Fermi-based methods in the warm dense matter regime. By comparing our first-principles EOS to the LEOS 5112 model for CH, we validate the specific heat assumptions in this model but suggest that the Grüneisen parameter is too large at low temperatures. Based on our first-principles EOSs, we predict the principal Hugoniot curve of polystyrene to be 2%-5% softer at maximum shock compression than that predicted by orbital-free density functional theory and SESAME 7593. By investigating the atomic structure and chemical bonding of hydrocarbons, we show a drastic decrease in the lifetime of chemical bonds in the pressure interval from 0.4 to 4 megabar. We find the assumption of linear mixing to be valid for describing the EOS and the shock Hugoniot curve of hydrocarbons in the regime of partially ionized atomic liquids. We make predictions of the shock compression of glow-discharge polymers and investigate the effects of oxygen content and C:H ratio on its Hugoniot curve. Our full suite of first-principles simulation results may be used to benchmark future theoretical investigations pertaining to hydrocarbon EOSs and should be helpful in guiding the design of future experiments on hydrocarbons in the gigabar regime.

UEfficient optimization of perturbative gadgets

Perturbative gadgets are general techniques for reducing many-body spin interactions to two-body ones using perturbation theory. This allows for potential realization of effective many-body interactions using more physically viable two-body ones. In parallel with prior work (arXiv:1311.2555 [quant-ph]), here we consider minimizing the physical resource required for implementing the gadgets initially proposed by Kempe, Kitaev and Regev (arXiv:quant-ph/0406180) and later generalized by Jordan and Farhi (arXiv:0802.1874v4). The main innovation of our result is a set of methods that efficiently compute tight upper bounds to errors in the perturbation theory. We show that in cases where the terms in the target Hamiltonian commute, the bounds produced by our algorithm are sharp for arbitrary order perturbation theory. We provide numerics which show orders of magnitudes improvement over gadget constructions based on trivial upper bounds for the error term in the perturbation series. We also discuss further improvement of our result by adopting the Schrieffer-Wolff formalism of perturbation theory and supplement our observation with numerical results.

Exchange-mediated dynamic screening in the integer quantum Hall effect regime

We study many-body interaction effects in the spatially resolved filling factor (ν) distribution for higher Landau levels (LLs) via self-consistent Hartree-Fock simulations in the integer quantum Hall (IQH) regime. Our results indicate a strong, interaction-induced tendency to avoid the simultaneous existence of partially filled spin-up and spin-down LLs. Rather, we find that such partially filled LLs consist of coexisting regions of full and empty LLs. At the boundaries between the regions of full and empty LLs, we observe edge stripes of nearly constant ν close to half-filling. This suggests that the exchange interaction induces a behavior similar to a Hund's rule for the occupation of the spin split LLs. The screening of the disorder and edge potential appears significantly reduced as compared to static Thomas-Fermi screening (Chklovskii D. B. et al., Phys. Rev. B, 46 (1992) 4026). Our results are consistent with a local, lateral ν-dependence of the exchange-enhanced spin splitting. Hence, on quantum-coherent length scales as probed here, the electron system of the IQH effect behaves similarly to a non-interacting single particle system —not because of the absence, but rather due to the dominance of many-body effects.

An investigation of the deposition of ceria on silica by quartz crystal microbalance: Observations on the effect of many body interactions

Interaction of ceria and silica is important in a wide variety of applications such as Chemical Mechanical Polish (CMP). In this work, a Quartz Crystal Microbalance (QCM) and mathematical modeling and simulations are used to investigate the dynamics of the attachment of ceria and silica under different flow rate, colloid concentration and pH conditions. The well-known DLVO forces of interaction were successfully used to distinguish between favorable and unfavorable conditions for the attachment of ceria and silica, and to develop a cleaning procedure to obtain reproducible data when using QCM for studying charged particles. The use of the optimized cleaning procedure as well as careful selection of the experimental conditions such as flow rate and concentration led to the important experimental observation that deposition of ceria particles onto silica surface is affected by many body particle–particle interactions in the bulk that can hinder the ability of the ceria particles to successfully attach to the silica surface. These experimental results could not be predicted by using a convection-diffusion-migration model which incorporates a single particle-surface DLVO interaction energy but ignores particle–particle interactions in the bulk. This is an important finding that permits the determination of the proper conditions for applications such as CMP and post CMP cleaning.

Bond order via light-induced synthetic many-body interactions of ultracold atoms in optical lattices

We show how bond order emerges due to light mediated synthetic interactions in ultracold atoms in optical lattices in an optical cavity. This is a consequence of the competition between both short- and long-range interactions designed by choosing the optical geometry. Light induces effective many-body interactions that modify the landscape of quantum phases supported by the typical Bose–Hubbard model. Using exact diagonalization of small system sizes in one-dimension, we present the many-body quantum phases the system can support via the interplay between the density and bond (or matter–wave coherence) interactions. We find numerical evidence to support that dimer phases due to bond order are analogous to valence bond states. Different possibilities of light-induced atomic interactions are considered that go beyond the typical atomic system with dipolar and other intrinsic interactions. This will broaden the Hamiltonian toolbox available for quantum simulation of condensed matter physics via atomic systems.

Quantum simulators based on the global collective light-matter interaction

We show that coupling ultracold atoms in optical lattices to quantized modes of an optical cavity leads to quantum phases of matter, which at the same time posses properties of systems with both short- and long-range interactions. This opens perspectives for novel quantum simulators of finite-range interacting systems, even though the light-induced interaction is global (i.e. infinitely long range). This is achieved by spatial structuring of the global light-matter coupling at a microscopic scale. Such simulators can directly benefit from the collective enhancement of the global light-matter interaction and constitute an alternative to standard approaches using Rydberg atoms or polar molecules. The system in the steady state of light induces effective many-body interactions that change the landscape of the phase diagram of the typical Bose-Hubbard model. Therefore, the system can support non-trivial superfluid states, bosonic dimer, trimers, etc. states and supersolid phases depending on the choice of the wavelength and pattern of the light with respect to the classical optical lattice potential. We find that by carefully choosing the system parameters one can investigate diverse strongly correlated physics with the same setup, i.e., modifying the geometry of light beams. In particular, we present the interplay between the density and bond (or matter-wave coherence) interactions. We show how to tune the effective interaction length in such a hybrid system with both short-range and global interactions.

Many-body, Pauli blocking and carrier-impurity interaction effects on the band gap of aluminum doped zinc oxide thin films: A new method to evaluate both hole and electron effective masses of degenerate semiconductors

In this study, the role of many-body interactions, Pauli blocking, and impurity-carrier interactions on the band gap variation were used to analyze the observed blue shift of the Al doped ZnO thin films. The optical constants and the band gap variations of our samples were obtained using spectroscopic ellipsometry method. Tauc-Lorentz model and effective medium approximation were used to model the dielectric function. We demonstrated for the first time that the combination of the Tauc formalism for absorption and the proposed model by Jain et al. for the band gap variation, can be used to determine both electron and hole effective masses. Reasonable results were obtained for carriers effective mass and their dependence on carrier concentration.

Beyond packing of hard spheres: The effects of core softness, non-additivity, intermediate-range repulsion, and many-body interactions on the glass-forming ability of bulk metallic glasses.

When a liquid is cooled well below its melting temperature at a rate that exceeds the critical cooling rate Rc, the crystalline state is bypassed and a metastable, amorphous glassy state forms instead. Rc (or the corresponding critical casting thickness dc) characterizes the glass-forming ability (GFA) of each material. While silica is an excellent glass-former with small Rc < 10(-2) K/s, pure metals and most alloys are typically poor glass-formers with large Rc > 10(10) K/s. Only in the past thirty years have bulk metallic glasses (BMGs) been identified with Rc approaching that for silica. Recent simulations have shown that simple, hard-sphere models are able to identify the atomic size ratio and number fraction regime where BMGs exist with critical cooling rates more than 13 orders of magnitude smaller than those for pure metals. However, there are a number of other features of interatomic potentials beyond hard-core interactions. How do these other features affect the glass-forming ability of BMGs? In this manuscript, we perform molecular dynamics simulations to determine how variations in the softness and non-additivity of the repulsive core and form of the interatomic pair potential at intermediate distances affect the GFA of binary alloys. These variations in the interatomic pair potential allow us to introduce geometric frustration and change the crystal phases that compete with glass formation. We also investigate the effect of tuning the strength of the many-body interactions from zero to the full embedded atom model on the GFA for pure metals. We then employ the full embedded atom model for binary BMGs and show that hard-core interactions play the dominant role in setting the GFA of alloys, while other features of the interatomic potential only change the GFA by one to two orders of magnitude. Despite their perturbative effect, understanding the detailed form of the intermetallic potential is important for designing BMGs with cm or greater casting thickness.

Simulations of dipolar fluids using effective many-body isotropic interactions.

The partition function of a system with pairwise-additive anisotropic dipole-dipole interactions is equal to that of a hypothetical system with many-body isotropic interactions [G. Stell, Phys. Rev. Lett. 32, 286 (1974)]. The effective many-body interactions contain n-body contributions of all orders. Each contribution is known as an expansion in terms of the particle-particle distances r, and the coefficients are temperature dependent. The leading-order two-body term is the familiar -r(-6) attraction, and the leading-order three-body term is equivalent to the Axilrod-Teller interaction. In this work, a fluid of particles with the leading-order two-body and three-body interactions is compared to an equivalent dipolar soft-sphere fluid. Molecular simulations are used to determine the conditions under which the effective many-body interactions reproduce the fluid-phase structures of the dipolar system. The effective many-body interaction works well at moderately high temperatures but fails at low temperatures where particle chaining is expected to occur. It is shown that an adjustment of the coefficients of the two-body and three-body terms leads to a good description of the structure of the dipolar fluid even in the chaining regime, due primarily to the ground-state linear configuration of the three-body Axilrod-Teller interaction. The vapor-liquid phase diagrams of systems with different Axilrod-Teller contributions are determined. As the strength of the three-body interaction is increased, the critical temperature and density both decrease and disappear completely above a threshold strength, where chaining eventually suppresses the condensation transition.

Plasmon signature in Dirac-Weyl liquids

We consider theoretically as a function of temperature the plasmon mode arising in three-dimensional Dirac liquids, i.e., systems with linear chiral relativistic single-particle dispersion, within the random phase approximation. We find that whereas no plasmon mode exists in the intrinsic (undoped) system at zero temperature, there is a well-defined finite-temperature plasmon with superlinear temperature dependence, rendering the plasmon dispersion widely tunable with temperature. The plasmon dispersion contains a logarithmic correction due to the ultraviolet-logarithmic renormalization of the electron charge, manifesting a fundamental many-body interaction effect as in quantum electrodynamics. The plasmon dispersion of the extrinsic (doped) system displays a minimum at finite temperature before it crosses over to the superlinear intrinsic behavior at higher temperature, implying that the high-temperature plasmon is a universal feature of Dirac liquids irrespective of doping. This striking characteristic temperature dependence of intrinsic Dirac plasmons along with the logarithmic renormalization is a unique manifestation of the three-dimensional relativistic Dirac nature of quasiparticle excitations and serves as an experimentally observable signature of three-dimensional Dirac materials.

Effective interactions in active Brownian suspensions.

Active colloids exhibit persistent motion, which can lead to motility-induced phase separation (MIPS). However, there currently exists no microscopic theory to account for this phenomenon. We report a first-principles theory, free of fit parameters, for active spherical colloids, which shows explicitly how an effective many-body interaction potential is generated by activity and how this can rationalize MIPS. For a passively repulsive system the theory predicts phase separation and pair correlations in quantitative agreement with simulation. For an attractive system the theory shows that phase separation becomes suppressed by moderate activity, consistent with recent experiments and simulations, and suggests a mechanism for reentrant cluster formation at high activity.

Thermodynamics of small alkali metal halide cluster ions: comparison of classical molecular simulations with experiment and quantum chemistry.

We evaluate the ability of selected classical molecular models to describe the thermodynamic and structural aspects of gas-phase hydration of alkali metal halide ions and the formation of small water clusters. To understand the effect of many-body interactions (polarization) and charge penetration effects on the accuracy of a force field, we perform Monte Carlo simulations with three rigid water models using different functional forms to account for these effects: (i) point charge nonpolarizable SPC/E, (ii) Drude point charge polarizable SWM4-DP, and (iii) Drude Gaussian charge polarizable BK3. Model predictions are compared with experimental Gibbs free energies and enthalpies of ion hydration, and with microscopic structural properties obtained from quantum DFT calculations. We find that all three models provide comparable predictions for pure water clusters and cation hydration but differ significantly in their description of anion hydration. None of the investigated classical force fields can consistently and quantitatively reproduce the experimental gas-phase hydration thermodynamics. The outcome of this study highlights the relation between the functional form that describes the effective intermolecular interactions and the accuracy of the resulting ion hydration properties.

Colloquium: Graphene spectroscopy

Spectroscopic studies of electronic phenomena in graphene are reviewed. A variety of methods and techniques are surveyed, from quasiparticle spectroscopies (tunneling, photoemission) to methods probing density and current response (infrared optics, Raman) to scanning probe nanoscopy and ultrafast pump-probe experiments. Vast complimentary information derived from these investigations is shown to highlight unusual properties of Dirac quasiparticles and many-body interaction effects in the physics of graphene.

Impurity cyclotron resonance of anomalous Dirac electrons in graphene

We have investigated a new feature of impurity cyclotron resonances common to various localized potentials of graphene. A localized potential can interact with a magnetic field in an unexpected way in graphene. It can lead to formation of anomalous boundstates that have a sharp peak with a width R in the probability density inside the potential and a broad peak of size magnetic length ℓ outside the potential. We investigate optical matrix elements of anomalous states and find that they are unusually small and depend sensitively on the magnetic field. The effect of many-body interactions on their optical conductivity is investigated using a self-consistent time-dependent Hartree–Fock approach. For a completely filled Landau level we find that an excited electron-hole pair, originating from the optical transition between two anomalous impurity states, is nearly uncorrelated with other electron-hole pairs, although it displays substantial exchange self-energy effects. This absence of correlation is a consequence of a small vertex correction in comparison to the difference between renormalized transition energies computed within the one electron-hole pair approximation. However, an excited electron-hole pair originating from the optical transition between a normal and an anomalous impurity state can be substantially correlated with other electron-hole states with a significant optical strength.

Effective field theory, three-loop perturbative expansion, and their experimental implications in graphene many-body effects

Many-body electron-electron interaction effects are theoretically considered in monolayer graphene from a continuum effective field-theoretic perspective by going beyond the standard leading-order perturbative renormalization group (RG) analysis. Given that the bare fine structure constant in graphene is of order unity, which is neither small to justify a perturbative expansion nor large enough for strong-coupling theories to be applicable, the problem is a difficult one, with some similarity to 2+1-dimensional strong-coupling quantum electrodynamics (QED). In this work, we take a systematic and comprehensive analytical approach, working primarily at the Dirac point (intrinsic graphene), by going up to three loops in the diagrammatic expansion to both ascertain the validity of perturbation theory and to estimate quantitatively higher-order renormalization effects. While no direct signatures for non-Fermi liquid behavior at the Dirac point have yet been observed experimentally, there is ample evidence for the interaction-induced renormalization of the graphene velocity as the carrier density approaches zero. We provide a critical comparison between theory and experiment, using both bare- and screened-interaction (RPA) calculations. We find that while the one-loop RG analysis gives reasonable agreement with the experimental data, especially when screening and finite-density effects are included through the RPA, the two-loop analysis reveals a strong-coupling critical point in suspended graphene, signifying either a quantum phase transition or a breakdown of the weak-coupling RG approach. We show that the latter is more likely by adapting Dyson's argument for the breakdown of perturbative QED to the case of graphene. We propose future experiments and theoretical directions to make further progress on this important and difficult problem.

Towards gauge-independent treatment of radiative capture in nuclear reactions: Applications to low-energy cluster-cluster collisions

Our departure point in describing electromagnetic (EM) interactions with nuclei (in general, bound systems of charged particles) is to use the Fock-Weyl criterion and a generalization of the Siegert theorem. It is shown how one can meet the gauge invariance principle (GIP) in all orders in the charge and construct the corresponding EM interaction operators in case of nuclear forces arbitrarily dependent on velocity. Along the guideline we have derived the conserved current density operator for a dicluster system (more precisely, the system of two finite-size clusters with many-body interaction effects included). In the context, we are addressing the current clusterization as a first step when accounting for possible cluster excitations. Being expressed through the electric and magnetic field strengths and matrix elements of the so-called generalized electric and magnetic dipole moments of the system, associated with the conserved current, the single-photon transition amplitude attains a manifestly gauge-independent (GI) form. The latter is essentially simplified at low energies.

Dirac materials

A wide range of materials, like d-wave superconductors, graphene, and topological insulators, share a fundamental similarity: their low-energy fermionic excitations behave as massless Dirac particles rather than fermions obeying the usual Schrodinger Hamiltonian. This emergent behavior of Dirac fermions in condensed matter systems defines the unifying framework for a class of materials we call "Dirac materials''. In order to establish this class of materials, we illustrate how Dirac fermions emerge in multiple entirely different condensed matter systems and we discuss how Dirac fermions have been identified experimentally using electron spectroscopy techniques (angle-resolved photoemission spectroscopy and scanning tunneling spectroscopy). As a consequence of their common low-energy excitations, this diverse set of materials shares a significant number of universal properties in the low-energy (infrared) limit. We review these common properties including nodal points in the excitation spectrum, density of states, specific heat, transport, thermodynamic properties, impurity resonances, and magnetic field responses, as well as discuss many-body interaction effects. We further review how the emergence of Dirac excitations is controlled by specific symmetries of the material, such as time-reversal, gauge, and spin-orbit symmetries, and how by breaking these symmetries a finite Dirac mass is generated. We give examples of how the interaction of Dirac fermions with their distinct real material background leads to rich novel physics with common fingerprints such as the suppression of back scattering and impurity-induced resonant states.