CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
SEARCH · Engineering Papers
Results for “quantum mechanics and quantum computing”
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.
First-Principles-Based Study of the Decomposition of Phenol and Hydroquinone on Pt(111) Combined with Quantitative Information from XPS Spectra to Address the Impact of Coverage and Number of Hydroxyl Functional Groups
A combined first-principles-based and experimental X-ray photoelectron spectroscopy approach was used to investigate the thermal decomposition of two model biofuel compounds, phenol and hydroquinone, on Pt(111) at both low and high coverages. The DFT-based approach yields adsorption geometries and energies, activation barriers and core-level binding energy shifts for C 1s and O 1s. Increasing the coverage in the theoretical model leads to slight shifts in core-level binding energies─toward higher values for C 1s and lower values for O 1s. It also alters the energy profiles of the decomposition reaction pathway, resulting in weaker adsorption energies and changes in both reaction and activation barriers. At low temperatures, we observe a multilayer for phenol and hydroquinone upon adsorption, with desorption occurring at 200 and 270 K, respectively. Following desorption of the multilayer, decomposition proceeds via initial O–H bond scission, followed by two parallel pathways involving either C–H or C–C bond scission, whereby in the case of phenol C–H bond scission occurs first. Here, we further provide characteristic core level binding energies by theoretical calculations that are subsequently used in experimental analyses, establishing a reference database for key spectra of phenolic functionalities applicable to a range of catalytic reactions.
09 BIOMASS FUELS↗
How Solvent Mediates the Catalytic Fate of 2-Methoxyphenol in Hydrodeoxygenation Catalysis?
Kinetic and isotopic assessments combined with theoretical calculations shed light on the mechanistic differences in terms of the number of hydrogen addition steps and the electronic charges of reactive hydrogen in the catalytic sequence leading to the Ar−OCH 3 scission of 2-methoxyphenol occurring at various interfaces of vapor−Ru, cyclohexane−Ru, and water−Ru. At vapor−Ru and cyclohexane−Ru interfaces, the initial hydrogen adatom (H*) additions on the phenyl ring disrupt its aromaticity, as a step required to occur before kinetically relevant Ar−OCH 3 scission on an ensemble of Ru sites. The number of H* addition events prerequisite to the formation of the kinetically relevant Ar−OCH 3 scission transition state, however, is larger at the cyclohexane−Ru interface than those at the vapor−Ru and water−Ru interfaces. The preferential solvation of 2-methoxyphenol derived reaction precursors and H* weakens their adsorption on Ru sites, which leads to a decrease in the activation barrier for the third H* addition step and in turn a larger number of H* additions. At the water−Ru interface, water as a polar protic solvent promotes the Ar−OCH 3 scission by enabling a mechanistically distinct catalytic route an adjacent water molecule from the solvent layer assists with the intramolecular proton shuttling from the hydroxyl group (−OH) to the vicinal methoxy group (−OCH 3 ) of 2-methoxyphenol, evolving a highly charged transition state that is further stabilized by the water solvent layer, thereby leading to Ar−OCH 3 scission rate constants that are much higher than those at the vapor−Ru interface and, when comparing at the same H 2 chemical potential, higher Ar−OCH 3 scission turnovers. Furthermore, these distinct mechanistic details of the Ar−OCH 3 scission at various reaction interfaces, inferred from experiment and theory, illustrate how a solvent either solvates reaction intermediates and transition states or directly participates in the formation of kinetically relevant transition states, leading to changes in the free energy landscape, thus rewriting the catalytic fates in hydrodeoxygenation reactions.
09 BIOMASS FUELS↗
Tuning Catalytic Reactivity via Wetting Control through Oxygen Vacancies: Ru Clusters on Anatase TiO 2 and CeO 2 Supports
The shape of supported metal particles regulates their catalytic reactivity and is determined by the degree of wetting between the metal particle and the support surface. Flattened particles that wet support surfaces were reported in various catalytic systems, particularly in the subnanometer size regime. Such consequential metal–support wetting phenomena are poorly understood, and methods to study them on powder catalysts under realistic conditions are lacking. Here, we investigate the size-dependent wetting behaviors of Ru particles on two reducible-oxide supports, anatase TiO 2 (TiO 2 -A) and CeO 2 , under reducing catalytic conditions. X-ray absorption spectroscopy (XAS), low-energy ion scattering (LEIS), and density functional theory (DFT) are combined to determine the shape of Ru particles. Ru particles remain three-dimensional without wetting the TiO 2 -A support within the coverage range studied (0.06–0.98 Ru nm –2 ). In contrast, at low coverages (<0.25 Ru nm –2 ), Ru wets the CeO 2 support to form flat, disordered structures. The higher wettability of CeO 2 than TiO 2 -A is attributed to oxygen vacancies in the near-surface region. The shape difference between small Ru particles or clusters on the two supports leads to drastically contrasting catalytic reactivities in polyolefin hydrogenolysis, despite similar diameters. As a result, this work highlights the implications of metal–support wetting, or cluster shape, on catalytic behaviors of small metal clusters, while establishing the foundation for future systematic studies of such a phenomenon in realistic systems, by delivering a multitechnique methodology and revealing governing fundamental principles.
42 ENGINEERING↗
Phonon thermal Hall effect in a metallic spin ice
It has become common knowledge that phonons can generate thermal Hall effect in a wide variety of materials, although the underlying mechanism is still controversial. We study longitudinal κ xx and transverse κ xy thermal conductivity in Pr 2 Ir 2 O 7 , which is a metallic analog of spin ice. Despite the presence of mobile charge carriers, we find that both κ xx and κ xy are dominated by phonons. A T/H scaling of κ xx unambiguously reveals that longitudinal heat current is substantially impeded by resonant scattering of phonons on paramagnetic spins. Upon cooling, the resonant scattering is strongly affected by a development of spin ice correlation and κ xx deviates from the scaling in an anisotropic way with respect to field directions. Strikingly, a set of the κ xx and κ xy data clearly shows that κ xy correlates with κ xx in its response to magnetic field including a success of the T/H scaling and its failure at low temperature. This remarkable correlation provides solid evidence that an indispensable role is played by spin-phonon scattering not only for hindering the longitudinal heat conduction, but also for generating the transverse response.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Three-State Majority-vote Model on Scale-Free Networks and the Unitary Relation for Critical Exponents
Abstract We investigate the three-state majority-vote model for opinion dynamics on scale-free and regular networks. In this model, an individual selects an opinion equal to the opinion of the majority of its neighbors with probability 1 − q , and different to it with probability q . The parameter q is called the noise parameter of the model. We build a network of interactions where z neighbors are selected by each added site in the system, a preferential attachment network with degree distribution k − λ , where λ = 3 for a large number of nodes N . In this work, z is called the growth parameter. Using finite-size scaling analysis, we obtain that the critical exponents $$\beta /\bar{\nu }$$ β / ν ¯ and $$\gamma /\bar{\nu }$$ γ / ν ¯ associated with the magnetization and the susceptibility, respectively. Using Monte Carlo simulations, we calculate the critical noise parameter q c as a function of z for the scale-free networks and obtain the phase diagram of the model. We find that the critical exponents add up to unity when using a special volumetric scaling, regardless of the dimension of the network of interactions. We verify this result by obtaining the critical noise and the critical exponents for the two and three-state majority-vote model on cubic lattice networks.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
New proposed ITER divertor design using carbon insert on tungsten to mitigate ELMs and secondary radiation effects on nearby components
Abstract Building a successful device for the magnetic fusion energy production is a great challenge. ITER is an international project of the tokamak based magnetic fusion design being developed for the demonstration of the feasibility of thermonuclear technologies for future realization of successful commercial fusion energy. A key obstacle to a successful magnetic fusion energy production is however, the performance during abnormal events including plasma disruptions and edge-localized modes (ELMs). A credible reactor design must tolerate at least a few of these transient events without serious consequences such as melting of the structure. This paper investigates and compares the performance of the current ITER tokamak design during two types of transient events, i.e., ELMs occurring at normal operation and disruptions during abnormal operation. We simulated the divertor components response using our integrated 3D HEIGHTS package. The simulations include self-consistent modeling of the interaction of the released core plasma particles with the initial solid divertor material, energy deposition processes, vaporization of divertor material, secondary plasma formation and MHD evolution, incident core particles collisions and scattering from this dense secondary plasma, photon radiation of secondary plasma, and the resulting heat loads on nearby components. Our simulations showed that using a small carbon insert around the strike point can significantly reduce the overall expected damage on the tungsten dome structure, reflector plates, and prevent tungsten vaporization and its potential core plasma contamination.
47 OTHER INSTRUMENTATION↗
Upstreamness and downstreamness in input–output analysis from local and aggregate information
Abstract Ranking sectors and countries within global value chains is of paramount importance to estimate risks and forecast growth in large economies. However, this task is often non-trivial due to the lack of complete and accurate information on the flows of money and goods between sectors and countries, which are encoded in input–output (I–O) tables. In this work, we show that an accurate estimation of the role played by sectors and countries in supply chain networks can be achieved without full knowledge of the I–O tables, but only relying on local and aggregate information, e.g., the total intermediate demand per sector. Our method, based on a rank-1 approximation to the I–O table, shows consistently good performance in reconstructing rankings (i.e., upstreamness and downstreamness measures for countries and sectors) when tested on empirical data from the world input–output database. Moreover, we connect the accuracy of our approximate framework with the spectral properties of the I–O tables, which ordinarily exhibit relatively large spectral gaps. Our approach provides a fast and analytically tractable framework to rank constituents of a complex economy without the need of matrix inversions and the knowledge of finer intersectorial details.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Does boron or nitrogen substitution affect hydrogen physisorption on open carbon surfaces?
Incorporation of heteroatoms in carbon materials is commonly expected to influence their physical or chemical properties. Furthermore, contrary to previous results for methane adsorption, no technologically significant effect was identified for the hydrogen physisorption energies (measured 4.1–4.6 kJ mol –1 and calculated q st = –ΔH ads = 4.1 ± 0.7 kJ mol –1 using a comprehensive set of levels of theory) as a function of B- and N-substitution of a mid-plane C-site on open carbon surfaces.
08 HYDROGEN↗
The multivariate interaction between Au and TiO 2 colloids: the role of surface potential, concentration, and defects
The interaction between a Au particle and a defective TiO 2 surface was simulated over 6000 unique conditions, and the relative importance of the parameters was determined by logistic regression.
36 MATERIALS SCIENCE↗
Deformation of Lie–Poisson algebras and chirality
Linearization of a Hamiltonian system around an equilibrium point yields a set of Hamiltonian symmetric spectra: If λ is an eigenvalue of the linearized generator, -λ and $\barλ$ (hence, -$\barλ$) are also eigenvalues—the former implies a time-reversal symmetry, while the latter guarantees the reality of the solution. However, linearization around a singular equilibrium point (which commonly exists in noncanonical Hamiltonian systems) works out differently, resulting in breaking of the Hamiltonian symmetry of spectra; time-reversal asymmetry causes chirality. This interesting phenomenon was first found in analyzing the chiral motion of the rattleback, a boat-shaped top having misaligned axes of inertia and geometry [Z. Yoshida et al., Phys. Lett. A 381, 2772–2777 (2017)]. To elucidate how chiral spectra are generated, we study the three-dimensional Lie–Poisson systems and classify the prototypes of singularities that cause symmetry breaking. The central idea is the deformation of the underlying Lie algebra; invoking Bianchi’s list of all three-dimensional Lie algebras, we show that the so-called class-B algebras, which are produced by asymmetric deformations of the simple algebra s o ( 3 ) , yield chiral spectra when linearized around their singularities. The theory of deformation is generalized to higher dimensions, including the infinite-dimensional Poisson manifolds relevant to fluid mechanics.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Dynamics of skyrmion contraction and expansion in a magnetic film
Contraction and expansion of skyrmions in ferromagnetic films are investigated. In centrosymmetric systems, the dynamics of a collapsing skyrmion is driven by dissipation. The collapse time has a minimum on the damping constant. In systems with broken inversion symmetry, the evolution of skyrmions toward equilibrium size is driven by the Dzyaloshinskii–Moriya interaction. As a result, expressions describing the time dependence of the skyrmion size are derived and their implications for skyrmion-based information processing are discussed.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Lowest-dimensional portals to SU( N ) exotics
New matter fields charged under the strong nuclear force would have dramatic phenomenological implications. In this work, we systematically explore how these new states, which we postulate belong to some representation of SU(3) of quantum chromo-dynamics, could interact with standard model fields: We analyze all lowest-dimensional “portal” operators for any SU(3) representation, and motivated by grand unification, we extend our results to SU( N ) for N > 3. We provide a publicly available python code, “tessellation,” that automatically constructs said lowest-dimensional portal operators for any new exotic matter field charged under SU( N ) for any N .
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Resonant amplification of enzymatic chemical oscillations by oscillating flow
Using theory and simulation, we analyzed the resonant amplification of chemical oscillations that occur due to externally imposed oscillatory fluid flows. The chemical reactions are promoted by two enzyme-coated patches located sequentially on the inner surface of a pipe that transports the enclosed chemical solution. In the case of diffusion-limited systems, the period of oscillations in chemical reaction networks is determined by the rate of the chemical transport, which is diffusive in nature and, therefore, can be effectively accelerated by the imposed fluid flows. In this work, we first identify the natural frequencies of the chemical oscillations in the unperturbed reaction–diffusion system and, then, use the frequencies as a forcing input to drive the system to resonance. We demonstrate that flow-induced resonance can be used to amplify the amplitude of the chemical oscillations and to synchronize their frequency to the external forcing. In particular, we show that even 10% perturbations in the flow velocities can double the amplitude of the resulting chemical oscillations. Particularly, effective control can be achieved for the two-step chemical reactions where during the first half-period, the fluid flow accelerates the chemical flux toward the second catalytic patch, while during the second half-period, the flow amplifies the flux to the first patch. The results can provide design rules for regulating the dynamics of coupled reaction–diffusion processes and can facilitate the development of chemical reaction networks that act as chemical clocks.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Finding the order in complexity: The electronic structure of 14-1-11 zintl compounds
Yb 14 MnSb 11 and Yb 14 MgSb 11 have rapidly risen to prominence as high-performing p-type thermoelectric materials. However, the fairly complex crystal structure of A 14 MX 11 Zintl compounds renders the interpretation of the electronic band structure obscure, making it difficult to chemically guide band engineering and optimization efforts. In this work, we delineate the valence-balanced Zintl chemistry of A 14 MX 11 compounds using the molecular orbital theory. By analyzing the electronic band structures of Yb 14 MgSb 11 and Yb 14 AlSb 11 , we show that the conduction band minimum is composed of either an antibonding molecular orbital originating from the (Sb 3 ) 7– trimer or a mix of atomic orbitals of A, M, and X. The singly degenerate valence band is comprised of non-bonding Sb pz orbitals primarily from the Sb atoms in the (MSb 4 ) m– tetrahedra and of isolated Sb atoms distributed throughout the unit cell. Such a chemical understanding of the electronic structure enables strategies to engineer electronic properties (e.g., the bandgap) of A 14 MX 11 compounds.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A multi-dimensional Child–Langmuir law for any diode geometry
While prior theoretical studies of multi-dimensional space-charge limited current (SCLC) assumed emission from a small patch on infinite electrodes, none have considered emission from an entire finite electrode. In this paper, we apply variational calculus (VC) and conformal mapping, which have previously been used to derive analytic solutions for SCLC density (SCLCD) for nonplanar one-dimensional geometries, to obtain mathematical relationships for any multi-dimensional macroscopic diode with finite cathode and anode. We first derive a universal mathematical relationship between space-charge limited potential and vacuum potential for any diode and apply this technique to determine SCLCD for an eccentric spherical diode. We then apply VC and the Schwartz–Christoffel transformation to derive an exact equation for SCLCD in a general two-dimensional planar geometry with emission from a finite emitter. Particle-in-cell simulations using VSim agreed within 4%–13% for a range of ratios of emitter width to gap distance using the thinnest electrodes practical for the memory constraints of our hardware, with the difference partially attributed to the theory's assumption of infinitesimally thin electrodes. After generalizing this approach to determine SCLCD for any orthogonal diode as a function of only the vacuum capacitance and vacuum potential, we derive an analytical formulation of the three-dimensional Child–Langmuir law for finite parallel rectangular and disk geometries. These results demonstrate the utility for calculating SCLCD for any diode geometry using vacuum capacitance and vacuum potential, which are readily obtainable for many diode geometries, to guide experiment and simulation development.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Time-series forecasting using manifold learning, radial basis function interpolation, and geometric harmonics
We address a three-tier numerical framework based on nonlinear manifold learning for the forecasting of high-dimensional time series, relaxing the “curse of dimensionality” related to the training phase of surrogate/machine learning models. At the first step, we embed the high-dimensional time series into a reduced low-dimensional space using nonlinear manifold learning (local linear embedding and parsimonious diffusion maps). Then, we construct reduced-order surrogate models on the manifold (here, for our illustrations, we used multivariate autoregressive and Gaussian process regression models) to forecast the embedded dynamics. Finally, we solve the pre-image problem, thus lifting the embedded time series back to the original high-dimensional space using radial basis function interpolation and geometric harmonics. The proposed numerical data-driven scheme can also be applied as a reduced-order model procedure for the numerical solution/propagation of the (transient) dynamics of partial differential equations (PDEs). In conclusion, we assess the performance of the proposed scheme via three different families of problems: (a) the forecasting of synthetic time series generated by three simplistic linear and weakly nonlinear stochastic models resembling electroencephalography signals, (b) the prediction/propagation of the solution profiles of a linear parabolic PDE and the Brusselator model (a set of two nonlinear parabolic PDEs), and (c) the forecasting of a real-world data set containing daily time series of ten key foreign exchange rates spanning the time period 3 September 2001–29 October 2020.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Investigating boosted decision trees as a guide for inertial confinement fusion design
Inertial confined fusion experiments at the National Ignition Facility have recently entered a new regime approaching ignition. Improved modeling and exploration of the experimental parameter space were essential to deepening our understanding of the mechanisms that degrade and amplify the neutron yield. The growing prevalence of machine learning in fusion studies opens a new avenue for investigation. Here in this paper, we have applied the Gradient-Boosted Decision Tree machine-learning architecture to further explore the parameter space and find correlations with the neutron yield, a key performance indicator. We find reasonable agreement between the measured and predicted yield, with a mean absolute percentage error on a randomly assigned test set of 35.5%. This model finds the characteristics of the laser pulse to be the most influential in prediction, as well as the hohlraum laser entrance hole diameter and an enhanced capsule fabrication technique. We used the trained model to scan over the design space of experiments from three different campaigns to evaluate the potential of this technique to provide design changes that could improve the resulting neutron yield. While these data-driven model cannot predict ignition without examples of ignited shots in the training set, it can be used to indicate that an unseen shot design will at least be in the upper range of previously observed neutron yields.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗