Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Correlation Functions”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Short-time dynamics of 2-thiouracil in the light absorbing S{sub 2}(ππ{sup ∗}) state

Ultrahigh quantum yields of intersystem crossing to the lowest triplet state T{sub 1} are observed for 2-thiouracils (2TU), which is in contrast to the natural uracils that predominantly exhibit ultrafast internal conversion to the ground state upon excitation to the singlet excited state. The intersystem crossing mechanism of 2TU has recently been investigated using second-order perturbation methods with a high-level complete-active space self-consistent field. Three competitive nonadiabatic pathways to the lowest triplet state T{sub 1} from the initially populated singlet excited state S{sub 2} were proposed. We investigate the initial decay dynamics of 2TU from the light absorbing excited states using resonance Raman spectroscopy, time-dependent wave-packet theory in the simple model, and complete-active space self-consistent field (CASSCF) and time dependent-Becke’s three-parameter exchange and correlation functional with the Lee-Yang-Parr correlation functional (TD-B3LYP) calculations. The obtained short-time structural dynamics in easy-to-visualize internal coordinates were compared with the CASSCF(16,11) predicted key nonadiabatic decay routes. Our results indicate that the predominant decay pathway initiated at the Franck-Condon region is toward the S{sub 2}/S{sub 1} conical intersection point and S{sub 2}T{sub 3} intersystem crossing point, but not toward the S{sub 2}T{sub 2} intersystem crossing point.

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

Exploring the NΛ–NΣ coupled system with high precision correlation techniques at the LHC

The interaction of Λ and Σ hyperons (Y) with nucleons (N) is strongly influenced by the coupled-channel dynamics. Due to the small mass difference of the NΛ and NΣ systems, the sizable coupling strength of the NΣ↔NΛ processes constitutes a crucial element in the determination of the NΛ interaction. In this letter we present the most precise measurements on the interaction of pΛ pairs, from zero relative momentum up to the opening of the NΣ channel. The correlation function in the relative momentum space for pΛ ⊕$\bar{p}$$\bar{Λ}$ pairs measured in high-multiplicity triggered pp collisions at $\sqrt{s}$=13 TeV at the LHC is reported. The opening of the inelastic NΣ channels is visible in the extracted correlation function as a cusp-like structure occurring at relative momentum k*=289 MeV/c. This represents the first direct experimental observation of the NΣ↔NΛ coupled channel in the pΛ system. The correlation function is compared with recent chiral effective field theory calculations, based on different strengths of the NΣ↔NΛ transition potential. A weaker coupling, as possibly supported by the present measurement, would require a more repulsive three-body NNΛ interaction for a proper description of the Λ in-medium properties, which has implications on the nuclear equation of state and for the presence of hyperons inside neutron stars.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Cyber–Physical Correlation Effects in Defense Games for Large Discrete Infrastructures

In certain critical infrastructures, correlations between cyber and physical components can be exploited to launch strategic attacks, so that disruptions to one component may affect others and possibly the entire infrastructure. Such correlations must be explicitly taken into account in ensuring the survival of the infrastructure. For large discrete infrastructures characterized by the number of cyber and physical components, we characterize the cyber–physical interactions at two levels: (i) the cyber–physical failure correlation function specifies the conditional survival probability of the cyber sub-infrastructure given that of the physical sub-infrastructure (both specified by their marginal probabilities), and (ii) individual survival probabilities of both sub-infrastructures are characterized by first-order differential conditions expressed in terms of their multiplier functions. We formulate an abstract problem of ensuring the survival probability of a cyber–physical infrastructure with discrete components as a game between the provider and attacker, whose utility functions are composed of infrastructure survival probability terms and cost terms, both expressed in terms of the number of components attacked and reinforced. We derive Nash equilibrium conditions and sensitivity functions that highlight the dependence of infrastructure survival probability on cost terms, correlation functions, multiplier functions, and sub-infrastructure survival probabilities. We apply these analytical results to characterize the defense postures of simplified models of metro systems, cloud computing infrastructures, and smart power grids.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1975

Broadband noise generated by turbulent inflow to rotor or stator blades in an annular duct

The Green's function relating the radiated pressure field to the fluctuating forces on rotor or stator blades is developed in the presence of dissipation due to turbulent velocity fluctuations and sound speed fluctuations. The resonances in the output power spectrum which would occur at the cut-off frequencies in the absence of dissipation should be removed and smeared out by the incorporation of dissipation. Wave number dependence is developed for an effective eddy viscosity due to the aforementioned fluctuations in the background medium. The space-time correlation function for blade-normal velocity fluctuations on a single or on two different blades is developed in terms of the velocity correlation tensor for the inflow under the assumptions of isotropy and (Taylor) frozen behavior. The correlation function is then simplified under certain approximations and the behavior of the blade-force correlation function is inferred.

Lnae, F.↗
DOE OSTI2016

Relativistic effects in Lyman-α forest

We present the calculation of the Lyman-alpha (Lyman-α) transmitted flux fluctuations with full relativistic corrections to the first order. Even though several studies exist on relativistic effects in galaxy clustering, this is the first study to extend the formalism to a different tracer of underlying matter at unique redshift range (z=2−5). Furthermore, we show a comprehensive application of our calculations to the Quasar-Lyman-α cross-correlation function. Our results indicate that the signal of relativistic effects are sizeable at Baryonic Acoustic Oscillation (BAO) scale mainly due to the large differences in density bias factors of our tracers. We construct an observable, the anti-symmetric part of the cross-correlation function, that is dominated by the relativistic signal and offers a new way to measure the relativistic terms at relatively small scales. The analysis shows that relativistic effects are important when considering cross-correlations between tracers with very different biases, and should be included in the data analysis of the current and future surveys. Moreover, the idea presented in this paper is highly complementary to other techniques and observables trying to isolate the effect of the relativistic corrections and thus test the validity of the theory of gravity beyond the Newtonian regime.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Angular systematics-free cosmological analysis of galaxy clustering in configuration space

ABSTRACT Galaxy redshift surveys are subject to incompleteness and inhomogeneous sampling due to the various constraints inherent to spectroscopic observations. This can introduce systematic errors on the summary statistics of interest, which need to be mitigated in cosmological analysis to achieve high accuracy. Standard practices involve applying weighting schemes based on completeness estimates across the survey footprint, possibly supplemented with additional weighting schemes accounting for density-dependent effects. In this work, we concentrate on pure angular systematics and describe an alternative approach consisting in analysing the galaxy two-point correlation function where angular modes are nulled. By construction, this procedure removes all possible known and unknown sources of angular observational systematics, but also part of the cosmological signal. We use a modified Landy–Szalay estimator for the two-point correlation function that relies on an additional random catalogue where angular positions are randomly drawn from the galaxy catalogue, and provide an analytical model to describe this modified statistic. We test the model by performing an analysis of the full anisotropic clustering in mock catalogues of luminous red and emission-line galaxies at 0.43 < z < 1.1. We find that the model fully accounts for the modified correlation function in redshift space, without introducing new nuisance parameters. The derived cosmological parameters from the analysis of baryon acoustic oscillations and redshift-space distortions display slightly larger statistical uncertainties, mostly for the growth rate of structure parameter fσ8 that exhibits a $50{{\ \rm per\ cent}}$ statistical error increase, but free from angular systematic error.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Giant Wilson loops and AdS 2 /dCFT 1

The 1/2-BPS Wilson loop in N = 4 supersymmetric Yang-Mills theory is an important and well-studied example of conformal defect. In particular, much work has been done for the correlation functions of operator insertions on the Wilson loop in the fundamental representation. In this paper, we extend such analyses to Wilson loops in the large-rank symmetric and antisymmetric representations, which correspond to probe D3 and D5 branes with AdS 2 × S 2 and AdS 2 × S 4 worldvolume geometries, ending at the AdS 5 boundary along a one-dimensional contour. We first compute the correlation functions of protected scalar insertions from supersymmetric localization, and obtain a representation in terms of multiple integrals that are similar to the eigenvalue integrals of the random matrix, but with important differences. Using ideas from the Fermi Gas formalism and the Clustering method, we evaluate their large N limit exactly as a function of the ’t Hooft coupling. The results are given by simple integrals of polynomials that resemble the Q-functions of the Quantum Spectral Curve, with integration measures depending on the number of insertions. Next, we study the correlation functions of fluctuations on the probe D3 and D5 branes in AdS. We compute a selection of three- and four-point functions from perturbation theory on the D-branes, and show that they agree with the results of localization when restricted to supersymmetric kinematics. We also explain how the difference of the internal geometries of the D3 and D5 branes manifests itself in the localization computation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Quantum Monte Carlo formalism for dynamical pions and nucleons

In most simulations of nonrelativistic nuclear systems, the wave functions found solving the many-body Schrodinger equations describe the quantum-mechanical amplitudes of the nucleonic degrees of freedom. In those simulations the pionic contributions are encoded in nuclear potentials and electroweak currents, and they determine the low-momentum behavior. In this work we present an alternative quantum Monte Carlo formalism in which both relativistic pions and nonrelativistic nucleons are explicitly included in the quantum-mechanical states of the system. We report the renormalization of the nucleon mass as a function of the momentum cutoff, a Euclidean time density correlation function that deals with the short-time nucleon diffusion, and the pion cloud density and momentum distributions. In the two-nucleon sector we show that the interaction of two static nucleons at large distances reduces to the one-pion exchange potential, and we fit the low-energy constants of the contact interactions to reproduce the binding energy of the deuteron and two neutrons in finite volumes. We show that the method can be readily applied to light-nuclei.

Madeira, Lucas↗
DOE OSTI2021

Quantum system dynamics with a weakly nonlinear Josephson junction bath

Here, we investigate the influence of a weakly nonlinear Josephson bath consisting of a chain of Josephson junctions on the dynamics of a small quantum system (LC oscillator). Focusing on the regime where the charging energy is the largest energy scale, we perturbatively calculate the correlation function of the Josephson bath to the leading order in the Josephson energy divided by the charging energy while keeping the cosine potential exactly. When the variation of the charging energy along the chain ensures fast decay of the bath correlation function, the dynamics of the LC oscillator that is weakly and capacitively coupled to the Josephson bath can be solved through the Markovian master equation. We establish a duality relation for the Josephson bath between the regimes of large charging and Josephson energies, respectively. The results can be applied to cases where the charging energy either is nonuniformly engineered or disordered in the chain. Furthermore, we find that the Josephson bath may become non-Markovian when the temperature is increased beyond the zero-temperature limit in that the bath correlation function gets shifted by a constant and does not decay with time.

1-dimensional systems↗
DOE OSTI2021

Entropy Pair Functional Theory: Direct Entropy Evaluation Spanning Phase Transitions

We prove that, within the class of pair potential Hamiltonians, the excess entropy is a universal, temperature-independent functional of the density and pair correlation function. This result extends Henderson’s theorem, which states that the free energy is a temperature dependent functional of the density and pair correlation. The stationarity and concavity of the excess entropy functional are discussed and related to the Gibbs–Bugoliubov inequality and to the free energy. We apply the Kirkwood approximation, which is commonly used for fluids, to both fluids and solids. Approximate excess entropy functionals are developed and compared to results from thermodynamic integration. The pair functional approach gives the absolute entropy and free energy based on simulation output at a single temperature without thermodynamic integration. We argue that a functional of the type, which is strictly applicable to pair potentials, is also suitable for first principles calculation of free energies from Born–Oppenheimer molecular dynamics performed at a single temperature. In conclusion, this advancement has the potential to reduce the evaluation the free energy to a simple modification to any procedure that evaluates the energy and the pair correlation function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A redefinition of the halo boundary leads to a simple yet accurate halo model of large-scale structure

ABSTRACT We present a model for the halo–mass correlation function that explicitly incorporates halo exclusion and allows for a redefinition of the halo boundary in a flexible way. We assume that haloes trace mass in a way that can be described using a single scale-independent bias parameter. However, our model exhibits scale-dependent biasing due to the impact of halo-exclusion, the use of a ‘soft’ (i.e. not infinitely sharp) halo boundary, and differences in the one halo term contributions to ξhm and ξmm. These features naturally lead us to a redefinition of the halo boundary that lies at the ‘by eye’ transition radius from the one-halo to the two-halo term in the halo–mass correlation function. When adopting our proposed definition, our model succeeds in describing the halo–mass correlation function with $\approx 2{{\ \rm per\ cent}}$ residuals over the radial range 0.1 h−1 Mpc < r < 80 h−1 Mpc, and for halo masses in the range 1013 h−1 M⊙ < M < 1015 h−1 M⊙. Our proposed halo boundary is related to the splashback radius by a roughly constant multiplicative factor. Taking the 87 percentile as reference we find rt/Rsp ≈ 1.3. Surprisingly, our proposed definition results in halo abundances that are well described by the Press–Schechter mass function with δsc = 1.449 ± 0.004. The clustering bias parameter is offset from the standard background-split prediction by $\approx 10{{\ \rm per\ cent}}\!-\!15{{\ \rm per\ cent}}$. This level of agreement is comparable to that achieved with more standard halo definitions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Constraining the $\overline{K}$N coupled channel dynamics using femtoscopic correlations at the LHC

The interaction of K – with protons is characterised by the presence of several coupled channels, systems like $\overline{K}$ 0 n and πΣ with a similar mass and the same quantum numbers as the K – p state. The strengths of these couplings to the K – p system are of crucial importance for the understanding of the nature of the Λ(1405) resonance and of the attractive K – p strong interaction. In this article, we present measurements of the K – p correlation functions in relative momentum space obtained in pp collisions at $\sqrt{s}$=13 Te, in p–Pb collisions at $\sqrt{s_{NN}}$= 5.02 Te, and (semi)peripheral Pb–Pb collisions at $\sqrt{s_{NN}}$= 5.02 Te. The emitting source size, composed of a core radius anchored to the K + p correlation and of a resonance halo specific to each particle pair, varies between 1 and 2 fm in these collision systems. The strength and the effects of the $\overline{K}$ 0 n and πΣ inelastic channels on the measured K – p correlation function are investigated in the different colliding systems by comparing the data with state-of-the-art models of chiral potentials. A novel approach to determine the conversion weights ω, necessary to quantify the amount of produced inelastic channels in the correlation function, is presented. In this method, particle yields are estimated from thermal model predictions, and their kinematic distribution from blast-wave fits to measured data. The comparison of chiral potentials to the measured K – p interaction indicates that, while the πΣ –K – p dynamics is well reproduced by the model, the coupling to the $\overline{K}$ 0 n channel in the model is currently underestimated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

When Do Band Gap Calculations Agree with Experiments in Monolayer-Protected Cu 14 and Au 20 Atomically Precise Nanoclusters? A (TD)-DFT Comparison of HOMO–LUMO, Fundamental, Optical, and Electrochemical Energy Gaps

In view of the tremendous progress in atomically precise metal nanoclusters where electrochemical and optical energetics are routinely supported by computations to establish structure−function correlations, we explore the relationship between different protocols for measuring and computing band gaps of two distinct organic ligand-protected nanoclusters: [Cu 14 H 10 (MBN) 3 (PPh 3 ) 8 ] + and Au 20 (TBBT) 16 . Through UV/visible spectroscopy and differential pulse voltammetry, we measure optical and electrochemical band gaps in those systems. We then compare these experimentally determined gaps to HOMO−LUMO gaps, fundamental gaps, vertical excitation energies, and E o ox − E o red potentials computed using different density functional theory (DFT) or time-dependent DFT (TDDFT) methods. Specifically, in both copper and gold nanoclusters, we test the effect of truncating inert ligands from the model and compare density functionals with varying degrees of Hartree−Fock (HF) exchange from 0 to 50%, range-separated hybrids with a varying long-range tuning parameter, different correlation functionals, basis sets, and (equilibrium and nonequilibrium) continuum solvation models. Despite having different frontier orbital characters (the copper nanocluster has a metal-to-ligand charge transfer character while the gold nanocluster has metal-centered frontier orbitals), both nanoclusters display a similar sensitivity of the HOMO−LUMO gap to the HF exchange that is partially mitigated when computing the fundamental, optical, and electrochemical gaps. Other factors, such as the nature of the correlation functional, basis set, and geometry relaxation, have a considerably smaller effect on computed band gaps in these systems. Overall, this work provides guidelines for factors of varied importance for correlating computed and experimental band gap values.

Chemical calculations↗
DOE OSTI2022

Detection of spatial clustering in the 1000 richest SDSS DR8 redMaPPer clusters with nearest neighbor distributions

ABSTRACT Distances to the k-nearest-neighbor (kNN) data points from volume-filling query points are a sensitive probe of spatial clustering. Here, we present the first application of kNN summary statistics to observational clustering measurement, using the 1000 richest redMaPPer clusters (0.1 ≤ z ≤ 0.3) from the SDSS DR8 catalog. A clustering signal is defined as a difference in the cumulative distribution functions (CDFs) of kNN distances from fixed query points to the observed clusters versus a set of unclustered random points. We find that the k = 1, 2-NN CDFs of redMaPPer deviate significantly from the randoms’ across scales of 35 to 155 Mpc, which is a robust signature of clustering. In addition to kNN, we also measure the two-point correlation function for the same set of redMaPPer clusters versus random points, which shows a noisier and less significant clustering signal within the same radial scales. Quantitatively, the χ2 distribution for both the kNN-CDFs and the two-point correlation function measured on the randoms peak at χ2 ∼ 50 (null hypothesis), whereas the kNN-CDFs (χ2 ∼ 300, p = 1.54 × 10−36) pick up a much more significant clustering signal than the two-point function (χ2 ∼ 100, p = 1.16 × 10−6) when measured on redMaPPer. Finally, the measured 3NN and 4NN CDFs deviate from the predicted k = 3, 4-NN CDFs assuming an ideal Gaussian field, indicating a non-Gaussian clustering signal for redMaPPer clusters, although its origin might not be cosmological due to observational systematics. Therefore, kNN serves as a more sensitive probe of clustering complementary to the two point correlation function, providing a novel approach for constraining cosmology and galaxy–halo connection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Symmetry Breaking with the SCAN Density Functional Describes Strong Correlation in the Singlet Carbon Dimer

The SCAN (strongly constrained and appropriately normed) meta-generalized gradient approximation (meta-GGA), which satisfies all 17 exact constraints that a meta-GGA can satisfy, accurately describes equilibrium bonds that are normally correlated. With symmetry breaking, it also accurately describes some sd equilibrium bonds that are strongly correlated. While sp equilibrium bonds are nearly always normally correlated, the C 2 singlet ground state is known from correlated wave function theory to be a rare case of strong correlation in an sp equilibrium bond. Earlier work that calculated atomization energies of the molecular sequence B 2 , C 2 , O 2 , and F 2 in the local spin density approximation (LSDA), the Perdew–Burke–Ernzerhof (PBE) GGA, and the SCAN meta-GGA, without symmetry breaking in the molecule, found that only SCAN was accurate enough to reveal an anomalous under-binding for C 2 . Here, this work shows that spin symmetry breaking in singlet C 2 , which involves the appearance of net up- and down-spin densities on opposite sides (not ends) of the bond, corrects that underbinding, with a small SCAN atomization-energy error more like that of the other three molecules, suggesting that symmetry breaking with an advanced density functional might reliably describe strong correlation. This article also discusses some general aspects of symmetry breaking and the insights into strong correlation that symmetry breaking can bring. The normally correlated low-lying triplet excited state has the right vertical excitation energy in SCAN but not in LSDA or PBE, where the triplet is a false ground state. Fractional occupation numbers are found only for the symmetry-unbroken singlet and only in LSDA and PBE GGA.

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

Dark Energy Survey Year 3 results: Measurement of the baryon acoustic oscillations with three-dimensional clustering

The three-dimensional correlation function offers an effective way to summarize the correlation of the large-scale structure even for imaging galaxy surveys. We have applied the projected three-dimensional correlation function, ξ p to measure the baryonic acoustic oscillations (BAO) scale on the first-three years Dark Energy Survey data. The sample consists of about 7 million galaxies in the redshift range 0.6< z p < 1.1 over a footprint of 4108 deg 2 . Our theory modeling includes the impact of realistic true redshift distributions beyond Gaussian photo-z approximation. ξ p is obtained by projecting the three-dimensional correlation to the transverse direction. To increase the signal-to-noise of the measurements, we have considered a Gaussian stacking window function in place of the commonly used top-hat. ξ p is sensitive to D M (z eff )/r s , the ratio between the comoving angular diameter distance and the sound horizon. Using the full sample, D M (z eff )/r s is constrained to be19.00 ± 0.67 (top-hat) and 19.15 ± 0.58 (Gaussian) at z eff = 0.835. The constraint is weaker than the angular correlation w constraint 18.84 ± 0.50), and we trace this to the fact that the BAO signals are heterogeneous across redshift. While ξ p responds to the heterogeneous signals by enlarging the error bar, w can still give a tight bound on D M s in this case. When a homogeneous BAO-signal subsample in the range 0.7 < z p <1.0 (z eff = 0.845) is considered, ξp yields 19.80 ± 0.67 (top-hat) and 19.84 ± 0.53 (Gaussian). The latter is mildly stronger than the w constraint (19.86 ± 0.55). We find that the ξ p results are more sensitive to photo-z because ξ p keeps the three-dimensional clustering information causing it to be more prone to photo-z noise. The Gaussian window gives more robust results than the top-hat as the former is designed to suppress the low signal modes. ξ p and the angular statistics such as w have their own pros and cons, and they serve an important crosscheck with each other.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Block Lanczos algorithm for lattice QCD spectroscopy and matrix elements

Recent work introduced a new framework for analyzing correlation functions with improved convergence and signal-to-noise properties, as well as rigorous quantification of excited-state effects, based on the Lanczos algorithm and spurious eigenvalue filtering with the Cullum-Willoughby test. Here, we extend this framework to the analysis of correlation-function matrices built from multiple interpolating operators in lattice quantum chromodynamics (QCD) by constructing an oblique generalization of the block Lanczos algorithm, as well as a new physically motivated reformulation of the Cullum-Willoughby test that generalizes to block Lanczos straightforwardly. The resulting block Lanczos method directly extends generalized eigenvalue problem (GEVP) methods, which can be viewed as applying a single iteration of block Lanczos. Block Lanczos provides qualitative and quantitative advantages over GEVP methods analogous to the benefits of Lanczos over the standard effective mass, including faster convergence to ground- and excited-state energies, explicitly computable two-sided error bounds, straightforward extraction of matrix elements of external currents, and asymptotically constant signal-to-noise. No fits or statistical inference are required. Proof-of-principle calculations are performed for noiseless mock-data examples as well as two-by-two proton correlation-function matrices in lattice QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Adaptively Compressed Exchange Operator for Large-Scale Hybrid Density Functional Calculations with Applications to the Adsorption of Water on Silicene

Density functional theory (DFT) calculations using hybrid exchange-correlation functionals have been shown to provide an accurate description of the electronic structures of nanosystems. However, such calculations are often limited to small system sizes due to the high computational cost associated with the construction and application of the Hartree-Fock (HF) exchange operator. Here, we demonstrate that the recently developed adaptively compressed exchange (ACE) operator formulation [J. Chem. Theory Comput. 2016, 12, 2242-2249] can enable hybrid functional DFT calculations for nanosystems with thousands of atoms. The cost of constructing the ACE operator is the same as that of applying the exchange operator to the occupied orbitals once, while the cost of applying the Hamiltonian operator with a hybrid functional (after construction of the ACE operator) is only marginally higher than that associated with applying a Hamiltonian constructed from local and semilocal exchange-correlation functionals. Therefore, this new development significantly lowers the computational barrier for using hybrid functionals in large-scale DFT calculations. We demonstrate that a parallel planewave implementation of this method can be used to compute the ground-state electronic structure of a 1000-atom bulk silicon system in less than 30 wall clock minutes and that this method scales beyond 8000 computational cores for a bulk silicon system containing about 4000 atoms. The efficiency of the present methodology in treating large systems enables us to investigate adsorption properties of water molecules on Ag-supported two-dimensional silicene. Our computational results show that water monomer, dimer, and trimer configurations exhibit distinct adsorption behaviors on silicene. In particular, the presence of additional water molecules in the dimer and trimer configurations induces a transition from physisorption to chemisorption, followed by dissociation on Ag-supported silicene. This is caused by the enhanced effect of hydrogen bonds on charge transfer and proton transfer processes. Lastly, such a hydrogen bond autocatalytic effect is expected to have broad applications for silicene as an efficient surface catalyst for oxygen reduction reactions and water dissociation.

36 MATERIALS SCIENCE↗