Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “liquid dynamics”

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 811 records · Page 45

Experimental Investigation of Mercury's Magma Ocean Viscosity: Implications for the Formation of Mercury's Cumulate Mantle, Its Subsequent Dynamic Evolution, and Crustal Petrogenesis

Mercury has a compositionally heterogeneous surface that was produced by different periods of igneous activity during Mercury’s history, perhaps suggesting heterogeneous mantle sources. Furthermore, understanding the structure of Mercury’s mantle formed during the planet’s magma ocean stage could help in developing a petrologic model for Mercury, and thus, its dynamic history in the context of crustal petrogenesis. We present results of falling sphere viscometry experiments on late stage Mercurian magma ocean analogue compositions. Owing to the presence of sulfur on the surface of Mercury, two compositions were tested, one with sulfur and one without. The liquids have viscosities of 0.6-3.9 (sulfur-bearing) and 0.6-10.9 Pa·s (sulfur-free), similar to an andesite at the tested experimental conditions (1.4-6.2 GPa, 1600-2000°C). We present new viscosity models that enable extrapolation beyond the experimental conditions and evaluate grain growth and the potential for crystal entrainment in a cooling, convecting magma ocean. We consider scenarios with and without a graphite flotation crust, which suggests two possible endmember outcomes for Mercury’s mantle structure. With a graphite flotation crust, crystallization of the mantle would be fractional with negatively buoyant minerals sinking to form a stratified cumulate pile according to the crystallization sequence. Without a flotation crust, crystals would largely remain entrained in the convecting liquid during solidification, producing a homogeneous mantle. In the context of these endmember models, the chemically heterogeneous surface could result from dynamical stirring or mixing of a mantle that was initially compositionally stratified, or from different extents of melting of a homogeneous mantle.

58 GEOSCIENCES↗

Unraveling Anisotropic and Pulsating Etching of ZnO Nanorods in Hydrochloric Acid via Correlative Electron Microscopy

Despite much technical progress achieved so far, the exact surface and shape evolution during wet chemical etching is less unraveled, especially in ionically bonded ceramics. Herein, by using in situ liquid cell transmission electron microscopy, a repeated two-stage anisotropic and pulsating periodic etching dynamic is discovered during the pencil shape evolution of a single crystal ZnO nanorod in aqueous hydrochloric acid. Specifically, the nanopencil tip shrinks at a slower rate along [0001̅] than that along the $\langle$101̅0$\rangle$ directions, resulting in a sharper ZnO pencil tip. Afterward, rapid tip dissolution happens due to accelerated etching rates along various crystal directions. Concurrently, the vicinal base region of the original nanopencil tip emerges as a new tip followed by the repeated sequence of tip shrinking and removal. The high-index surfaces, such as {101̅m} (m = 0, 1, 2, or 3) and {21̅1̅n} (n = 0, 1, 2, or 3), are found to preferentially expose in different ratios. Our 3D electron tomography, convergent beam electron diffraction, middle-angle bright-field STEM, and XPS results indicate the dissociative Cl – species were bound to the Zn-terminated tip surfaces. Furthermore, DFT calculation suggests the preferential Cl – passivation over the {101̅1} and (0001) surfaces of lower energy than others, leading to preferential surface exposures and the oscillatory variation of different facet etching rates. The boosted reactivity due to high-index nanoscale surface exposures is confirmed by comparatively enhanced chemical sensing and CO 2 hydrogenation activity. In conclusion, these findings provide an in-depth understanding of anisotropic wet chemical etching of ionic nanocrystals and offer a design strategy for advanced functional materials.

36 MATERIALS SCIENCE↗

Using Machine Learning to Greatly Accelerate Path Integral Ab Initio Molecular Dynamics

Ab initio molecular dynamics (AIMD) has become one of the most popular and robust approaches for modeling complicated chemical, liquid, and material systems. However, the formidable computational cost often limits its widespread application in simulations of the largest-scale systems. Furthermore, the situation becomes even more severe in cases where the hydrogen nuclei may be better described as quantized particles using a path integral representation. Here, we present a computational approach that combines machine learning with recent advances in path integral contraction schemes, and we achieve a 2 orders of magnitude acceleration over direct path integral AIMD simulation while at the same time maintaining its accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chapter 14: Dynamic Flux Analysis: An Experimental Approach of Fluxomics

Metabolic flux analysis represents an essential perspective to understand cellular physiology and offers quantitative information to guide pathway engineering. A valuable approach for experimental elucidation of metabolic flux is dynamic flux analysis, which estimates the relative or absolute flow rates through a series of metabolic intermediates in a given pathway. It is based on kinetic isotope labeling experiments, liquid chromatography-mass spectrometry (LC-MS), and computational analysis that relate kinetic isotope trajectories of metabolites to pathway activity. Herein, we illustrate the mathematic principles underlying the dynamic flux analysis and mainly focus on describing the experimental procedures for data generation. This protocol is exemplified using cyanobacterial metabolism as an example, for which reliable labeling data for central carbon metabolites can be acquired quantitatively. This protocol is applicable to other microbial systems as well and can be readily adapted to address different metabolic processes.

BASIC BIOLOGICAL SCIENCES↗

Extending JuTrack’s capabilities to the FRIB accelerator to enhance online modeling

JuTrack is a Julia-based accelerator modeling and tracking package that utilizes compiler-level automatic differentiation (AD) to enable fast and accurate derivative calculations. While JuTrack provides a solid foundation for beam dynamics simulations, its capabilities must be extended to support the Facility for Rare Isotopes (FRIB) linac. This includes modeling heavy-ion linac accelerator components such as the liquid-lithium charge stripper, which facilitates efficient acceleration by remove electrons from heavy isotopes, and incorporating multi-charge state acceleration tracking, which allows for charge-dependent beam dynamics. These extensions address challenges such as the beam matching and optimization of multi charge state through various accelerating structures and beam-material interaction modeling while maintaining the auto differentiation capability. This work focuses on adapting JuTrack to incorporate these elements, enhancing its online modeling abilities. We present modifications to JuTrack’s framework and demonstrate their performance in FRIB simulations.

Accelerator Physics↗

The total neutron cross section of liquid and solid ammonia

Ammonia is a material of interest for future neutron moderators at high-power sources due to its high hydrogen density, low melting point, and resistance to polymerization in an intense radiation field. Its performance in such applications cannot currently be calculated due to the absence of suitable computer models for the interaction of neutrons with ammonia under relevant conditions. In an effort to develop suitable scattering kernels for computer simulations of moderator performance, we have conducted a series of Density Functional Theory and Molecular Dynamics calculations of the molecular-level thermal properties of ammonia at various temperatures within both the solid and liquid phases. In this paper, we compare computer calculations for the energy-dependent total neutron cross section of ammonia, based on these models, to experimental measurements of those cross sections at temperatures of 221 K, 180 K, and 35 K. The experimental data were collected over an energy range from 0.1 meV to 10 eV using time-of-flight techniques at the Low Energy Neutron Source (LENS) facility at Indiana University. This comparison provides a first validation in the development of thermal scattering libraries for Monte Carlo source design simulations based on liquid and solid ammonia. In conclusion, we also provide some insights into where additional development of tools for creating such models may be needed.

Ammonia↗

Interfacial Issues and Modification of Solid Electrolyte Interphase for Li Metal Anode in Liquid and Solid Electrolytes

Abstract The high energy density required for the next generation of lithium batteries will likely be enabled by a shift toward lithium metal anode from the conventional intercalation‐based anode such as graphite. However, several critical challenges for Li metal originate from its highly reactive nature and the hostless reaction of deposition and stripping impede the practical use of Li metal as an anode. The role of the solid electrolyte interphase (SEI) is very important for the Li metal anode where the SEI must protect the dynamically changing surface of the Li metal. Since the SEI‐generating reaction mechanisms for the two different electrolyte systems, liquid and solid, are considerably different, the SEI layers formed between the Li metal and the electrolytes in the two electrolyte systems have substantially different properties, causing different interfacial issues. Inhibition of the interfacial problems requires different strategies to reinforce the SEI layer for each of the electrolyte systems. However, the differences in the two electrolyte systems have not been clearly compared in the prior literature. In this report, the interfacial issues for the two different electrolyte systems are compared and different strategies for SEI modification are provided to overcome the issues.

25 ENERGY STORAGE↗

Characterization of plasma in RF jet interacting with water: Thomson scattering versus spectral line broadening

Here, we carry out a detailed characterization of an Ar radiofrequency plasma jet interacting with liquid. The focus of the paper is measurement of the electron density by examining the Stark broadening of hydrogen Balmer α line and Thomson scattering (TS). Supporting diagnostics are done to investigate the channel evolution and movement, and gas temperature. The comparison of the two approaches shows the caveats and the advantages that should gain attention in future applications. In particular, the plasma channel dynamics have a significant impact on the TS signal and may result in physical phenomena being missed due to channel dislocations. The conclusions of the work elucidate the pitfalls for interpreting the results of TS when the discharge is a dynamic filamentary plasma. This work establishes the temporal evolution of the plasma and the gas parameters in a plasma–liquid system employed for investigation of plasma-induced electrodeless electrochemistry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microphysics of liquid water in sub-10 nm ultrafine aerosol particles

Abstract. Ultrafine aerosol particles with sizes smaller than 50 nm have been shown in recent studies to serve as a large source of cloud condensation nuclei that can promote additional cloud droplet formation under supersaturation conditions. Knowledge of the microphysics of liquid water in these droplets remains limited, particularly in the sub-10 nm particle size range, due to experimental and theoretical challenges associated with the complexity of aerosol components and the small length scales of interest (e.g., difficulty of precisely sampling the liquid–air interface, questionable validity of mean-field theoretical representations). Here, we carried out molecular dynamics (MD) simulations of aerosol particles with diameters between 1 and 10 nm and characterized atomistic-level structure and water dynamics in well-mixed and phase-separated systems with different particle sizes, NaCl salinities, and pimelic acid (PML) organic surface loadings as a function of distance from the time-averaged Gibbs dividing interface or instantaneous water–air interface. We define a sphericity factor (ϕ) that can shed light on the phase-mixing state of nanodroplets, and we reveal an unexpected dependence of mixing state on droplet size. Our results also evidence an ion concentration enhancement in ultrafine aerosols, which should modulate salt nucleation kinetics in sub-10 nm droplets, and provide detailed characterization of the influence of droplet size on surface tension and on water self-diffusivity near the interface. Analysis of water evaporation free energy and water activity demonstrates the validity of the Kelvin equation and Köhler theory at droplet sizes larger than 4 nm under moderate salinities and organic loadings and the need for further extension to account for ion concentration enhancement in sub-10 nm aerosols, droplet-size-dependent phase separation effects, and a sharp decrease in the cohesiveness of liquid water in sub-4 nm droplets. Finally, we show that an idealized fractional surface coating factor (fs) can be used to categorize and reconcile water accommodation coefficients (α*) observed in MD simulations and experimental results in the presence of organic coatings, and we resolve the droplet size dependence of α*.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

In Situ Liquid Electron Microscope Cells Strongly Attenuate Electrochemical Behavior

In situ electrochemical Scanning Transmission Electron Microscopy (ec-S/TEM) probes the dynamics of electrochemical processes in real time and at high spatial resolution. This method utilizes microfabricated electrochemical cells, with enclosed dimensions of 800 μm wide, 800 μm long, and 1–2 μm tall, with electron-transparent silicon nitride windows. Working, counter, and reference electrodes in the configuration studied are deposited Pt with a planar surface area of about 1 μm 2 and a height of 50 nm. The cell confines the electrolyte within a 1 μm tall channel, and the microfabricated electrodes have a non-standard geometry, are coplanar, and are not front-facing. As such, standard assumptions of electrochemical experiments do not apply. COMSOL modeling was used to compare and contrast differences in electrochemical behavior between conventional experimental setups and the in situ ec-S/TEM cells. Cell height strongly affected voltammogram measurements: peak magnitudes and shifts correlate directly to cell height. Reactant concentration above the electrode quickly decreases during deposition, leading to a planar diffusion-dominated regime. The concentration gradient prioritizes particle growth along the thin edges of the electrode and not on the planar part of the electrode. Experiments done in the ec-S/TEM ex situ of a TEM verify the model’s accuracy.

Wittman, Reed M. [Oak Ridge National Laboratory (O↗

Facile Measurement of the Rotation of a Single Optically Trapped Nanoparticle Using the Diagonal Ratio of a Quadrant Photodiode

Optical tweezers are a powerful tool for exploring physical properties of particles in various environments through analysis of their dynamics in a trapping potential. Analysis of the trapped particle's position is often done using interferometric back-focal-plane detection microscopy in which interference between the trapping laser and the light interacting with the trapped particle is projected on a quadrant photo-diode (QPD). The QPD measurements are almost universally the voltage ratios between the transverse components of the detector (left/right; up/down) that are calibrated to ascertain particle positions and orientations (i.e., rotation). In this paper, we demonstrate how measuring the diagonal ratio of the QPD allows tracking the spinning and orbital motion of optically trapped objects with enhanced sensitivity compared to the conventional transverse ratio and reveals additional information about particle dynamics. The diagonal ratio "balances" the signal from opposite sides of the QPD, allowing common mode noise reduction and increased measurement sensitivity of the rotating nanospheres at frequencies from below 200 Hz to above 15 kHz. We show that the variability in particle rotation rates is due to slight differences in the aspect ratio of the nanoparticles. We also measure in situ a gradual acceleration in the spinning frequency of trapped nanoparticles due to photothermal shedding of surface ligands and the resulting decrease of rotational friction in solution. The accelerating particles spin at rates exceeding 30 kHz and are eventually "optically printed" onto the glass coverslip on top of the liquid cell enclosure. The diagonal ratio is a simple, readily available, and powerful enabler for measuring rotational dynamics with a 5-fold sensitivity increase over conventional QPD measurements. This affords opportunities for controlling the orientational dynamics of trapped nanoparticles, which is crucial for various applications, including nanoscale mechanical motors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling graphene sheet growth and dynamical matrix calculations using molecular dynamics

Molecular dynamics (MD) has been an incredibly useful tool to model physical processes that were synthesized experimentally but not fully understood. MD, through the use of semi-empirical inter-atomic potentials, has allowed understanding of different physical processes in materials science. Yet as well as providing useful insights into materials science, molecular dynamics has a wider range of usability. In this report, I will be detailing how MD can be used to study graphene formation from a carbon liquid which requires high temperatures and pressures. Beyond this, I will describe the usefulness of MD for understanding the physics for phonon transport quantum sensors. To do this, MD was employed to determine the dynamical matrix by treating atoms as coupled oscillators. An accurate understanding of the dynamical matrix of a system is required to calculate the non-equilibrium Green’s function used to describe the phonon transport within phonon wave-guides. I found that, across multiple pressures and temperatures, randomly placed carbon atoms will show evidence of pent-first formation with semi-empirical models. Density functional theory (DFT), on the other hand, was too computationally expensive to use for full scale MD simulations, but we have the possibility of training a machine learned interatomic potential to approximate DFT for carbon in the environments being studied for pent-first graphene sheet formation.

36 MATERIALS SCIENCE↗

Examination of the cerium α-ε phase transition under dynamic loading with x-ray diffraction

Recent examination of the cerium Hugoniot with pyrometry and x-ray diffraction (XRD) has revealed a narrow solid-liquid two-phase region. It has been suggested that nonequilibrium melting may be occurring along the Ce Hugoniot, with either melt kinetics or a sluggish α-ε transition impeding the transition. In particular, the kinetics of the α-ε is unknown and the location of the phase boundary is in dispute. Static measurements suggest a nearly vertical phase boundary that intersects the Hugoniot at 6–7 GPa. This lies in direct conflict with dynamic measurements along the Hugoniot observing α-Ce through incipient melt. This work presents dynamic experiments using XRD to examine the behavior of the α-ε phase transition. The results show that the α-ε phase transition occurs through a tetragonal distortion, with the transition beginning at temperatures below the solid Hugoniot. Following the initial deviation from an ideal fcc structure, the c/a ratio is found to gradually increase with no steady value observed in the ε phase within the range of these experiments (below 17 GPa). Multiple diffraction patterns captured during the peak stress state show no significant change in c/a ratio prior to uniaxial release, upon which Ce reverts to an fcc structure. Furthermore, the results indicate that the α-ε transition occurs rapidly, both on loading and release. An examination of the c/a ratio with increasing temperatures suggests 11.5 GPa as a lower bound for the location of the α-ε-liquid triple point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Computational fluid dynamics modeling and analysis of silica nanoparticle synthesis in a flame spray pyrolysis reactor

Flame Spray Pyrolysis (FSP) is a method for large-scale production of nanoparticles and nanoscale powders employed in a wide range of industrial applications. Particle size and morphology are complex functions of the physicochemical phenomena occurring in the FSP reactor. An extensive study of FSP-related phenomena can be utilized to develop effective strategies for achieving desired particle size/morphology and scaling up the overall yield of an FSP system. In this work, a computational fluid dynamics (CFD) model of an FSP reactor is developed to simulate the coupling of key phenomena involved in the particle synthesis process: liquid spray breakup and evaporation, mixing, combustion, and particle formation/growth of silica nanoparticles. Herein, the particle sizes and their distributions from the CFD simulations are validated against experimental data. Subsequently, the simulations are utilized to investigate the impact of process parameters on the resultant flame dynamics and particle growth. Firstly, the CFD results show that the particle sizes are strongly correlated with the precursor concentration in the solvent. At lower precursor concentrations, the spread of the distribution is relatively insensitive to the value of the concentration. At higher concentrations, the spread is higher as the collision probability between particles is higher. Secondly, increasing the pilot flow rate increases the length of the pilot flames impacting the local ignition location of the spray flame. Lastly, it is shown that the dispersion gas flow rate strongly influences the spray flame shape. This shape can be used for control of particle growth as it helps determine the regions of high temperature and the residence time of the particles in the high temperature region enabling the design and process optimization of the FSP reactor.

42 ENGINEERING↗

Interfacial Hydrogen-Bond Dynamics in Transition Metal Compounds

Understanding how water behaves when confined within atomic layers of active transition-metal carbides, nitrides, and carbonitrides is essential for uncovering the fundamental principles needed to engineer solid–liquid interfaces at the atomic scale. Yet, how lattice element chemistry and surface termination groups collectively regulate the structure and mobility of such interlayer water remains poorly understood. Here, we present a composition-controlled investigation of interlayer water dynamics in layered transition-metal nitride, carbide, and carbonitride systems using a systematic integration of quasi-elastic neutron scattering (QENS), ab initio molecular dynamics (AIMD) simulations, and density functional theory (DFT) calculations. QENS measurements show that nitride-rich systems host mobile, translationally diffusing water with thermally activated self-diffusion coefficients on the order of 10 –10 m 2 s –1 , whereas mixed C/N lattices confine water to localized, nontranslational motion that is insensitive to temperature. AIMD and DFT reveal that lattice C/N chemistry and surface functional group composition reshape the first hydration layer by modulating the surface electronic structure and termination-dependent hydrogen-bond networks, leading to pronounced differences in water ordering and thermal resilience. On the other hand, fully carbide systems exhibit intermediate behavior, highlighting that water mobility is not primarily controlled by the hydration level alone but by the coupling between lattice composition and surface chemistry. Overall, this study establishes how surface chemistry and lattice composition jointly control interfacial hydrogen bond dynamics, offering a mechanistic framework for designing transition-metal layered materials with tailored interfacial transport properties.

Hydration↗

CO 2 adsorption mechanisms at the ZIF-8 interface in a Type 3 porous liquid

Porous liquids (PLs) are an attractive material for gas separation and carbon sequestration due to their permanent internal porosity and high adsorption capacity. PLs that contain zeolitic imidazole frameworks (ZIFs), such as ZIF-8, form PLs through exclusion of aqueous solvents from the framework pore due to its hydrophobicity. The gas adsorption sites in ZIF-8 based PLs are historically unknown; gas molecules could be captured in the ZIF-8 pore or adsorb at the ZIF-8 interface. To address this question, ab initio molecular dynamics was used to predict CO 2 binding sites in a PL composed of a ZIF-8 particle solvated in a water, ethylene glycol, and 2-methylimidazole solvent system. Further, the results show that CO 2 energetically prefers to reside inside the ZIF-8 pore aperture due to strong van der Waals interactions with the terminal imidazoles. However, the CO 2 binding site can be blocked by larger solvent molecules that have greater adsorption interactions. CO 2 molecules were unable to diffuse into the ZIF-8 pore, with CO 2 adsorption occurring due to binding with the ZIF-8 surface. Therefore, future design of ZIF-based PLs for enhanced CO 2 adsorption should be based on the strength of gas binding at the solvated particle surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of Dynamical Peculiarities in Nanoalloys at Subsystems Level: Dynamical Degrees of Freedom, Temperature Differences, and the Chameleon Effect

Abstract A novel analysis of the dynamical behavior of nanoalloy systems, as represented by model Ni/Al 13‐atom clusters, over a broad range of energies that cover the stage‐wise transition of the systems from their solid‐like to liquid‐like state is presented. Conceptually, the analysis is rooted in partitioning the systems into judiciously chosen subsystems and characterizing the latter in terms of subsystem‐specific dynamical descriptors that include dynamical degrees of freedom , root‐mean‐square bond‐length fluctuation , and element‐specific subsystem temperature . The analysis reveals a host of intriguing new peculiarities in the dynamical behavior of the Ni/Al 13‐mers, among which are what we call the chameleon effect and the difference in the temperatures of the Ni and Al subsystems at high energies, a difference that strongly depends on the cluster composition and also changes with energy. These do not have an analog in pure Ni 13 and Al 13 and are explained in terms of the coupled effects of the difference between the masses of the Ni and Al atoms (the mass effect) and of the difference in the anharmonicity of the overall interaction potential as experienced by the Ni and Al subsystems of the clusters (the potential effect).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-resolution infrared spectroscopy of jet cooled cyclobutyl in the α-CH stretch region: large-amplitude puckering dynamics in a 4-membered ring radical

We report gas-phase cyclobutyl radical (c-C 4 H 7 ) is generated at a rotational temperature of T rot = 26(1) K in a slit-jet discharge mixture of 70% Ne/30% He and 0.5–0.6% cyclobromobutane (c-C 4 H 7 Br). A fully rovibrationally resolved absorption spectrum of the α-CH stretch fundamental band between 3062.9 cm −1 to 3075.7 cm −1 is obtained and analyzed, yielding first precision structural and dynamical information for this novel radical species. The α-CH stretch band origin is determined to be 3068.7887(4) cm −1 , which implies only a modest (≈0.8 cm −1 ) blue shift from rotationally unresolved infrared spectroscopic studies of cyclobutyl radicals in liquid He droplets [ A. R. Brown, P. R. Franke and G. E. Douberly, J. Phys. Chem. A , 2017, 121 , 7576–7587]. Of particular dynamical interest, a one-dimensional potential energy surface with respect to the ring puckering coordinate is computed at CCSD(T)/ANO2 level of theory and reveals a double minimum C s puckered geometry, separated by an exceedingly shallow planar C 2v transition state barrier (E barr ≈ 1 cm −1 ). Numerical solutions on this double minimum potential yield a zero-point energy for the ground state (E zero-point ≈ 27 cm −1 ) greatly in excess of the interconversion barrier. This is indicative of highly delocalized, large amplitude motion of the four-membered ring structure, for which proper vibrationally averaging of the moment of inertia tensor reproduces the experimentally determined inertial defect remarkably well. Finally, intensity alternation in the experimental spectrum due to nuclear spin statistics upon exchange of three indistinguishable H atom pairs (IH = ½) matches K a + K c = even : odd = 36 : 28 predictions, implying that the unpaired electron in the radical center lies in an out-of-plane p π orbital. Thus, high-resolution infrared spectroscopy provides first experimental confirmation of a shallow double minimum ring puckering potential with a highly delocalized ground state wave function peaked at a planar C 2v transition state geometry consistent with a cyclobutyl π radical.

Chemistry↗