Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERACTION RANGE”

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 289 records · Page 16

Manipulating meso-scale solvent structure from Pd nanoparticle deposits in deep eutectic solvents

Deep Eutectic Solvents (DESs) are complex solutions that present unique challenges compared to traditional electrolytes. Unlike most electrolytes and ionic liquids, DESs have delicate hydrogen bond networks that are responsible for their highly sensitive compositional dependence on the melting point. Prior work has demonstrated unique nanoscale structure both experimentally and theoretically that brings both challenges and opportunities to their adoption in traditional electrochemical processes. In this study, we use in-situ sample-rotated ultra-small angle X-ray scattering to resolve the near-interface solvent structure after electrodepositing Pd nanoparticles onto a glassy carbon electrode in choline chloride:urea and choline chloride:ethylene glycol DESs. Our results indicate that solvent structure can be observed on the meso-scale in the choline chloride:urea and choline chloride:ethylene glycol systems. Importantly, this extended solvent structure peaks at -0.5 V (vs. Ag/AgCl) and decreases slightly with higher overpotentials during Pd electrodeposition. Experimentally, the nature of this structure is more pronounced in the ethylene glycol system, as evidenced by both the X-ray scattering as well as the electrochemical impedance spectroscopy (EIS). Here, molecular dynamics simulations and dipolar orientation analysis reveal that chloride delocalization near the Pd interface and long-range interactions between the choline and each hydrogen bond donor are very different and qualitatively consistent with the experimental data. These results show how the long-range solvent-deposit interactions can be tuned by changing the hydrogen bond donor in the DES and the applied potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Beam Dump eXperiment

Hadronic matter makes about 14% of the known universe. The remaining 86% is Dark Matter (DM). Since it does not interact with the ordinary matter via electromagnetic force, DM is not visible and, to date, it escaped detection. The search for Dark Matter (DM) is one of the hottest topic in modern physics. Despite the increasing number of astrophysical and cosmological observations proving the existence, so far no particle physics experiment has detected DM yet. Up to now, most of the experimental efforts have been focused on the so called WIMP (Weakly Interacting Massive Particle) paradigm, which predicts heavy DM (10 GeV-10 TeV mass range) interacting with Standard Model (SM) particles via the weak force mediators (W or Z bosons). More recently, due to the lack of clear evidence of WIMPs, other models of DM gained the interest of the physics community. These models consider Light DM particles (LDM), in MeV-GeV mass range. Among LDM theories, the Dark Photon theory predicts the existence of a Dark Sector interacting with SM particles via a new massive vector boson (Dark Photon, Heavy Photon or A0), mediator of a new force. This scenario, despite being theoretically well motivated, is remarkably experimentally unexplored. The Beam Dump eXperiment (BDX), is an approved experiment at Jefferson Lab (JLab), aiming to discover the DM predicted witihn the Dark Photon theory. The experiment uses the CEBAF (Continuous Electron Beam Accelerator Facility) 11 GeV electron beam, impinging on the JLab Hall-A beam-dump, to produce a beam of DM particles, detected by a ~ 1 m3 detector made of thallium doped cesium iodide (CsI(Tl)) crystals located ~ 20 m downstream. In order to achieve excellent background rejection, the experimental setup includes active vetos and passive shielding surrounding CsI(Tl) crystals. Given the weakness of the DM-SM interaction, the scattering of a DM particle in the BDX detector is a rare event. Therefore, the characterization of the expected background is a critical aspect of the experiment. Both cosmogenic and penetrating SM particles produced by the beam interaction in the dump contribute to the background of the experiment. While the cosmogenic contribution can be measured during the experiment when the beam is off, beam-related background can only be estimated via Monte Carlo (MC) simulations. Therefore, a careful assessment of possible systematics introduced by the MC is needed. The beam-correlated background characterization was made by measuring the muon flux produced by the interaction of the CEBAF beam with Hall-A dump in a dedicated experimental campaign in spring 2018 [1]. The comparison to the ?ux predicted by MC allowed us to validate the BDX simulation framework. The reach of BDX, i.e. the region that the experiment can probe in the LDM theory parameter space, depends critically on the background rejection capability and signal detection effciency. For this reason, the detector setup was fine-tuned through a dedicated study. The response to LDM and background events was evaluated for different setups and selection cuts. As a result of this procedure, the configuration resulting in the best sensitivity was selected [2]. The reach calculation was performed taking into account the different LDM production mechanisms, including the contribution of secondary particles produced in the beam-dump. In particular, the effect of the secondary positrons annihilation was found to be extremely significant. In the majority of the sensitivity studies, this contribution is neglected, but recent results [3] [4] demonstrated that this process significantly enhances the sensitivity of lepton beam-dump experiments. Currently, the BDX collaboration is focused on the deployment and operation of a small detector, called BDX-MINI, built to perform a preliminary physics measurement searching for LDM at JLab. This test will pave the way to the realization of the full BDX experiment. The measurement is currently ongoing but results are expected by the end of this year. During my PhD I was involved in all aspects of the BDX experiment: design, simulation, prototyping and data analysis. The main results of my work are reported in this thesis. This manuscript is organized as follows: the first Chapter provides an introduction to the theory of LDM, with particular attention to the Dark Photon paradigm; the second Chapter illustrates the BDX experimental setup, the LDM production and detection mechanisms and the expected backgrounds; Chapter 3 and 4 describe, respectively, the BDX-HODO measurement, with a detailed description of the simulations, and the BDX experimental setup and analysis cuts optimization. Chapter 5 reports about BDXMINI detector characterization, calibration and sensitivity estimate. Finally, Chapter 6 describes in detail the calculation of the secondary positron annihilation contribution to the sensitivity of BDX and other electron-beam thick-target experiment.

Marsicano, Luca↗

On the determination of the cosmic infrared background radiation from the high-energy spectrum of extragalactic gamma-ray sources

In a recent paper Stecker, De Jager, & Salamon have suggested using the observed approximately MeV to TeV spectra of extragalactic gamma-ray sources as probes of the local density of the cosmic infrared background radiation (CIBR) and have subsequently claimed a first possible measurement of the CIBR from the analysis of the gamma-ray spectrum of Mrk 421 (De Jager, Stecker, & Salamon). The CIBR from normal galaxies consists of two components: a stellar emission component (CIBRs), and a thermal dust emission component (CIBRd). Photons with energies in the approximately 0.1-2 TeV range interact primarily with the CIBRs, whereas interactions with CIBRd dominate the absorption of photons in the approximately 2-100 TeV energy range. SDS 92 and DSS94 considered only the interaction of the gamma-rays with the dust emission component of the CIBR. We present here an improved analysis of the absorption of extragalactic TeV gamma rays by the CIBR, taking the dual nature of its origin into account. Applying the analysis to the observed gamma-ray spectrum of Mrk 421, a BL Lac object at z = 0.031, we find agreement with DSS94 tentative evidence for absorption by the CINRs. Our analysis therefore limits the detection of the CIBR to the approximately 15-40 micron wavelength regime which, considering the uncertainties in the highest energy (greater than 4 TeV) data and ion the possibility of absorption inside the source, many turn out to be an upper limit on its energy density. At shorter wavelengths (lambda approximately = 1-15 microns), where the gamma-ray interactions are dominated by the CIBRs, our analysis definitely yields only an upper limit on the energy density of the CIBR. In contrast, DSS94 have claimed a possible first measurement of the CIBR over the entire 1-120 micron wavelength region. The upper limit on the CIBRs and tentative detection of the CIBRd are consistent with normal galaxies contributing most of the energy to the CIBR, and constrain the contribution of some exotic sources. With careful modeling of infrared foreground emissions, these constraints on the CIBR are above the values measurable by the DIRBE experiment on board the Cosmic Background Explorer (COBE) satellite.

Dwek, Eli↗

Three-dimensional checkerboard spin structure on a breathing pyrochlore lattice

The standard approach to realize a spin-liquid state is through magnetically frustrated states, relying on ingredients such as the lattice geometry, dimensionality, and magnetic interaction type of the spins. While Heisenberg spins on a pyrochlore lattice with only antiferromagnetic nearest-neighbor interactions are theoretically proven disordered, spins in real systems generally include longer-range interactions. The spatial correlations at longer distances typically stabilize a long-range order rather than enhancing a spin-liquid state. Both states can, however, be destroyed by short-range static correlations introduced by chemical disorder. Here, using disorder-free specimens with a clear long-range antiferromagnetic order, we refine the spin structure of the Heisenberg spinel ZnFe 2 O 4 through neutron magnetic diffraction. The unique wave vector (1, 0, $\frac{1}{2}$) leads to a spin structure that can be viewed as alternatively stacked ferromagnetic and antiferromagnetic tetrahedra in a three-dimensional checkerboard form. Stable coexistence of these opposing types of clusters is enabled by the bipartite breathing pyrochlore crystal structure, leading to a second-order phase transition at 10 K. The diffraction intensity of ZnFe 2 O 4 is an exact complement to the inelastic scattering intensity of several chromate spinel systems which are regarded as model classical spin liquids. Our results challenge this attribution, and suggest instead of the six-spin ring mode, spin excitations in chromate spinels are closely related to the (1, 0, $\frac{1}{2}$) type of spin order and the four-spin ferromagnetic cluster locally at one tetrahedron.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Long-Term Aging in Miscible Polymer Nanocomposites

Here we find that the initial, solvent-cast state of nanoparticles (NPs) in a polymer matrix temporally evolves during thermal annealing such that, at steady state, NPs maximize their distance from each other subject to mass balance constraints. The observed timescales for this unexpected structural reorganization, as probed by small-angle X-ray scattering, are temperature-dependent and can be prohibitively large, especially at temperatures around and below 1.2Tg. X-ray photon correlation spectroscopy measurements during reorganization reveal that the collective NP dynamics slow down with annealing at constant temperature; this is accompanied by changes in the low-frequency regime in macroscopic viscoelastic measurements in equilibrated materials. By ruling out other potential sources for these effects (i.e., electrostatic interactions, adsorbed layers), we attribute these results to a long-ranged repulsive force between the NPs caused by fluctuations in the polymer phase, i.e., the "anti-Casimir" effect proposed by Obhukhov and Semenov [Long-range interactions in polymer melts: The anti-Casimir effect. Phys Rev Lett 2005, 95 (3), 038305]. Thus, our results highlight the important role of long-term, slow NP reorganization on the structure and, subsequently, the properties of polymer nanocomposites (PNCs), even in the case of nominally miscible polymer nanoparticle hybrids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Onset of Scrambling as a Dynamical Transition in Tunable-Range Quantum Circuits

In a fast-scrambling many-body quantum system, information is spread and entanglement is built up on a time scale that grows logarithmically with the system size. This is of fundamental interest in understanding the dynamics of many-body systems, as well as in efficiently producing entangled resource states and error-correcting codes. In this work, we identify a dynamical transition marking the onset of scrambling in quantum circuits with different levels of long-range connectivity. In particular, we show that as a function of the interaction range for circuits of different structures, the tripartite mutual information exhibits a scaling collapse around a critical point between two clearly defined regimes of different dynamical behavior. We study this transition analytically in a related long-range Brownian-circuit model and show how the transition can be mapped onto the statistical mechanics of a long-range Ising model in a particular region of parameter space. This mapping predicts mean-field critical exponents v=–1/(1 + s c ), which are consistent with the critical exponents extracted from Clifford-circuit numerics. In addition to systems with conventional power-law interactions, we identify the same phenomenon in deterministic sparse circuits that can be realized in experiments with neutral-atom arrays.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interactions of particles with “continuous spin” fields

Powerful general arguments allow only a few families of long-range interactions, exemplified by gauge field theories of electromagnetism and gravity. However, all of these arguments presuppose that massless fields have zero spin scale (Casimir invariant) and hence exactly boost invariant helicity. This misses the most general behavior compatible with Lorentz symmetry. We present a Lagrangian formalism describing interactions of matter particles with bosonic “continuous spin” fields with arbitrary spin scale ρ. Remarkably, physical observables are well approximated by familiar theories at frequencies larger than ρ, with calculable deviations at low frequencies and long distances. For example, we predict specific ρ-dependent modifications to the Lorentz force law and the Larmor formula, which lay the foundation for experimental tests of the photon’s spin scale. We also reproduce known soft radiation emission amplitudes for nonzero ρ. The particles’ effective matter currents are not fully localized to their worldlines when ρ ≠ 0, which motivates investigation of manifestly local completions of our theory. Our results also motivate the development of continuous spin analogues of gravity and non-Abelian gauge theories. Given the correspondence with familiar gauge theory in the small ρ limit, we conjecture that continuous spin particles may in fact mediate known long-range forces, with testable consequences.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Three-Dimensional Tracking of Interfacial Hopping Diffusion

Theoretical predictions have suggested that molecular motion at interfaces—which influences processes including heterogeneous catalysis, (bio)chemical sensing, lubrication and adhesion, and nanomaterial self-assembly—may be dominated by hypothetical “hops” through the adjacent liquid phase, where a diffusing molecule readsorbs after a given hop according to a probabilistic “sticking coefficient.” Here, we use three-dimensional (3D) single-molecule tracking to explicitly visualize this process for human serum albumin at solid-liquid interfaces that exert varying electrostatic interactions on the biomacromolecule. Following desorption from the interface, a molecule experiences multiple unproductive surface encounters before readsorption. An average of approximately seven surface collisions is required for the repulsive surfaces, decreasing to approximately two and a half for surfaces that are more attractive. The hops themselves are also influenced by long-range interactions, with increased electrostatic repulsion causing hops of longer duration and distance. These findings explicitly demonstrate that interfacial diffusion is dominated by biased 3D Brownian motion involving bulk-surface coupling and that it can be controlled by influencing short- and long-range adsorbate-surface interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ab initio calculation of electron–phonon coupling in monoclinic β-Ga{sub 2}O{sub 3} crystal

The interaction between electrons and vibrational modes in monoclinic β-Ga{sub 2}O{sub 3} is theoretically investigated using ab-initio calculations. The large primitive cell of β-Ga{sub 2}O{sub 3} gives rise to 30 phonon modes all of which are taken into account in transport calculation. The electron-phonon interaction is calculated under density functional perturbation theory and then interpolated using Wannier–Fourier interpolation. The long-range interaction elements between electrons and polar optical phonon (POP) modes are calculated separately using the Born effective charge tensor. The direction dependence of the long-range POP coupling in a monoclinic crystal is explored and is included in the transport calculations. Scattering rate calculations are done using the Fermi golden rule followed by solving the Boltzmann transport equation using the Rode's method to estimate low field mobility. A room temperature mobility of 115 cm{sup 2}/V s is observed. Comparison with recent experimentally reported mobility is done for a wide range of temperatures (30 K–650 K). It is also found that the POP interaction dominates the electron mobility under low electric field conditions. The relative contribution of the different POP modes is analyzed and the mode 21 meV POP is found to have the highest impact on low field electron mobility at room temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Unified superfluid dark sector

We present a novel theory of a unified dark sector, where late-time cosmic acceleration emerges from the dark matter superfluid framework. The system is described by a superfluid mixture consisting of two distinguishable states with a small energy gap, such as the ground state and an excited state of dark matter. Given their contact in the superfluid, interaction between those states can happen, converting one state into the other. This long range interaction within the superfluid couples the two superfluid phonon species through a cosine potential motivated by Josephson/Rabi interactions. As a consequence of this potential, a new dynamics of late-time accelerated expansion emerges in this system, without the need of dark energy, coming from a universe containing only this two-state DM superfluid. Because the superfluid species are non-relativistic, their sound speeds remain suitably small throughout the evolution. We calculate the expansion history and growth of linear perturbations, and compare the results to ΛCDM cosmology. For the fiducial parameters studied here, the predicted expansion and growth function are close to those of ΛCDM, but the difference in the predicted growth rate is significant at late times. The present theory nicely complements the recent proposal of dark matter superfluidity to explain the empirical success of MOdified Newtonian Dynamics (MOND) on galactic scales, thus offering a unified framework for dark matter, dark energy, and MOND phenomenology.

79 ASTRONOMY AND ASTROPHYSICS↗

Analyzing the Fierz rearrangement freedom for local chiral two-nucleon potentials

Chiral effective field theory is a framework to derive systematic nuclear interactions. It is based on the symmetries of quantum chromodynamics and includes long-range pion physics explicitly, while shorter-range physics is expanded in a general operator basis. The number of low-energy couplings at a particular order in the expansion can be reduced by exploiting the fact that nucleons are fermions and therefore obey the Pauli exclusion principle. The antisymmetry permits the selection of a subset of the allowed contact operators at a given order. When local regulators are used for these short-range interactions, however, this “Fierz rearrangement freedom” is violated. In this paper, we investigate the impact of this violation at leading order (LO) in the chiral expansion. We construct LO and next-to-leading order (NLO) potentials for all possible LO-operator pairs and study their reproduction of phase shifts, the 4 He ground-state energy, and the neutron-matter energy at different densities. We demonstrate that the Fierz rearrangement freedom is partially restored at NLO where subleading contact interactions enter. Here we also discuss implications for local chiral three-nucleon interactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Radio frequency wave interactions with a plasma sheath: The role of wave and plasma sheath impedances

Radio frequency (RF) sheaths form near surfaces where plasma and strong RF fields coexist. The effect of these RF sheaths on wave propagation near the boundary can be characterized by an effective sheath impedance that includes both resistive and capacitive contributions describing RF sheath rectification and RF power absorption in the sheath [J. R. Myra and D. A. D'Ippolito, Phys. Plasmas 22, 062507 (2015)]. We define a dimensionless parameter, the ratio of incoming wave impedance to the sheath impedance, which determines the characteristics of the interaction, ranging from quasi-conducting to quasi-insulating, or in the case of matched impedances, to either perfect absorption or a sheath-plasma resonance. A semi-analytical analysis is carried out for electrostatic slow waves in the ion cyclotron range of frequencies. For the propagating slow wave case, where the incident wave is partially reflected, the fraction of power dissipated in the sheath is calculated. For the evanescent slow wave case, which admits a sheath-plasma resonance, an amplification factor is calculated. Using the impedance ratio approach, RF sheath interactions are characterized for a range of RF wave and plasma parameters including plasma density, magnetic field angle with respect to the surface, wave frequency, and wave-vector components tangent to the surface. For a particularly interesting example case, results are compared with the rfSOL code [H. Kohno and J. R. Myra, Comput. Phys. Commun. 220, 129 (2017)]. Finally, electromagnetic effects, absent from the semi-analytical analysis, are assessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Proton radiation effects in indium oxide using cascade molecular dynamics simulations

Metal oxide (MO) semiconductors, characterized by their wide band gaps and notable charge transport properties, are promising candidates for electronic applications in extreme environments, including near-Earth space. However, atomistic simulations of radiation–matter interactions in MOs remain challenging due to the limitations of existing interatomic potentials, which often fail to capture both the short-range repulsive forces essential for radiation damage modeling and the long-range electrostatic effects governing defect evolution. In this work, we develop a customized interatomic potential tailored for radiation damage simulations in indium oxide (In 2 O 3 ) as a model system, a representative MO material. Our potential integrates the Ziegler-Biersack-Littmark potential to accurately describe short-range interactions with Buckingham and Coulombic potentials to account for long-range forces. We perform molecular dynamics simulations of low-energy proton irradiation using this custom potential. We employ the primary knock-on atom (PKA) cascade method to study atomic displacements and primary defect formation. Simulations were conducted for 1 keV proton irradiation in four randomly chosen directions, and PKA-driven defect analyses at 5, 10, and 15 keV to examine the effects of direction and energy level on damage generation. Our results provide insight into the impact of irradiation direction and energy level on the cascade evolution and defect formation mechanisms.

Atomistic simulations↗

Hall conductance and the statistics of flux insertions in gapped interacting lattice systems

We study charge transport for zero-temperature infinite-volume gapped lattice systems in two dimensions with short-range interactions. We show that the Hall conductance is locally computable and is the same for all systems that are in the same gapped phase. Here, we provide a rigorous version of Laughlin’s flux-insertion argument, which shows that for short-range entangled systems, the Hall conductance is an integer multiple of e 2 / h . We show that the Hall conductance determines the statistics of flux insertions. For bosonic short-range entangled systems, this implies that the Hall conductance is an even multiple of e 2 / h . Finally, we adapt a proof of quantization of the Thouless charge pump to the case of infinite-volume gapped lattice systems in one dimension.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Melting phase diagram of bubble phases in high Landau levels

A low-disorder, two-dimensional electron system (2DES) subjected to a large perpendicular magnetic field and cooled to very low temperatures provides a rich platform for studies of many-body quantum phases. The magnetic field quenches the electrons' kinetic energy and quantizes the energy into a set of Landau levels, allowing the Coulomb interaction to dominate. In excited Landau levels, the fine interplay between short- and long-range interactions stabilizes bubble phases, Wigner crystals with more than one electron per unit cell. In this work, we present the screening properties of bubble phases, probed via a simple capacitance technique where the 2DES is placed between a top and a bottom gate and the electric field penetrating through the 2DES is measured. The bubbles formed at very low temperatures screen the electric field poorly as they are pinned by the residual disorder potential, allowing a large electric field to reach the top gate. As the temperature is increased, the penetrating electric field decreases and, surprisingly, exhibits a pronounced minimum at a temperature that appears to coincide with the melting temperature of the bubble phase. We deduce a quantitative phase diagram, as a function of Landau level filling factor (ν) and temperature, for the transition from the bubble to liquid phases for 4 ≤ ν ≤ 5.

2-dimensional systems↗

Demonstration of long-range correlations via susceptibility measurements in a one-dimensional superconducting Josephson spin chain

Abstract Spin chains have long been considered an effective medium for long-range interactions, entanglement generation, and quantum state transfer. In this work, we explore the properties of a spin chain implemented with superconducting flux circuits, designed to act as a connectivity medium between two superconducting qubits. The susceptibility of the chain is probed and shown to support long-range, cross-chain correlations. In addition, interactions between the two end qubits, mediated by the coupler chain, are demonstrated. This work has direct applicability in near term quantum annealing processors as a means of generating long-range, coherent coupling between qubits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Grassmann phase space methods for fermions. II. Field theory

In both quantum optics and cold atom physics, the behaviour of bosonic photons and atoms is often treated using phase space methods, where mode annihilation and creation operators are represented by c-number phase space variables, with the density operator equivalent to a distribution function of these variables. The anti-commutation rules for fermion annihilation, creation operators suggests the possibility of using anti-commuting Grassmann variables to represent these operators. However, in spite of the seminal work by Cahill and Glauber and a few applications, the use of Grassmann phase space methods in quantum-atom optics to treat fermionic systems is rather rare, though fermion coherent states using Grassmann variables are widely used in particle physics. This paper presents a phase space theory for fermion systems based on distribution functionals, which replace the density operator and involve Grassmann fields representing anti-commuting fermion field annihilation, creation operators. It is an extension of a previous phase space theory paper for fermions (Paper I) based on separate modes, in which the density operator is replaced by a distribution function depending on Grassmann phase space variables which represent the mode annihilation and creation operators. This further development of the theory is important for the situation when large numbers of fermions are involved, resulting in too many modes to treat separately. Here Grassmann fields, distribution functionals, functional Fokker–Planck equations and Ito stochastic field equations are involved. Typical applications to a trapped Fermi gas of interacting spin 1/2 fermionic atoms and to multi-component Fermi gases with non-zero range interactions are presented, showing that the Ito stochastic field equations are local in these cases. For the spin 1/2 case we also show how simple solutions can be obtained both for the untrapped case and for an optical lattice trapping potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Describing the Diverse Geometries of Gold from Nanoclusters to Bulk—A First-Principles-Based Hybrid Bond-Order Potential

Molecular dynamics simulations using empirical force fields (EFFs) are crucial for gaining fundamental insights into atomic structure and long time scale dynamics of Au nanoclusters with far-reaching applications in energy and devices. This approach is thwarted by the failure of currently available EFFs in describing the size-dependent dimensionality and diverse geometries exhibited by Au clusters (e.g., planar structures, hollow cages, tubes, pyramids, space-filled structures). Here, we mitigate this issue by introducing a new hybrid bond-order potential (HyBOP), which accounts for (a) short-range interactions via Tersoff-type BOP terms that accurately treat bond directionality and (b) long-range dispersion effects by a scaled Lennard–Jones term whose contribution depends on the local atomic density. We optimized the independent parameters for our HyBOP using a global optimization scheme driven by genetic algorithms. Moreover, to ensure good transferability of these parameters across different length scales, we used an extensive training data set that encompasses structural and energetic properties of 1000 13-atom Au clusters, surface energies, as well as bulk polymorphs, obtained from density functional theory (DFT) calculations. Our newly developed HyBOP has been found to accurately describe (a) global minimum energy configurations at different cluster sizes as well as order of stability of various cluster configurations at any size, (b) critical size of transition from planar to globular clusters, (c) evolution of structural motifs with cluster size, and (c) thermodynamics, structure, elastic properties, and energetic ordering of bulk condensed phases as well as surfaces, in excellent agreement with DFT calculations and spectroscopic experiments. This makes our newly developed HyBOP a valuable, computationally robust but inexpensive tool for investigating a wide range of materials phenomena occurring in Au at the atomistic level.

36 MATERIALS SCIENCE↗