Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quantum chemistry methods”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

CO adsorption on MnO(100): Experimental benchmarks compared to DFT

In this study, CO adsorption on the MnO(100) surface was studied using temperature programmed desorption (TPD) and density functional theory (DFT). TPD results show that CO is weakly-bound on MnO(100), with an experimental adsorption energy of -35.6 ± 2.1 kJ/mol at terrace sites in the low coverage limit. PBE simulations suggest that CO adsorption causes an implausible (2 × 2) surface reconstruction. PBE+U simulations show no signs of surface reconstruction, and provide an accurate estimate of the adsorption energy (-36.4 kJ/mol) when combined with the DFT-D3 method with Becke-Jonson damping to correct for van der Waals interactions. This simulation also shows that CO adsorbs C-down onto the Mn 2+ terrace site in a tilted geometry, which is also observed experimentally and computationally on the similarly-structured NiO(100) transition metal oxide surface. TPD results for large doses show a plateauing of the coverage at about 0.4 monolayers of CO at 85 K, with a defect coverage equivalent to 0.08 monolayers. Adsorption associated with defect sites is indicated by a high-temperature desorption tail which is not satisfactorily explained by DFT simulations of simple step or oxygen vacancy defects.

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

AFM-IR and s-SNOM-IR measurements of chemically addressable monolayers on Au nanoparticles

The performance of catalysts depends on their nanoscale properties, and local variations in structure and composition can have a dramatic impact on the catalytic reactivity. Therefore, probing the localized reactivity of catalytic surfaces using high spatial resolution vibrational spectroscopy, such as infrared (IR) nanospectroscopy and tip-enhanced Raman spectroscopy, is essential for mapping their reactivity pattern. Two fundamentally different scanning probe IR nanospectroscopy techniques, namely, scattering-type scanning near-field optical microscopy (sSNOM) and atomic force microscopy-infrared spectroscopy (AFM-IR), provide the capabilities for mapping the reactivity pattern of catalytic surfaces with a spatial resolution of ~20 nm. Herein, we compare these two techniques with regard to their applicability for probing the vibrational signature of reactive molecules on catalytic nanoparticles. For this purpose, we use chemically addressable self-assembled molecules on Au nanoparticles as model systems. We identified significant spectral differences depending on the measurement technique, which originate from the fundamentally different working principles of the applied methods. While AFM-IR spectra provided information from all the molecules that were positioned underneath the tip, the s-SNOM spectra were more orientation-sensitive. Due to its field-enhancement factor, the s-SNOM spectra showed higher vibrational signals for dipoles that were perpendicularly oriented to the surface. The s-SNOM sensitivity to the molecular orientation influenced the amplitude, position, and signal-to-noise ratio of the collected spectra. Ensemble-based IR measurements verified that differences in the localized IR spectra stem from the enhanced sensitivity of s-SNOM measurements to the adsorption geometry of the probed molecules.

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

Identification of Active Sites of Pure and Nitrogen-Doped Carbon Materials for Oxygen Reduction Reaction Using Constant-Potential Calculations

Nitrogen-doped carbon materials are promising metal-free catalysts for the electrochemical oxygen reduction reaction (ORR). A better theoretical understanding on the nature of the active site(s) would help further optimization of their ORR activity. Although quantum mechanical calculations have been widely employed to elucidate the active sites over various catalysts, these calculations are typically done assuming constantcharge conditions rather than the experimentally relevant constant-potential conditions. In this study, we employ the double-reference method to simulate the energetics of the ORR over pure and N-doped carbon materials under constant-potential conditions. We demonstrate that constant-potential calculations enable more accurate theoretical predictions, comparing well with existing experiments. Our key findings are (1) the zigzag edge of pure graphite is highly active for ORR, (2) the pyridinic N-doped armchair edge is highly active for ORR in alkaline media but not in acid, and (3) graphitic N can donate electrons to pyridinic N to enhance the ORR activity. Furthermore, these fundamental insights provide guidelines for the design of better carbonbased ORR catalysts.

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

Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory

Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe 2+ , Fe 3+ , and Co 3+ , including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. Furthermore, we find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method—in terms of computational cost and in terms of accuracy on monatomic systems—and therefore we chose this method for molecular benchmarks. Then Kohn–Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF 2 , FeCl 2 , and CoF 2 . We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.

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

Active learning for robust, high-complexity reactive atomistic simulations

Machine learned reactive force fields based on polynomial expansions have been shown to be highly effective for describing simulations involving reactive materials. Nevertheless, the highly flexible nature of these models can give rise to a large number of candidate parameters for complicated systems. In these cases, reliable parameterization requires a well-formed training set, which can be difficult to achieve through standard iterative fitting methods. In this paper, we present an active learning approach based on cluster analysis and inspired by Shannon information theory to enable semi-automated generation of informative training sets and robust machine learned force fields. The use of this tool is demonstrated for development of a model based on linear combinations of Chebyshev polynomials explicitly describing up to four-body interactions, for a chemically and structurally diverse system of C/O under extreme conditions. We show that this flexible training database management approach enables development of models exhibiting excellent agreement with Kohn–Sham density functional theory in terms of structure, dynamics, and speciation.

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

High-level coupled-cluster energetics by Monte Carlo sampling and moment expansions: Further details and comparisons

We recently proposed a novel approach to converging electronic energies equivalent to high-level coupled-cluster (CC) computations by combining the deterministic CC(P;Q) formalism with the stochastic configuration interaction (CI) and CC Quantum Monte Carlo (QMC) propagations. This article extends our initial study [J. E. Deustua, J. Shen, and P. Piecuch, Phys. Rev. Lett. 119, 223003 (2017)], which focused on recovering the energies obtained with the CC method with singles, doubles, and triples (CCSDT) using the information extracted from full CI QMC and CCSDT-MC, to the CIQMC approaches truncated at triples and quadruples. It also reports our first semi-stochastic CC(P;Q) calculations aimed at converging the energies that correspond to the CC method with singles, doubles, triples, and quadruples (CCSDTQ). The ability of the semi-stochastic CC(P;Q) formalism to recover the CCSDT and CCSDTQ energies, even when electronic quasi-degeneracies and triply and quadruply excited clusters become substantial, is illustrated by a few numerical examples, including the F–F bond breaking in F 2 , the automerization of cyclobutadiene, and the double dissociation of the water molecule.

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

Analysis of the Amorphous and Interphase Influence of Comononomer Loading on Polymer Properties toward Forwarding Bioadvantaged Copolyamides

In this paper we present an approach for selectively modifying the properties of semicrystalline polymers by introducing “bioadvantaged” counits. With this approach, the unique functionality of biomass can be leveraged to tailor the properties of the amorphous phase of semicrystalline polymers with minimal impact on crystallinity and thermomechanical properties. As a model case, PA 6,6 copolyamides were produced using the bioadvantaged monomer trans-3-hexenedioic acid (t3HDA). The analogous structure of t3HDA to adipic acid, a PA 6,6 monomer, allows for seamless integration. Screening over the entire composition range identified the t3HDA loading (20 mol%) beyond which properties deviate appreciably from Nylon 6,6. Once identified, copolyamides of suitable compositions were upgraded to commercial quality and fully characterized to assess the influence of counit loading and polymer structure on thermal and mechanical properties. Samples were characterized using gel permeation chromatography (GPC), proton nuclear magnetic resonance spectroscopy (1H NMR), heteronuclear single quantum coherence spectroscopy (HSQC), wide-angle X-ray scattering (WAXS), differential scanning calorimetry (DSC), thermogravimetric analysis (TGA), dynamic mechanical analysis (DMA), tensile testing, flexural testing, and water absorption testing. t3HDA units were shown to hydrate during the harsh polycondensation to 3-hydroxyhexanedioic acid (3HHDA) and fully incorporate into the polymer backbone. Loading levels up to 20% were shown to have comparable thermal and mechanical properties in the dry state, yet moisture absorption—a known method for improving the toughness, yield strain, and elongation of polyamides—was enhanced by over 100% at 20% loading. This case study on bioadvantaged copolymers elucidates the governing structure-function principles that can be leveraged to forward value-added renewable polymers.

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

Scaling Laws of Exciton Recombination Kinetics in Low Dimensional Halide Perovskite Nanostructures

Carrier recombination is a crucial process governing the optical properties of a semiconductor. Although various theoretical approaches have been utilized to describe carrier behaviors, a quantitative understanding of the impact of defects and interfaces in low dimensional semiconductor systems is still elusive. Here, we develop a model system consisting of chemically tunable, highly luminescent halide perovskite nanocrystals to illustrate the role of carrier diffusion and material dimensionality on the carrier recombination kinetics and luminescence efficiency. Our advanced synthetic methods provide a well-controlled colloidal system consisting of nanocrystals with different aspect ratios, halide compositions, and surface conditions. Using this system, we reveal the scaling laws of photoluminescence quantum yield and radiative lifetime with respect to the aspect ratio of nanocrystals. The scaling laws derived herein are not only a phenomenological observation but proved a powerful tool disentangling the carrier dynamics of microscopic systems in a quantitative and interpretable manner. The investigation of our model system and theoretical formulation bring to light the dimensionality, as a hidden constraint on carrier dynamics, and identify the diffusion length as an important parameter that distinguishes nanoscale and macroscale carrier behaviors. Furthermore, the conceptual distinction in carrier dynamics in different dimensionality regimes informs new design rules for optical devices where complex microstructures are involved.

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

Quantum stabilization and flat hydrogen-based bands of nitrogen-doped lutetium hydride

We explore the electronic and structural properties of Lu-H-N structures with specific N,H ordering due to their remarkably narrow hydrogen-based bands at the Fermi level. Although exhibits a classical instability persisting up to 17 GPa, it is anharmonically stable near ambient pressure when accounting for quantum nuclear effects. The presence of flat bands near is understood to arise from destructive quantum interference between and surrounding orbitals, with certain types of defects leaving the flat bands unaffected. The results suggest that if this system is superconducting, there is an optimal pressure near ambient where the superconducting is maximized by anharmonically stabilized low-frequency and nonadiabatically coupled high-frequency hydrogen modes. Despite the metastability of this structure, its dynamical stability when modeled beyond a classical harmonic approach suggests the narrow bands near do not cause structural distortions in certain Lu-H-N stoichiometries.

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

Competition between Hückel’s Rule and Jahn–Teller Distortion in Small Carbon Rings: A Quantum Monte Carlo Study

Structural properties and energetics of carbon rings are studied with the diffusion Monte Carlo (DMC) method. Our DMC-based geometry optimization reveals that both polyynic C 4n and cumulenic C 4n+2 rings exhibit bond length alternations for n >= 3, which is understood to be due to Jahn-Teller distortions. The bond length alternation even in a cumulenic (4n + 2) carbon ring was experimentally observed in a recently synthesized C 18 molecule. From a comparison of the DMC cohesive energies of C 4n with those of C 4n + 2 , we present a comprehensive picture of the competition between Huckel's rule and Jahn-Teller distortion in small carbon rings; the former is more dominant than the latter for n < 5 where C 4n + 2 rings are more stable than C 4n , while C 4n rings are as stable as C 4n + 2 for n < 5 where dimerization effects due to Jahn-Teller distortion are more important.

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

Optical properties of a diamond NV color center from capped embedded multiconfigurational correlated wavefunction theory

Diamond defects are among the most promising qubits. Modeling their properties through accurate quantum mechanical simulations can further their development into robust units of information. We use the recently developed capped density functional embedding theory (capped-DFET) with the multiconfigurational n-electron valence second-order perturbation theory to characterize the electronic excitation energies for different spin manifolds of the well-characterized negatively charged substitutional N defect adjacent to a vacancy (V C ) in diamond (N C V C − ). We successfully reproduce vertical excitation energies for both triplet and singlet states of N C V C − with errors < 0.1 eV. Unlike other embedding methods, capped-DFET exhibits robust predictions that are approximately independent of the embedded cluster size: it only requires a cluster to contain the defect atoms and their nearest neighbors (as small as a 40-atom capped cluster). Furthermore, our method is free from slowly converging Coulomb interactions between charged defects, and thus also only weakly dependent on supercell size.

Chemistry↗
DOE OSTI2023

Accelerating Quantum Materials Development with Advances in Transmission Electron Microscopy

Quantum materials are driving a technology revolution in sensing, communication, and computing, while simultaneously testing many core theories of the past century. Materials such as topological insulators, complex oxides, superconductors, quantum dots, color center-hosting semiconductors, and other types of strongly correlated materials can exhibit exotic properties such as edge conductivity, multiferroicity, magnetoresistance, superconductivity, single photon emission, and optical-spin locking. These emergent properties arise and depend strongly on the material’s detailed atomic-scale structure, including atomic defects, dopants, and lattice stacking. In this review, we describe how progress in the field of electron microscopy (EM), including in situ and in operando EM, can accelerate advances in quantum materials and quantum excitations. We begin by describing fundamental EM principles and operation modes. We then discuss various EM methods such as (i) EM spectroscopies, including electron energy loss spectroscopy (EELS), cathodoluminescence (CL), and electron energy gain spectroscopy (EEGS); (ii) four-dimensional scanning transmission electron microscopy (4D-STEM); (iii) dynamic and ultrafast EM (UEM); (iv) complementary ultrafast spectroscopies (UED, XFEL); and (v) atomic electron tomography (AET). We describe how these methods could inform structure–function relations in quantum materials down to the picometer scale and femtosecond time resolution, and how they enable precision positioning of atomic defects and high-resolution manipulation of quantum materials. For each method, we also describe existing limitations to solve open quantum mechanical questions, and how they might be addressed to accelerate progress. Among numerous notable results, our review highlights how EM is enabling identification of the 3D structure of quantum defects; measuring reversible and metastable dynamics of quantum excitations; mapping exciton states and single photon emission; measuring nanoscale thermal transport and coupled excitation dynamics; and measuring the internal electric field and charge density distribution of quantum heterointerfaces- all at the quantum materials’ intrinsic atomic and near atomic-length scale. Finally, we conclude by describing open challenges for the future, including achieving stable sample holders for ultralow temperature (below 10K) atomic-scale spatial resolution, stable spectrometers that enable meV energy resolution, and high-resolution, dynamic mapping of magnetic and spin fields. With atomic manipulation and ultrafast characterization enabled by EM, quantum materials will be poised to integrate into many of the sustainable and energy-efficient technologies needed for the 21st century.

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

Revisiting Competing Paths in Electrochemical CO 2 Reduction on Copper via Embedded Correlated Wavefunction Theory

We re-evaluate two key steps in the mechanism of CO 2 reduction on copper at a higher level of theory capable of correcting inherent errors in density functional theory (DFT) approximations, namely, embedded correlated wavefunction (ECW) theory. Here, we consider the CO reduction step on Cu(111), which is critical to understanding reaction selectivity. Here, we optimize embedding potentials at the periodic plane-wave DFT level using density functional embedding theory (DFET). All possible adsorption sites (adsites) for each adsorbate then are screened with ECW theory at the catalytically active site to refine the local electronic structure. Unsurprisingly, DFT and ECW theory predict different adsite preferences, largely because of DFT’s inability to properly situate the CO 2π* level. Differing preferred adsites suggest that different reaction pathways could emerge from DFT versus ECW theory. Starting from these preferred ECW theory adsites, we then obtain reaction pathways at the plane-wave DFT level using the climbing-image nudged elastic band method to determine minimum energy paths. Thereafter, we perform ECW calculations at the catalytically active site to correct the energetics at each interpolated structure (image) along the reaction pathways. Via this approach, we confirm that the first step in CO reduction via hydrogen transfer on Cu(111) is to form hydroxymethylidyne (*COH) instead of formyl (*CHO). Although the prediction to preferentially form *COH is consistent with that of DFT, the two theories predict quite different structural and mechanistic behaviors, suggesting that verification is needed for other parts of the mechanism of CO 2 reduction, which is the subject of ongoing work.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Transition-Metal-Related Quantum Emitters in Wurtzite AlN and GaN

Transition-metal centers exhibit a paramagnetic ground state in wide-bandgap semiconductors and are promising for nanophotonics and quantum information processing. Specifically, there is a growing interest in discovering prominent paramagnetic spin defects that can be manipulated using optical methods. Here, we investigate the electronic structure and magneto-optical properties of Cr and Mn substitutional centers in wurtzite AlN and GaN. We use state-of-the-art hybrid density functional theory calculations to determine level structure, stability, optical signatures, and magnetic properties of these centers. The excitation energies are calculated using the constrained occupation approach and rigorously verified with the complete active space configuration interaction approach. Our simulations of the photoluminescence spectra indicate that $Cr$$^{1+}_{Al}$ in AlN and $Cr$$^{1+}_{Ga}$ in GaN are responsible for the observed narrow quantum emission near 1.2 eV. We compute the zero-field splitting (ZFS) parameters and outline an optical spin polarization protocol for $Cr$$^{1+}_{Al}$ and $Cr$$^{1+}_{Ga}$. Our results demonstrate that these centers are promising candidates for spin qubits.

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

Controlling Molecular Structure and Spin with Multiconfigurational Quantum Chemistry (Final Technical Report)

For many first-row transition metal complexes, structure-property relationships can be obtained from high-level molecular geometry optimizations and subsequent electronic structure studies. The project funded under this award utilized newly implemented fully internally contracted (FIC) nuclear gradients for extended multi-state (XMS) complete active space second-order multireference perturbation theory (CASPT2) to explore geometry changes in first-row transition metal coordination complexes. At the start of the project, only one full geometry optimization using FIC-CASPT2 analytical gradients had been reported (J. Chem. Theory Comput., 2016, 12 (8), 3781), and many open-questions regarding the performance and achievable accuracy in applying such computations to larger complexes persisted. Key deliverables in the report include the implementation of the numerical Hessian and subsequent vibrational analysis in the BAGEL program package. This includes both full Hessian vibrational analysis (FHVA) and a partial Hessian vibrational analysis (PHVA). The latter of which allows us to reduce the significant cost of computing the vibrational frequencies in a large molecule by focusing on the modes of highest interest. The gradient code also proved useful in the development of an approach to improve the Hubbard-U correction by removing bias in predicting spin-splitting in Fe(II) complexes via plane-wave density functional theory (DFT). Finally, we showed for three families of complexes (spin-crossover (SCO) complexes, metallocorroles complexes, and those with metal-metal bonds) that CASPT2 can result in good molecular geometries and established best practices in undertaking this work. We are now using this approach in collaborative projects with experimental groups, which would not have been possible at the start of this project. A common theme also arose through this work that static and dynamic correlation must be recovered in a balanced way to yield quantitative results. By systematically varying how both effects were included, we were able to make key chemical insights. This work supported the training of two postdoctoral scholars, one graduate student, and one undergraduate student in applying high-level wavefunction based methods to challenging systems.

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

Tracking Thermal Transport in Colloidal Quantum Dot Films Using in Situ Time-Resolved X-ray Diffraction

Colloidal quantum dots (QDs) and their thin films are increasingly used in electronic and photonic devices replacing traditional bulk semiconductors. However, thermal properties of QDs remain underexplored relative to device development efforts. This study shows the use of time-resolved X-ray diffraction as a contact-free method to probe the thermal response of QDs in environments representative of the active layer in QD optoelectronic devices, providing in situ insights for future thermal management strategies. Through the extraction of Debye–Waller factors on a subnanosecond time scale, we directly capture the heating and cooling of core/shell CdSe/CdS QDs following pulsed optical excitation. In a QD thin film that actively provides optical gain, the thermal conductivity is found to be as low as 0.13 W m –1 K –1 , because of the poor heat flow between close-packed QD solids. For QDs dispersed in liquids, interfacial thermal conductance dominates thermal relaxation, with a conductance of 15 MW m –2 K –1 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Au as a Surrogate for F: The Case of UAu 6 vs UF 6

Here, anion photoelectron spectroscopy and first-principles quantum chemistry are used to demonstrate to what degree Au can act as a surrogate for F in UF6 and its anion. Unlike UF6, UAu6 exhibits strong ligand–ligand, i.e., Au–Au, interactions, resulting in three low-lying isomers, two of which are three-dimensional while the third isomer has a ring-like quasi two-dimensional structure. Additionally, all the UAu6 isomers have open-shell electrons, which in nearly all cases are localized on the central U atom. As a result, the adiabatic electron affinity and vertical detachment energy are measured to be 3.05 ± 0.05 and 3.28 ± 0.05 eV, respectively, and are in very good agreement with calculations.

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

Grand canonical inverse design of multicomponent colloidal crystals

Inverse design methods are powerful computational approaches for creating colloidal systems which self-assemble into a target morphology by reverse engineering the Hamiltonian of the system. Despite this, these optimization procedures tend to yield Hamiltonians which are too complex to be experimentally realized. An alternative route to complex structures involves the use of several different components, however, conventional inverse design methods do not explicitly account for the possibility of phase separation into compositionally distinct structures. We present an inverse design scheme for multicomponent colloidal systems by combining active learning with a method to directly compute their ground state phase diagrams. This explicitly accounts for phase separation and can locate stable regions of Hamiltonian parameter space which grid-based surveys are prone to miss. Using this we design low-density, binary structures with Lennard-Jones-like pairwise interactions that are simpler than in the single component case and potentially realizable in an experimental setting. This reinforces the concept that ground states of simple, multicomponent systems might be rich with previously unappreciated diversity, enabling the assembly of non-trivial structures with only few simple components instead of a single complex one.

36 MATERIALS SCIENCE↗