Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecule”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

Morphology of a columnar stack of coronene molecules encapsulated in a single-walled carbon nanotube

The morphology of stacked coronene molecules encapsulated in a single-walled carbon nanotube (SWCNT) is investigated using atomistic simulation. First, the minimum energy configuration of coronene molecules in a SWCNT is sought by means of conjugate gradient (CG) minimization. Secondly, encapsulation of coronene molecules into a SWCNT existing in a coronene atmosphere is simulated by means of molecular dynamics (MD). In both of the simulations, the diameter of the SWCNT ranges from 1.35 to 1.69 nm, and the final configurations of coronene molecules within a SWCNT are examined. In a thin SWCNT, coronene molecules tilt against the radial direction of the SWCNT and slide relative to each other, whereas in a thick SWCNT, they do not tilt but rotate relative to each other. In a SWCNT of the intermediate diameter, they tilt, slide, and rotate. For the SWCNT diameter less than or equal to 1.52 nm, the mean tilt angle of the stacked coronene molecules almost linearly decreases with increasing the diameter, whereas for the diameter above 1.52 nm, it is approximately 0{sup ∘}. To check the validity of the results, the MD simulations are performed changing the density of the coronene atmosphere and the length of the SWCNT; the results prove to be valid. Finally, the effects of temperature on the mean tilt angle and mean intermolecular distance of stacked coronene molecules are examined by a rather simplified simulation, which shows that both of them increase with increasing temperature.

36 MATERIALS SCIENCE↗

Fast optical transport of ultracold molecules over long distances

Abstract Optically trapped laser-cooled polar molecules hold promise for new science and technology in quantum information and quantum simulation. Large numerical aperture optical access and long trap lifetimes are needed for many studies, but these requirements are challenging to achieve in a magneto-optical trap (MOT) vacuum chamber that is connected to a cryogenic buffer gas beam source, as is the case for all molecule laser cooling experiments so far. Long distance transport of molecules greatly eases fulfilling these requirements as molecules are placed into a region separate from the MOT chamber. We realize a fast transport method for ultracold molecules based on an electronically focus-tunable lens combined with an optical lattice. The high transport speed is achieved by the 1D red-detuned optical lattice, which is generated by interference of a focus-tunable laser beam and a focus-fixed laser beam. Efficiency of 48(8)% is realized in the transport of ultracold calcium monofluoride (CaF) molecules over 46 cm distance in 50 ms, with a moderate heating from 32(2) μ K to 53(4) μ K. Positional stability of the molecular cloud allows for stable loading of an optical tweezer array with single molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗

Electrochemically modulated single-molecule localization microscopy for in vitro imaging cytoskeletal protein structures

A new concept of electrochemically modulated single-molecule localization super-resolution imaging is developed. Applications of single-molecule localization super-resolution microscopy have been limited due to insufficient availability of qualified fluorophores with favorable low duty cycles. The key for the new concept is that the “On” state of a redox-active fluorophore with unfavorable high duty cycle could be driven to “Off” state by electrochemical potential modulation and thus become available for single-molecule localization imaging. The new concept was carried out using redox-active cresyl violet with unfavorable high duty cycle as a model fluorophore by synchronizing electrochemical potential scanning with a single-molecule localization microscope. The two cytoskeletal protein structures, the microtubules from porcine brain and the actins from rabbit muscle, were selected as the model target structures for the conceptual imaging in vitro. The super-resolution images of microtubules and actins were obtained from precise single-molecule localizations determined by modulating the On/Off states of single fluorophore molecules on the cytoskeletal proteins via electrochemical potential scanning. Importantly, this method could allow more fluorophores even with unfavorable photophysical properties to become available for a wider and more extensive application of single-molecule localization microscopy.

electrochemical modulation↗

Single-Molecule Chemical Reactions Unveiled in Molecular Junctions

Understanding chemical processes at the single-molecule scale represents the ultimate limit of analytical chemistry. Single-molecule detection techniques allow one to reveal the detailed dynamics and kinetics of a chemical reaction with unprecedented accuracy. It has also enabled the discoveries of new reaction pathways or intermediates/transition states that are inaccessible in conventional ensemble experiments, which is critical to elucidating their intrinsic mechanisms. Thanks to the rapid development of single-molecule junction (SMJ) techniques, detecting chemical reactions via monitoring the electrical current through single molecules has received an increasing amount of attention and has witnessed tremendous advances in recent years. Research efforts in this direction have opened a new route for probing chemical and physical processes with single-molecule precision. This review presents detailed advancements in probing single-molecule chemical reactions using SMJ techniques. We specifically highlight recent progress in investigating electric-field-driven reactions, reaction dynamics and kinetics, host–guest interactions, and redox reactions of different molecular systems. Finally, we discuss the potential of single-molecule detection using SMJs across various future applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Guest Molecule Exchange Kinetics for the 2012 Ignik Sikumi Gas Hydrate Field Trial

A commercially viable technology for producing methane from natural gas hydrate reservoirs remains elusive. Short-term depressurization field tests have demonstrated the potential for producing natural gas via dissociation of the clathrate structure, but the long-term performance of the depressurization technology ultimately requires a heat source to sustain the dissociation. A decade of laboratory experiments and theoretical studies have demonstrated the exchange of pure CO2 and N2-CO2 mixtures with CH4 in sI gas hydrates, yielding critical information about molecular mechanisms, recoveries, and exchange kinetics. Findings indicated the potential for producing natural gas with little to no production of water and rapid exchange kinetics, generating sufficient interest in the guest-molecule exchange technology for a field test. In 2012 the U.S. DOE/NETL, ConocoPhillips Company, and Japan Oil, Gas and Metals National Corporation jointly sponsored the first field trial of injecting a mixture of N2-CO2 into a CH4-hydrate bearing formation beneath the permafrost on the Alaska North Slope. Known as the Ignik Sikumi #1 Gas Hydrate Field Trial, this experiment involved three stages: 1) the injection of a N2-CO2 mixture into a targeted hydrate-bearing layer, 2) a 4-day pressurized soaking period, and 3) a sustained depressurization and fluid production period. Data collected during the three stages of the field trial were made available after an extensive quality check. These data included continuous temperature and pressure logs, injected and recovered fluid compositions and volumes. The Ignik Sikumi #1 data set is extensive, but contains no direct evidence of the guest-molecule exchange process. This investigation is directed at using numerical simulation to provide an interpretation of the collected data. A numerical simulator, STOMP-HYDT-KE, was recently completed that solves conservation equations for energy, water, mobile fluid guest molecules, and hydrate guest molecules, for up to three gas hydrate guest molecules: CH4, CO2, and N2. The independent tracking of mobile fluid and hydrate guest molecules allows for the kinetic exchange of guest molecules between the mobile fluids and hydrate. The particular interest of this numerical investigation is to determine whether kinetic exchange parameters, determined from laboratory-scale experiments, are directly applicable to interpreting the Ignik Sikumi #1 data.

Alaska North Slope↗

Components of the Bond Energy in Polar Diatomic Molecules, Radicals, and Ions Formed by Group-1 and Group-2 Metal Atoms

While many transition metal complexes are known to have high multireference character, the multireference character of main-group closed-shell singlet diatomic molecules like BeF, CaO, and MgO has been less studied. Yet, many group-1 and group-2 diatomic molecules do have multireference character, and they provide informative systems for studying multireference character because they are simpler than transition metal compounds. The goal of the present work is to understand these multireference systems better so that, ultimately, we can apply what we learn to more complicated multireference systems and to the design of new exchange-correlation functionals for treating multireference systems more adequately. Fourteen main-group diatomic molecules and one triatomic molecule (including radicals, cations, and anions, as well as neutral closed-shell species) have been studied for this article. Eight of these molecules contain a group-1 element, and six contain a group-2 element. Seven of these molecules are multireference systems, and eight of them are single-reference systems. Fifty-three exchange-correlation functionals of 11 types [local spin-density approximation (LSDA), generalized gradient approximation (GGA), nonseparable gradient approximation (NGA), global-hybrid GGA, meta-GGA, meta-NGA, global-hybrid meta GGA, range-separated hybrid GGA, range-separated hybrid meta-GGA, range-separated hybrid meta-NGA, and DFT augmented with molecular mechanics damped dispersion (DFT-D)] and the Hartree–Fock method have been applied to calculate the bond distance, bond dissociation energy (BDE), and dipole moment of these molecules. All of the calculations are converged to a stable solution by allowing the symmetry of the Slater determinant to be broken. A reliable functional should not only predict an accurate BDE but also predict accurate components of the BDE, so each bond dissociation energy has been decomposed into ionization potential (IP) of the electropositive element, electron affinity of the electronegative bonding partner (EA), atomic excitation energy (EE) to prepare the valence states of the interacting partners, and interaction energy (IE) of the valence-prepared states. Adding Hartree–Fock exchange helps to obtain better results for atomic excitation energy, and this leads to improvements in getting the right answer for the right reason. The following functionals are singled out for reasonably good performance on all three of bond distance, BDE, and dipole moment: B97-1, B97-3, MPW1B95, M05, M06, M06-2X, M08-SO, N12-SX, O3LYP, TPSS, τ-HCTHhyb, and GAM; all but two (TPSS and GAM) of these functionals are hybrid functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Growth of Ammonium Bisulfate Clusters by Adsorption of Oxygenated Organic Molecules

Quantum chemical calculations were employed to model the interactions of the [(NH 4 + ) 4 (HSO 4 - ) 4 ] ammonium bisulfate cluster with one or more molecular products of monoterpene oxidation. A strong interaction was found between the bisulfate ion of the cluster and a carboxylic acid, aldehyde or ketone functionality of the organic molecule. Free energies of adsorption for carboxylic acids were in the -70 to -73 kJ/mol range, while those for aldehydes and ketones were in the -46 to -50 kJ/mol range. These values suggest that a small ambient ammonium bisulfate cluster, such as the [(NH 4 + ) 4 (SO 4 - ) 4 ] cluster, is able to adsorb an oxygenated organic molecule. Although adsorption of the first molecule is highly favorable, adsorption of subsequent molecules is not, suggesting that sustained uptake of organic molecules does not occur, and thus is not a pathway for continuing growth of the cluster. This result is consistent with ambient measurements showing that particles below ~1 nm grow slowly, while those above 1 nm grow at an increasing rate presumably due to a lower surface energy barrier enabling the uptake of organic molecules. This work provides insight into the molecular level interactions which affect sustained cluster growth by uptake of organic molecules.

54 ENVIRONMENTAL SCIENCES↗

Effects of Functional Groups in Redox-Active Organic Molecules: A High-Throughput Screening Approach

Nonaqueous redox flow batteries have attracted recent attention with their potential for high electrochemical storage capacity, with organic electrolytes serving as solvents with a wide electrochemical stability window. Organic molecules can also serve as electroactive species, where molecules with low reduction potentials or high oxidation potentials can provide substantial chemical energy. To identify promising electrolytes in a vast chemical space, high-throughput screening (HTS) of candidate molecules plays an important role, where HTS is used to calculate properties of thousands of molecules and identify a few organic molecules worthy of further attention in battery research. Here, in this work, we present reduction and oxidation potentials obtained from HTS of 4178 molecules. The molecules are composed of base groups of five- or six-membered rings with one or two functional groups attached, with the set of possible functional groups including both electron-withdrawing and electron-donating groups. In addition to observing the trends in potentials that result from differences in organic base groups and functional groups, we analyze the effects of molecular characteristics such as multiple bonds, Hammett parameters, and functional group position. In conclusion, this work provides useful guidance in determining how the identities of the base groups and functional groups are correlated with desirable reduction and oxidation potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Active learning enables generation of molecules that advance the known Pareto front

Although generative models hold promise for discovering molecules with optimized desired properties, they often fail to suggest synthesizable molecules that improve upon the properties of the structures represented in the training distribution. We find that this limitation arises not only from the molecule generation process itself, but also from the poor generalization capabilities of molecular property predictors. We address this challenge by creating a closed-loop molecule generation pipeline with iterative retraining on new quantum chemical simulation data. Compared against static, single-pass generative modeling approaches, only our closed-loop iterative workflow generates molecules with properties extending beyond the training distribution (up to 0.44 standard deviations beyond the original range) and achieves a 79% improvement in out-of-distribution molecule classification accuracy. Furthermore, by conditioning molecular generation on thermodynamic stability data obtained during the iterative loop, the proportion of stable and hence potentially synthesizable molecules generated is 3.5x higher than the next-best model.

Chemistry↗

Selective triplet exciton formation in a single molecule

The formation of excitons in organic molecules by charge injection is an essential process in organic light-emitting diodes (OLEDs). According to a simple model based on spin statistics, the injected charges form spin-singlet (S1) excitons and spin-triplet (T 1 ) excitons in a 1:3 ratio. After the first report of a highly efficient OLED based on phosphorescence, which is produced by the decay of T 1 excitons, more effective use of these excitons has been the primary strategy for increasing the energy efficiency of OLEDs. Another route to improving OLED energy efficiency is reduction of the operating voltage. Because T 1 excitons have lower energy than S 1 excitons (owing to the exchange interaction), use of the energy difference could—in principle—enable exclusive production of T 1 excitons at low OLED operating voltages. However, a way to achieve such selective and direct formation of these excitons has not yet been established. Here we report a single-molecule investigation of electroluminescence using a scanning tunnelling microscope and demonstrate a simple method of selective formation of T 1 excitons that utilizes a charged molecule. A 3,4,9,10-perylenetetracarboxylicdianhydride (PTCDA) molecule adsorbed on a three-monolayer NaCl film atop Ag(111) shows both phosphorescence and fluorescence signals at high applied voltage. In contrast, only phosphorescence occurs at low applied voltage, indicating selective formation of T 1 excitons without creating their S 1 counterparts. The bias voltage dependence of the phosphorescence, combined with differential conductance measurements, reveals that spin-selective electron removal from a negatively charged PTCDA molecule is the dominant formation mechanism of T 1 excitons in this system, which can be explained by considering the exchange interaction in the charged molecule. Our findings show that the electron transport process accompanying exciton formation can be controlled by manipulating an electron spin inside a molecule. We anticipate that designing a device taking into account the exchange interaction could realize an OLED with a lower operating voltage.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Characterization of Fe 2+ ions in Fe,H/SSZ-13 zeolites: FTIR spectroscopy of CO and NO probe molecules

The IR spectra of adsorbed CO and NO probe molecules were used to characterize the coordination chemistry of Fe2+ ions in solution ion exchanged Fe,H/SSZ-13 zeolites. The effects of Fe ion exchange levels, as well the sample pre-treatment conditions on the adsorption of these probe molecules were investigated. The ion exchange levels (in the range of the study) did not affect significantly the IR spectra of either probe molecule, the IR features and their intensity ratios were very similar. Experiments with both probe molecules substantiated the presence of two distinct types of Fe2+ ions in cationic positions. We assign these two Fe2+ ions to two distinct cationic positions: Fe2+ in 6R and 8R positions. NO initially adsorbs preferentially onto Fe2+ sites in the 6R position, and then populates sites in the 8R. Fe2+ ions in the 8R positions require the interaction of more than one NO molecule to move them out from their adsorbate-free cationic positions. As soon as they move from their stable positions, they are able to bind to multiple NO molecules, and form mostly tri-nitrosyls. These tri-nitrosyls, however, are only stable in the presence of gas phase NO; under dynamic vacuum they lose one of the NO molecules from their coordination sphere and form stable di-nitrosyls. The adsorption of CO is much weaker on Fe2+ sites than that of NO, and requires cryogenic sample temperatures to initiate CO adsorption. Under the conditions applied in this study, only mono-carbonyl formation was observed. Reduction in H2 at 773 K increased the number of Fe2+ adsorption sites, primarily in the 8R locations. Oxidation by N2O, on the other hand, selectively reduced the adsorption of both CO and NO on the Fe2+ sites in 8R positions. Adsorbed oxygen left behind from the decomposition of N2O at 573 K readily reacted with CO to produce CO2 even at 150 K. The authors gratefully acknowledge the US Department of Energy (DOE), Energy Efficiency and Renewable Energy, Vehicle Technologies Program for the support of this work. The research described in this paper was performed in the Environmental Molecular Sciences Laboratory (EMSL), a national scientific user facility sponsored by the DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory (PNNL). PNNL is operated for the US DOE by Battelle.

coordination chemistry of Fe2+ ions↗

“Freedom of design” in chemical compound space: towards rational in silico design of molecules with targeted quantum-mechanical properties

The rational design of molecules with targeted quantum-mechanical (QM) properties requires an advanced understanding of the structure–property/property–property relationships (SPR/PPR) that exist across chemical compound space (CCS). In this work, we analyze these fundamental relationships in the sector of CCS spanned by small (primarily organic) molecules using the recently developed QM7-X dataset, a systematic, extensive, and tightly converged collection of 42 QM properties corresponding to ≈4.2M equilibrium and non-equilibrium molecular structures containing up to seven heavy/non-hydrogen atoms (including C, N, O, S, and Cl). By characterizing and enumerating progressively more complex manifolds of molecular property space—the corresponding high-dimensional space defined by the properties of each molecule in this sector of CCS—our analysis reveals that one has a substantial degree of flexibility or “freedom of design” when searching for a single molecule with a desired pair of properties or a set of distinct molecules sharing an array of properties. To explore how this intrinsic flexibility manifests in the molecular design process, we used multi-objective optimization to search for molecules with simultaneously large polarizabilities and HOMO–LUMO gaps; analysis of the resulting Pareto fronts identified non-trivial paths through CCS consisting of sequential structural and/or compositional changes that yield molecules with optimal combinations of these properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discovery of hydrogen storage molecules using large language models and machine learning

Accelerating the discovery of new molecules with targeted properties is a central challenge in molecular design. In this contribution, we present an AI-driven molecular discovery framework that integrates Large Language Models (LLMs) for generative molecular design with Machine Learning (ML)-based screening to identify novel Liquid Organic Hydrogen Carrier (LOHC) candidates. Using the developed framework, LOHC molecules were systematically generated, evaluated, and refined iteratively, combining LLM-guided molecular generation and ML-predicted hydrogenation enthalpies (Δ H ), under physicochemical property constraints such as optimal melting points (MP), desired hydrogen storage capacity (wt% H 2 ), and synthetic accessibility (SA) scores. This approach enabled the discovery of 42 new LOHC candidates in two distinct campaigns, one seeded with experimentally known and another with previously computationally identified LOHCs, respectively. Although we began with different numbers of starting molecules (31 vs . 7 seed molecules), both runs yielded a comparable number of viable candidates, suggesting an influence of chemically intuitive seed molecule selection for success. Selected LOHC molecules, such as 3-methyl pyridine, 1-ethylnapthalene, 1,1-diphenylethane, and benzofuran, were experimentally tested and compared with benchmark LOHCs (toluene and 9-ethylcarbazole) for hydrogenation using a series of commercial supported metal catalysts. The order of conversion into fully hydrogenated products at 200 °C was 3-methyl pyridine (100%) > 9-ethyl carbazole (86.4%) > 2,3-benzofuran (74%) > 1,1-diphenylethane (66.9%) > 1-ethylnapthalene (66.7%) > toluene (57%), further validating the AI-guided molecular design. This study demonstrates promise of LLM-driven molecular design in conjunction with ML-based screening for accelerated discovery and design of molecules.

Harb, Hassan [Argonne National Laboratory (ANL), A↗

Comparison of light out-coupling enhancements in single-layer blue-phosphorescent organic light emitting diodes using small-molecule or polymer hosts

Single-layer blue phosphorescence organic light emitting diodes (OLEDs) with either small-molecule or polymer hosts are fabricated using solution process and the performances of devices with different hosts are investigated. The small-molecule device exhibits luminous efficiency of 14.7 cd/A and maximum power efficiency of 8.39 lm/W, which is the highest among blue phosphorescence OLEDs with single-layer solution process and small molecular hosts. Using the same solution process for all devices, comparison of light out-coupling enhancement, with brightness enhancement film (BEF), between small-molecule and polymer based OLEDs is realized. Due to different dipole orientation and anisotropic refractive index, polymer-based OLEDs would trap less light than small molecule-based OLEDs internally, about 37% better based simulation results. In spite of better electrical and spectroscopic characteristics, including ambipolar characteristics, higher carrier mobility, higher photoluminescence quantum yield, and larger triplet state energy, the overall light out-coupling efficiency of small molecule-based devices is worse than that of polymer-based devices without BEF. However, with BEF for light out-coupling enhancement, the improved ratio in luminous flux and luminous efficiency for small molecule based device is 1.64 and 1.57, respectively, which are significantly better than those of PVK (poly-9-vinylcarbazole) devices. In addition to the theoretical optical simulation, the experimental data also confirm the origins of differential light-outcoupling enhancement. The maximum luminous efficiency and power efficiency are enhanced from 14.7 cd/A and 8.39 lm/W to 23 cd/A and 13.2 lm/W, respectively, with laminated BEF, which are both the highest so far for single-layer solution-process blue phosphorescence OLEDs with small molecule hosts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Rotation mechanism of methylammonium molecules in organometal halide perovskite in cubic phase: An ab initio molecular dynamics study

Rotation of methylammonium (CH 3 NH 3 or MA) molecules is believed to govern the excellent transport properties of photocarriers in the MA lead iodide (MAPbI 3 ) perovskite. Of particular interest is its cubic phase, which exists in industrially important films at room temperature. In order to investigate the rotational behaviors of the MA molecules, we have performed ab initio molecular dynamics simulations of cubic-MAPbI 3 at room temperature. There are two types of rotational motions of MA molecules in a crystalline PbI 3 cage: reorientation of a whole molecule and intramolecular rotation around the C–N bond within MA molecules. Using a cubic symmetry-assisted analysis (CSAA), we found that the prominent orientation of the C–N bond is the crystalline ❬110❭directions, rather than the ❬100❭ and ❬111❭ directions. Rapid rotation around the C–N bond is also observed, which easily occurs when the rotational axis is parallel to the ❬110❭ directions according to the CSAA. To explain the atomistic mechanisms underlying these CSAA results, we have focused on the relation between H–I hydrogen bonds and the orientation of an MA molecule. Here, the hydrogen bonds were defined by population analysis, and it has been found that, while H atoms in the CH 3 group (HC) hardly interacts with I atoms, those in the NH3 group (HN) form at least one hydrogen bond with I atoms and their interatomic distances are in a wide range, 2.2–3.7 Å. Based on these findings, we have given a possible explanation to why the ❬110❭ directions are preferred. Namely, the atomic arrangement and interatomic distance between MA and surrounding I atoms are most suitable for the formation of hydrogen bonds. In addition to films, these results are potentially applicable to the rotational behaviors in bulk MAPbI 3 as well, considering that the atomistic structure and time constants regarding the rotation of MA molecules statistically agree with bulk experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The ratio of the number of states in asymmetric and symmetric ozone molecules deviates from the statistical value of 2

Accurate calculations of vibrational states in singly and doubly substituted ozone molecules are carried out, up to dissociation threshold. Analysis of these spectra reveals noticeable deviations from the statistical factor of 2 for the ratio between the number of states in asymmetric and symmetric ozone molecules. It is found that, for the lower energy parts of spectra, the ratio is less than 2 in the singly substituted ozone molecules, but it is more than 2 in the doubly substituted ozone molecules. However, the upper parts of spectra, just below dissociation thresholds, exhibit a different behavior. In this energy range, the singly and doubly substituted ozone molecules behave similar, with the ratio of states in asymmetric and symmetric ozone molecules being more than 2 in both cases. Here, this property may contribute to an explanation of the mysterious η-effect in the ozone forming reaction that favors the formation of the asymmetric ozone molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An evaluation for geometries, formation enthalpies, and dissociation energies of diatomic and triatomic (C, H, N, O), NO 3 , and HNO 3 molecules from the PAW DFT method with PBE and optB88-vdW functionals

The structural geometries, formation enthalpies, and dissociation energies of all diatomic and triatomic molecules consisting of the four basic elements C, H, N, and/or O are calculated using the projector augmented wave density functional theory (DFT) method with the Perdew–Burke–Ernzerhof and optB88-vdW exchange-correlation functionals. The calculations are also extended to two larger molecules NO 3 and HNO 3 , which consist of four and five atoms, respectively. In total, 82 molecules or isomers are considered in the calculations. The geometric parameters including 42 bond lengths and 15 bond angles of these molecules calculated using the planewave DFT method are highly satisfactory, relative to the available experimental data. The error analysis is also performed for 49 formation enthalpies and 138 dissociation energies (including 51 atomization energies as well as the corresponding bond dissociation energies). The results are also compared with the previous data from various atomic-orbital-based methods for molecules and from similar or different planewave DFT methods for various solids and other molecules. This provides an informative and instructive evaluation especially for calculating the large-size material systems containing these small molecules as well as for developing the DFT methods further.

74 ATOMIC AND MOLECULAR PHYSICS↗

Structure and dynamics of C60 molecules on Au(111)

Earlier studies of C60 adsorption on Au(111) reported many interesting and complex features. We have performed coordinated low-energy electron diffraction, scanning tunneling microscopy (STM), and density functional theory studies to elucidate some of the details of the monolayer commensurate (2√3 × 2√3)R30° phase. We have identified the adsorption geometries of the two states that image as dim and bright in STM. These consist of a C60 molecule with a hexagon side down in a vacancy (hex-vac) and a C60 molecule with a carbon-carbon 6:6 bond down on a top site (6:6-top), respectively. We have studied the detailed geometries of these states and find that there is little distortion of the C60 molecules, but there is a rearrangement of the substrate near the C60 molecules. The two types of molecules differ in height, by about 0.7 Å, which accounts for most of the difference in their contrast in the STM images. The monolayer displays dynamical behavior, in which the molecules flip from bright to dim, and vice versa. We interpret this flipping as the result of the diffusion of vacancies in the surface layers of the substrate. Our measurements of the dynamics of this flipping from one state to the other indicate that the activation energy is 0.66 ± 0.03 eV for flips that involve nearest-neighbor C60 molecules, and 0.93 ± 0.03 for more distant flips. Based on calculated activation energies for vacancies diffusing in Au, we interpret these to be a result of surface vacancy diffusion and bulk vacancy diffusion. These results are compared to the similar system of Ag(111)-(2√3 × 2√3)R30°-C60. In both systems, the formation of the commensurate C60 monolayer produces a large number of vacancies in the top substrate layer that are highly mobile, effectively melting the interfacial metal layer at temperatures well below their normal melting temperatures.

36 MATERIALS SCIENCE↗