Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COULOMB FIELD”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Mitigation of space charge effects using electron column at IOTA ring

We investigate a novel method to mitigate space charge effects of high intensity proton beams propagating in circular accelerators by means of trapping and controlling electrons generated from beam-induced residual gas ionization. This compensation method uses Coulomb repulsion force between a proton beam and electrons to mitigate self-space charge effects of the beam if it passes through a plasma column. The transverse electron-proton (e-p) instability in the plasma column is well controlled by the longitudinal magnetic field of a solenoid magnet and the bias voltages on electrodes. In this report, we will show simulation results how to control distributions of electrons and ions as well as that of the proton beam inside the column.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Probing the Boundary between Classical and Quantum Mechanics by Analyzing the Energy Dependence of Single-Electron Scattering Events at the Nanoscale

The relation between the energy-dependent particle and wave descriptions of electron–matter interactions on the nanoscale was analyzed by measuring the delocalization of an evanescent field from energy-filtered amplitude images of sample/vacuum interfaces with a special aberration-corrected electron microscope. The spatial field extension coincided with the energy-dependent self-coherence length of propagating wave packets that obeyed the time-dependent Schrödinger equation, and underwent a Goos–Hänchen shift. The findings support the view that wave packets are created by self-interferences during coherent–inelastic Coulomb interactions with a decoherence phase close to Δφ = 0.5 rad. Due to a strictly reciprocal dependence on energy, the wave packets shrink below atomic dimensions for electron energy losses beyond 1000 eV, and thus appear particle-like. Consequently, our observations inevitably include pulse-like wave propagations that stimulate structural dynamics in nanomaterials at any electron energy loss, which can be exploited to unravel time-dependent structure–function relationships on the nanoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Topological Kondo effect with spinful Majorana fermions

Motivated by the importance of studying topological superconductors beyond the mean-field approximation, we here investigate mesoscopic islands of time-reversal-invariant topological superconductors. We characterize the spectrum in the presence of strong order-parameter fluctuations in the presence of an arbitrary number of Kramers pairs of Majorana edge states and study the effect of coupling the Coulomb blockaded island to external leads. In the case of an odd fermionic parity on the island, we derive an unconventional Kondo Hamiltonian in which metallic leads couple to both topological Majorana degrees of freedom (which keep track of the parity in different leads) and the overall spin 1 2 in the island. For the simplest case of a single wire (two pairs of Majorana edge states), we demonstrate that anisotropies are irrelevant in the weak coupling renormalization group flow. This permits us to solve the Kondo problem in the vicinity of a Toulouse-type point using Abelian bosonization. We demonstrate a residual ground-state entropy of ln ( 2 ) , which is protected by spin-rotation symmetry, but reduced to ln ( 2 ) (as in the spinless topological Kondo effect) by symmetry-breaking perturbations. In the symmetric case, we further demonstrate the simultaneous presence of both Fermi-liquid and non-Fermi-liquid-like thermodynamics (depending on the observable) and derive charge and spin transport signatures of the Coulomb blockaded island. Published by the American Physical Society 2024

Materials Science↗
DOE OSTI2024

Doped moiré magnets: Renormalized flat bands and excitonic phases

Here, we explore the phase diagram of a twisted bilayer of strongly interacting electrons on a honeycomb lattice close to half filling using the slave-boson mean-field theory. Our analysis indicates that a variety of new phases can be realized as a function of chemical doping and twist angle. In particular, we find a nonmagnetic excitonic insulating phase that breaks the translational symmetry of the underlying moiré pattern. This phase results from the interplay of strong Coulomb interactions and the twist angle. In addition, we show that the features of the renormalized dispersion, such as the magic angles, depend significantly on the interactions. Our results highlight the rich physics arising in doped moiré superlattices of Mott insulators.

flat bands↗
DOE OSTI2024

A novel A-ϕ formulation for efficient electromagnetic computations

In numerical simulation of electromagnetic prob- lems, potential-based formulations have advantages over field- based formulations because of their immunity to the low- frequency breakdown catastrophe. Recently, an all-frequency stable formulation has been proposed to solve electromagnetic problems in a wide frequency band. This formulation employs magnetic vector and electric scalar potentials, as well as an auxiliary potential to enforce an inhomogeneous Coulomb gauge. However, such a three-potential formulation has encountered convergence issues at low frequencies. To avoid the convergence issue and improve computational efficiency, a novel formulation with only two potentials and an implicit Coulomb gauge is proposed in this paper. A numerical example is presented to show its solution accuracy and computational efficiency.

Mekonnen, Minyechil↗
DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Fokker-Planck simulations of fast ion ICRF and electron EC heating in a mirror plasma using CQL3D-m

The CQL3D-m continuum bounce-average Fokker-Planck code is adapted for magnetic mirror plasmas [1] and is now routinely used in no-free-parameter classical integrated modeling of mirror devices [2, 3]. In the present effort, we report on two RF methods of plasma heating in mirror machine. The fast ions (FI) are heated by Fast waves at 2nd-4th harmonic, where FIs originate from neutral beam injection at 45 degrees to the magnetic field. The scenario shows an efficient ion heating near the FI bouncing point. The electrons are heated by X-mode launched from the high magnetic field side towards the resonance. Different from the tokamak applications, CQL3D-m provides an evolving self-consistent ambipolar parallel electric field, which determines the shape of the loss cone and hence an accurate confinement time of both ions and electrons. Also, it includes a description of ion and electron sources and sinks (related to charge exchange and impact ionization) which are updated at every time step. CQL3D-m utilizes a fully nonlinear Coulomb collision operator that is important for the significantly non-Maxwellian ion distributions typically established in mirror plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nuclear-level effective theory of 𝜇→𝑒 conversion: Formalism and applications

Over the next decade new 𝜇→𝑒 conversion searches at Fermilab (Mu2e) and J-PARC (COMET, DeeMe) are expected to advance limits on charged lepton flavor violation (CLFV) by more than four orders of magnitude. Here, by considering the consequence of 𝑃 and 𝐶⁢𝑃 on elastic 𝜇→𝑒 conversion and the structure of possible charge and current densities, we show that rates are governed by six nuclear responses and a single scale, 𝑞/𝑚 𝑁 , where 𝑞 ≈ 𝑚 𝜇 is the momentum transferred from the leptons to the nucleus. To relate this result to microscopic formulations of CLFV, we construct in nonrelativistic effective theory (NRET) the CLFV nucleon-level interaction, pointing out the relevance of the dimensionless scales 𝑦=($^{𝑞⁢𝑏}_2$) 2 ⁢ > |$\overrightarrow{𝑣}_N$⁢|>|$\overrightarrow{𝑣}_μ$⁢|>|$\overrightarrow{𝑣}_T$| , where 𝑏 is the nuclear size, $\overrightarrow{𝑣}_N$ and $\overrightarrow{𝑣}_μ$ are the nucleon and muon intrinsic velocities, and $\overrightarrow{𝑣}_T$ is the target recoil velocity. We discuss previous work, noting the lack of a systematic treatment of the various small parameters. Because the parameter 𝑦 is not small, a proper calculation of 𝜇→𝑒 conversion requires a full multipole expansion of the nuclear response functions, an apparently daunting task with Coulomb-distorted electron partial waves. We demonstrate that the multipole expansion can be carried out to high precision by introducing a simplifying local momentum 𝑞 eff for the electron. Previous work has been limited to simple charge or spin interactions, thereby treating the nucleus effectively as a point particle. We show that such formulations are not compatible with the general form of the 𝜇→𝑒 conversion rate, failing to generate three of the six allowed nuclear response functions. The inclusion of the nucleon velocity $\overrightarrow{𝑣}_N$ yields an NRET with 16 operators and a rate of the general form. Consequently, in the current discovery era for CLFV, it provides the most sensible starting point for experimental analysis, defining what can and cannot be determined about CLFV from the highly exclusive process of 𝜇→𝑒 conversion. Finally, we expand the NRET operator basis to account for the effects of $\overrightarrow{𝑣}_μ$, associated with the muon's lower component, generating corrections to the CLFV coefficients of the point-nucleus response functions. Using advanced shell-model methods, we compute 𝜇→𝑒 conversion rates for a series of experimental targets, deriving bounds on the coefficients of the CLFV operators. These calculations are the first to include a general basis of CLFV operators, full evaluation of the associated nuclear response functions, and an accurate treatment of electron and muon Coulomb effects. We discuss target selection as an experimental “knob” that can be turned to probe the microscopic origins of CLFV. We describe two types of coherence that enhance certain CLFV operators and selection rules that blind elastic 𝜇→𝑒 conversion to others. We discuss the matching of the NRET onto higher level effective field theories, such as those constructed at the light quark level, noting opportunities to build on existing work in direct detection of dark matter. We discuss the relation of 𝜇→𝑒 conversion to 𝜇→𝑒+𝛾 and 𝜇→3⁢𝑒, showing how MEG II and Mu3e results will complement those of Mu2e and COMET. Finally we describe a accompanying script—in Mathematica and Python versions—that can be used to compute 𝜇→𝑒 conversion rates in various nuclear targets for the full set of NRET operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Stress Field Dynamics and Fault Slip Potential in the Paradox Basin

Abstract The Paradox Basin, straddling Utah, Colorado, Arizona, and New Mexico is characterized by an intricate amalgamation of evaporites and clastic layers and is dominated by prominent salt walls and related subsurface structures. Our research offers a new examination of the stress distribution across the basin, deriving from continuous and discrete stress measurements conducted in boreholes in the region and focal mechanism analysis, emphasizing variations over salt structures. Integrating Coulomb failure criteria with probabilistic methods, we assess potential fault movements resulting from fluid pressure alterations. Our approach provides a comprehensive understanding of the Paradox Basin's state of stress, showing a continuous change of the maximum horizontal stress orientation from N‐S at the Wasatch Fault Zone to WNW‐ESE in the northern part of the Paradox Basin and to WSW‐ENE in the southern part of the basin. Further East, into the Colorado Plateau and the Uncompahgre Uplift, the S H max orientation becomes E‐W. Decoding stress orientation dynamics has enabled critical insights into fault slip potential, especially in the basin's northern region. The salt wall faults are less likely to slip, and the Paradox Formation's evaporite and clastic rock sequence can serve as a potential low seismic risk target for carbon storage and hydrocarbon extraction.

Geochemistry & Geophysics↗
DOE OSTI2024

Dark fluxes from electromagnetic cascades

We study dark sector production in electromagnetic (EM) cascades. This problem requires accurate simulations of Standard Model (SM) and dark sector processes, both of which impact angular and energy distributions of emitted particles that ultimately determine flux predictions in a downstream detector. We describe the minimal set of QED processes which must be included to faithfully reproduce a SM cascade, and identify a universal algorithm to generate a dark sector flux given a Monte-Carlo simulation of a SM shower. We provide a new tool, "Image missing", which simulates EM cascades with associated dark vector production, and compare it against existing literature and “off the shelf” tools. The signal predictions at downstream detectors can strongly depend on the nontrivial interplay (and modelling) of SM and dark sector processes, in particular multiple Coulomb scattering and positron annihilation. We comment on potential impacts of these effects for realistic experimental setups.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Solid State Quantum Refrigeration Superconducting, Absorption and Measurement Based (Final Technical Report)

During this DOE grant, DE-SC0017890, in place for the past six years, all proposed research was carried out and published in peer-reviewed papers, as well as other projects that emerged during the research. In that effort the research team accomplished all proposed research, as well as many closely related research projects discovered and conceived of during the grant. These works included “Efficient Quantum Measurement Engines”, a work published in Physical Review Letters, giving a theory of quantum measurement-based engines, which uses quantum measurement as a resource. These engines are designed to efficiently convert energy from the stochastic quantum measurement process into useful work. Further publications include “Experimental Realization of a Quantum Dot Energy Harvester”, a joint theory and experimental work in collaboration with the group of Charles Smith in Cambridge, UK, as well as long time theoretical collaborators, Rafael Sánchez and Björn Sothmann. This work, featured as an Editor’s Suggestion in Physical Review Letters, realized an earlier theoretical proposal of ours, whereby two resonant tunneling quantum dots are connected to a central electronic cavity that is heated by a hot energy source. We also published “Superconducting Quantum Refrigerator: Breaking and Rejoining Cooper Pairs with Magnetic Field Cycles” a work done in collaboration with experimentalist Francesco Giazotto from ENS Pisa, Italy, which also resulted in a patent. This paper, published in Phys. Rev. Applied, advanced the concept of a cyclic fridge based on the normal/superconducting phase transition together with layered materials separated by tunnel junctions. We also completed the proposed research on a heat transistor, publishing “Thermal transistor and thermometer based on Coulomb-coupled conductors”, carried out as a collaboration between my group and theorists Splettstoesser (Lund U., Sweden), Sothmann (U. Duisburg-Essen, Germany), and Sánchez (U. Autónoma de Madrid, Spain). We carried out an analysis of a quantum coupled to a quantum point contact as a sensitive thermometer and heat transistor. We found the optimal statistical estimator for the temperature and compared it with experiments on the same type of devices. We also investigated autonomous quantum absorption refrigerators using quantum dots to cool by using a very hot thermal reservoir to drive heat between two other reservoirs. In the article “Quantifying the quantum heat contribution from a driven superconducting circuit”, we demonstrated that for a driven superconducting circuit, we showed heat flow provided by a hot source to the qubit can be switched on and off by varying external parameters, the frequency and the intensity of the driving. In the work “Stochastic thermodynamic cycles of a mesoscopic thermoelectric engine”, we reconsidered the autonomous thermoelectric heat engine in terms of underlying cycles. Rather than periodic behavior, the cycles were stochastic in nature. Nevertheless, by undertaking a graph theoretical analysis of the elementary transport processed, great quantitative and qualitative insight could be found. We also considered the quantum measurement process and showed that a quantum version of Maxwell’s demon could be related to the work extraction of a quantum system, closely related to arrow-of-time measures for quantum measurement, as described in our article “Thermodynamics of quantum measurement and Maxwell's demon's arrow of time”. This work was selected in Phys. Rev. A as an Editor’s Suggestion. A recent preprint titled “Cyclic Superconducting Quantum Refrigerators Using Guided Fluxon Propagation” accomplished an important piece of this grant: to propose a new kind of quantum refrigerator using the dynamics of fluxons in a type II superconductor. This invention envisioned a race-track type geometry where fluxons are confined. By applying a gradient of magnetic field together with electric current in a Corbino geometry, the circulating fluxons can actively cool a cold reservoir, realizing a new type of cyclic superconducting refrigerator. We also investigated the possibility of thermal control from different points of view. The application of quantum measurement to the system gives a new kind of control on the system of interest – we have pioneered this approach and shown that measurement can boost the thermal power of quantum engines as described in “Continuous measurement boosted adiabatic quantum thermal machines”. The ability to have heat flows on demand is an outstanding challenge, and we have provided new solutions to this problem in Thermal control across a chain of electronic nanocavities” for a chain of electron cavities using gating voltage control. The control methods using qubit/qubit coupling to create absorption fridges at their most fundamental level have also been developed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Multinode Multi-GPU Two-Electron Integrals: Code Generation Using the Regent Language

The computation of two-electron repulsion integrals (ERIs) is often the most expensive step of integral-direct self-consistent field methods. Formally it scales as O(N 4 ), where N is the number of Gaussian basis functions used to represent the molecular wave function. In practice, this scaling can be reduced to O(N 2 ) or less by neglecting small integrals with screening methods. The contributions of the ERIs to the Fock matrix are of Coulomb (J) and exchange (K) type and require separate algorithms to compute matrix elements efficiently. We previously implemented highly efficient GPU-accelerated J-matrix and K-matrix algorithms in the electronic structure code TeraChem. Although these implementations supported the use of multiple GPUs on a node, they did not support the use of multiple nodes. This presents a key bottleneck to cutting-edge ab initio simulations of large systems, e.g., excited state dynamics of photoactive proteins. We present our implementation of multinode multi-GPU J- and K-matrix algorithms in TeraChem using the Regent programming language. Regent directly supports distributed computation in a task-based model and can generate code for a variety of architectures, including NVIDIA GPUs. We demonstrate multinode scaling up to 45 GPUs (3 nodes) and benchmark against hand-coded TeraChem integral code. Finally, we also outline our metaprogrammed Regent implementation, which enables flexible code generation for integrals of different angular momenta.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations

We present the development and verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations. We apply the corrected binary collision model originally developed for charged-particles collisions to all considered scattering channels, including Coulomb collisions, elastic neutral–neutral and charged–neutral collisions, ionization, excitation, and fusion. The model's implementation is described and verified through a series of simulations, including charged- and neutral-particle thermal equilibration, slowing of electrons in warm solid-density aluminum, collisional damping of a Langmuir wave, helium gas breakdown in an applied electric field, and thermonuclear and beam–target fusion. Then, we demonstrate the model within simulations of hydrogen plasma formation in the Princeton Field-Reversed Configuration as well as of the burning of aneutronic fusion fuel p-11B. The latter includes measurement of the fusion power density in a low-density plasma and fusion production due to the stopping of a proton ignitor beam in a compressed boron target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Glassy correlated state induced by disorder in the frustrated antiferromagnet Tb 2 Zr 2 O 7

Here, we study the low-temperature thermomagnetic properties of the polycrystalline material Tb 2 Zr 2 O 7 through ac susceptibility and specific heat measurements. This zirconate displays a defect-fluorite structure in which the magnetic Tb 3+ and nonmagnetic Zr 4+ cations sit randomly on the same metal sublattice. No long-range magnetic order is found down to 100 mK, although dominant antiferromagnetic interactions are observed and the spins remain dynamic down to the lowest temperatures investigated. We observed a frequency-dependent peak around 2.5 K which is well described by many models of a canonical spin-glass transition. In-field specific heat measurements and the recovered entropy of the system R ln(4) suggest a two doublets ground state separated by 7 K. Comparisons to the pyrochlores Tb 2 Ti 2 O 7 and Tb 2 Hf 2 O 7 reinforce the collective spin-glass-type behavior and the opening for discussions of an exotic Coulomb phase in this material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Doping-driven resistive collapse of the Mott insulator in a minimal model for VO 2

Here we study the resistive collapse of the Mott insulator state in the dimer Hubbard model. This minimal model has been used to describe the physics of VO 2 , and should be relevant to other strongly correlated materials. It incorporates the physics of correlated dimers and the explicit competition between on-site Coulomb repulsion and magnetic exchange interactions. Our results unveil that between the Mott insulator at half filling and the Fermi liquid metal at high doping there is an intermediate bad metallic phase with exotic features such as a pseudogap, orbital selectivity, and a first-order metal-metal transition. The model is solved within dynamical mean field theory by means of quantum Monte Carlo, which provides the numerically exact solution of the model in the limit of large lattice dimensionality. This model can be considered as a minimal one that captures exotic phenomena associated to the physics of a doped Mott insulator, shading light on their basic physical mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nanoshell quantum dots: Quantum confinement beyond the exciton Bohr radius

Nanoshell quantum dots (QDs) represent a novel class of colloidal semiconductor nanocrystals (NCs), which supports tunable optoelectronic properties over the extended range of particle sizes. Traditionally, the ability to control the bandgap of colloidal semiconductor NCs is limited to small-size nanostructures, where photoinduced charges are confined by Coulomb interactions. A notorious drawback of such a restricted size range concerns the fact that assemblies of smaller nanoparticles tend to exhibit a greater density of interfacial and surface defects. This presents a potential problem for device applications of semiconductor NCs where the charge transport across nanoparticle films is important, as in the case of solar cells, field-effect transistors, and photoelectrochemical devices. Furthermore, the morphology of nanoshell QDs addresses this issue by enabling the quantum-confinement in the shell layer, where two-dimensional excitons can exist, regardless of the total particle size. Such a geometry exhibits one of the lowest surface-to-volume ratios among existing QD architectures and, therefore, could potentially lead to improved charge-transport and multi-exciton characteristics. The expected benefits of the nanoshell architecture were recently demonstrated by a number of reports on the CdS bulk /CdSe nanoshell model system, showing an improved photoconductivity of solids and increased lifetime of multi-exciton populations. Along these lines, this perspective will summarize the recent work on CdS bulk /CdSe nanoshell colloids and discuss the possibility of employing other nanoshell semiconductor combinations in light-harvesting and lasing applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Automatic Time Step Control to Resolve Hydromechanically Driven Fault Reactivation, Spontaneous Nucleation, and Seismic Arrest

Abstract A physical understanding of the progression from flow‐driven (quasi‐static) poromechanical deformation to dynamic fault rupture is critical to the resilient operations of several engineering systems. These processes are bridged by a progression from fault reactivation to the spontaneous nucleation of unstable sliding. Toward addressing this challenge, novel automatic time step size control methods are developed to enable accurate and efficient simulation of these dynamics and transitions from the first principles. The controllers combine local models for discretization error and Coulomb failure conditions to automatically adjust the time step size across several orders of magnitude. The methods do not require additional empirical or theoretical input and can resolve the pre‐rupture, interseismic, and seismic periods to the allowed accuracy. The computational results reveal that the proposed methods automatically capture the onset of reactivation and nucleation for homogeneous and heterogeneous fields. Hydrodynamic and structural heterogeneity lead to disparate critical nucleation sizes compared to those predicted by theory. The results highlight its potential in predicting induced seismicity in realistic subsurface engineering systems and at practical scales.

Environmental Sciences & Ecology↗
DOE OSTI2021

Direct observation of kink evolution due to Hund’s coupling on approach to metal-insulator transition in NiS 2- x Se x

Understanding characteristic energy scales is a fundamentally important issue in the study of strongly correlated systems. In multiband systems, an energy scale is affected not only by the effective Coulomb interaction but also by the Hund’s coupling. Direct observation of such energy scale has been elusive so far in spite of extensive studies. Here, we report the observation of a kink structure in the low energy dispersion of NiS 2- x Se x and its characteristic evolution with x , by using angle resolved photoemission spectroscopy. Dynamical mean field theory calculation combined with density functional theory confirms that this kink originates from Hund’s coupling. We find that the abrupt deviation from the Fermi liquid behavior in the electron self-energy results in the kink feature at low energy scale and that the kink is directly related to the coherence-incoherence crossover temperature scale. Our results mark the direct observation of the evolution of the characteristic temperature scale via kink features in the spectral function, which is the hallmark of Hund’s physics in the multiorbital system.

36 MATERIALS SCIENCE↗