SEARCH · Engineering Papers
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.
Molecules with HOMO-LUMO gap calculated using DFTB. The molecules are a subset of a dataset generated for previous work (https://www.biorxiv.org/content/10.1101/2021.12.10.471928v1), which augmented the Enamine REAL database (https://enamine.net/). Molecules may be available for purchase through Enamine.
In view of the tremendous progress in atomically precise metal nanoclusters where electrochemical and optical energetics are routinely supported by computations to establish structure−function correlations, we explore the relationship between different protocols for measuring and computing band gaps of two distinct organic ligand-protected nanoclusters: [Cu 14 H 10 (MBN) 3 (PPh 3 ) 8 ] + and Au 20 (TBBT) 16 . Through UV/visible spectroscopy and differential pulse voltammetry, we measure optical and electrochemical band gaps in those systems. We then compare these experimentally determined gaps to HOMO−LUMO gaps, fundamental gaps, vertical excitation energies, and E o ox − E o red potentials computed using different density functional theory (DFT) or time-dependent DFT (TDDFT) methods. Specifically, in both copper and gold nanoclusters, we test the effect of truncating inert ligands from the model and compare density functionals with varying degrees of Hartree−Fock (HF) exchange from 0 to 50%, range-separated hybrids with a varying long-range tuning parameter, different correlation functionals, basis sets, and (equilibrium and nonequilibrium) continuum solvation models. Despite having different frontier orbital characters (the copper nanocluster has a metal-to-ligand charge transfer character while the gold nanocluster has metal-centered frontier orbitals), both nanoclusters display a similar sensitivity of the HOMO−LUMO gap to the HF exchange that is partially mitigated when computing the fundamental, optical, and electrochemical gaps. Other factors, such as the nature of the correlation functional, basis set, and geometry relaxation, have a considerably smaller effect on computed band gaps in these systems. Overall, this work provides guidelines for factors of varied importance for correlating computed and experimental band gap values.
Topological phases in laterally confined low-dimensional nanographenes have emerged as versatile design tools that can imbue otherwise unremarkable materials with exotic band structures ranging from topological semiconductors and quantum dots to intrinsically metallic bands. The periodic boundary conditions that define the topology of a given lattice have thus far prevented the translation of this technology to the quasi-zero-dimensional (0D) domain of small molecular structures. Here, we describe the synthesis of a polycyclic aromatic hydrocarbon (PAH) featuring two localized zero modes (ZMs) formed by the topological junction interface between a trivial and nontrivial phase within a single molecule. First-principles density functional theory calculations predict a strong hybridization between adjacent ZMs that gives rise to an exceptionally small HOMO–LUMO gap. Scanning tunneling microscopy and spectroscopy corroborate the molecular structure of 9/7/9-double quantum dots and reveal an experimental quasiparticle gap of 0.16 eV, corresponding to a carbon-based small molecule long-wavelength infrared (LWIR) absorber.
Abstract Graph Convolutional Neural Network (GCNN) is a popular class of deep learning (DL) models in material science to predict material properties from the graph representation of molecular structures. Training an accurate and comprehensive GCNN surrogate for molecular design requires large-scale graph datasets and is usually a time-consuming process. Recent advances in GPUs and distributed computing open a path to reduce the computational cost for GCNN training effectively. However, efficient utilization of high performance computing (HPC) resources for training requires simultaneously optimizing large-scale data management and scalable stochastic batched optimization techniques. In this work, we focus on building GCNN models on HPC systems to predict material properties of millions of molecules. We use HydraGNN, our in-house library for large-scale GCNN training, leveraging distributed data parallelism in PyTorch. We use ADIOS, a high-performance data management framework for efficient storage and reading of large molecular graph data. We perform parallel training on two open-source large-scale graph datasets to build a GCNN predictor for an important quantum property known as the HOMO-LUMO gap. We measure the scalability, accuracy, and convergence of our approach on two DOE supercomputers: the Summit supercomputer at the Oak Ridge Leadership Computing Facility (OLCF) and the Perlmutter system at the National Energy Research Scientific Computing Center (NERSC). We present our experimental results with HydraGNN showing (i) reduction of data loading time up to 4.2 times compared with a conventional method and (ii) linear scaling performance for training up to 1024 GPUs on both Summit and Perlmutter.
Tetrapyrrole macrocycles serve a multitude of roles in biological systems, including oxygen transport by heme and light harvesting and charge separation by chlorophylls and bacteriochlorophylls. Synthetic tetrapyrroles are utilized in diverse applications ranging from solar-energy conversion to photomedicine. Nevertheless, students beginning tetrapyrrole research, as well as established practitioners, are often puzzled when comparing properties of related tetrapyrroles. Questions arise as to why optical spectra of two tetrapyrroles often shift in wavelength/energy in a direction opposite to that predicted by common chemical intuition based on the size of a π-electron system. Gouterman's four-orbital model provides a framework for understanding these optical properties. Similarly, it can be puzzling as to why the oxidation potentials differ significantly when comparing two related tetrapyrroles, yet the reduction potentials change very little or shift in the opposite direction. In order to understand these redox properties, it must be recognized that structural/electronic alterations affect the four frontier molecular orbitals (HOMO, LUMO, HOMO-1 and LUMO+1) unequally and in many cases the LUMO+1, and not the LUMO, may track the HOMO in energy. This perspective presents a fundamental framework concerning tetrapyrrole electronic properties that should provide a foundation for rational molecular design in tetrapyrrole science.
We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.
A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.
Complex organic–inorganic interfaces are important for device and sensing applications. Charge transfer doping is prevalent in such applications and can affect the interfacial energy level alignments (ELA), which are determined by many-body interactions. We develop an approximate ab initio many-body GW approach that can capture many-body interactions due to interfacial charge transfer. The approach uses significantly less resources than a regular GW calculation but gives excellent agreement with benchmark GW calculations on an F4TCNQ/graphene interface. We find that many-body interactions due to charge transfer screening result in gate-tunable F4TCNQ HOMO–LUMO gaps. We further predict the ELA of a large system of experimental interest—4,4'-bis(dimethylamino)bipyridine (DMAP-OED) on monolayer MoS 2 , where charge transfer screening results in an ~1 eV reduction of the molecular HOMO–LUMO gap. Comparison with a two-dimensional electron gas model reveals the importance of explicitly considering the intraband transitions in determining the charge transfer screening in organic–inorganic interface systems.
Organic π-conjugated polymers with a triplet ground state have been the focus of recent research for their interesting and unique electronic properties, arising from the presence of the two unpaired electrons. These compounds are usually built from alternating electron-donating and electron-accepting monomer pairs which lower the HOMO–LUMO gap and yield a triplet state instead of the typical singlet ground state. In this paper, we use density functional theory calculations to explore the design rules that govern the creation of a ground-state triplet conjugated polymer and find that a small HOMO–LUMO gap in the singlet state is the best predictor for the existence of a triplet ground state, compared to previous use of a pro-quinoidal bonding character. Furthermore, this work can accelerate the discovery of new stable triplet materials by reducing the computational resources needed for electronic-state calculations and the number of potential candidates for synthesis.
We revisit the problem of unphysical charge density delocalization/fractionalization induced by the self-interaction error of common approximate Kohn–Sham (KS) density functional theory functionals on simulation of small to medium-sized proteins in a vacuum. Aside from producing unphysical electron densities and total energies, the vanishing of the HOMO–LUMO gap associated with the unphysical charge delocalization leads to an unphysical low-energy spectrum and catastrophic failure of most popular solvers for the KS self-consistent field (SCF) problem. We apply a robust quasi-Newton SCF solver to obtain solutions for some of these difficult cases. The anatomy of the charge delocalization is revealed by the natural deformation orbitals obtained from the density matrix difference between the Hartree–Fock and KS solutions; the charge delocalization not only can occur between charged fragments (such as in zwitterionic polypeptides) but also involves neutral fragments. The vanishing-gap phenomenon and troublesome SCF convergence are both attributed to the unphysical KS Fock operator eigenspectra of molecular fragments (e.g., amino acids or their side chains). Analysis of amino acid pairs suggests that the unphysical charge delocalization can be partially ameliorated by the use of some range-separated hybrid functionals but not by semilocal or standard hybrid functionals. Last, we demonstrate that solutions without the unphysical charge delocalization can be located even for semilocal KS functionals highly prone to such defects, but such solutions have non-Aufbau character and are unstable with respect to mixing of the non-overlapping “frontier” orbitals. Caution should be exercised when unexpectedly small (or vanishing) HOMO–LUMO gaps and atypical SCF convergence patterns (e.g., oscillatory) are observed in KS DFT simulations in any context (bio or otherwise).
Among C 60 ’s diverse functionalities, its potential application in CO 2 sequestration has gained increasing interest. However, the processes involved are sensitive to the molecule’s electronic structure, aspects of which remain debated and require greater precision. To address this, we performed structural optimization of fullerene C 60 using the QM MP2/6–31G* method. The nonplanarity of the optimized icosahedron is characterized by two types of dihedral angles: 138° and 143°. The 120 dihedrals of 138° occur between two hexagons intersecting at C–C bonds of 1.42 Å, while the 60 dihedrals of 143° are observed between hexagons and pentagons at C–C bonds of 1.47 Å. NBO analysis reveals less pyramidal sp 1.78 hybridization for carbons at the 1.42 Å bonds and more pyramidal sp 2.13 hybridization for the 1.47 Å bonds. Electrostatic potential charges range from -0.04 a.u. to 0.04 a.u. on the carbon atoms. Second-order perturbation analysis indicates that delocalization interactions in the C–C bonds of 1.42 Å (143.70 kcal/mol) and 1.47 Å (34.98 kcal/mol) are 22% and 38% higher, respectively, than those in benzene. MP2/Def2SVP calculations yield a correlation energy of 13.49 kcal/mol per electron for C 60 , slightly higher than the 11.68 kcal/mol for benzene. However, the results from HOMO-LUMO calculations should be interpreted with caution. This study may assist in the rational design of fullerene C 60 derivatives for CO 2 reduction systems.
Density functional theory (DFT) and beyond-DFT methods are often used in combination with photoelectron spectroscopy to obtain physical insights into the electronic structure of molecules and solids. The Kohn–Sham eigenvalues are not electron removal energies except for the highest occupied orbital. The eigenvalues of the highest occupied molecular orbitals often underestimate the electron removal or ionization energies due to the self-interaction (SI) errors in approximate density functionals. Here, we adapt and implement the density-consistent effective potential method of Kohut, Ryabinkin, and Staroverov (2014) to obtain SI-corrected local effective potentials from the SI-corrected Fermi–Löwdin orbitals and density in the Fermi–Löwdin orbital self-interaction correction scheme. The implementation is used to obtain the density of states (photoelectron spectra) and HOMO–LUMO gaps for a set of molecules and polyacenes. Good agreement with experimental values is obtained compared to a range of SI uncorrected density functional approximations.
First-principles calculations for electrochemistry require accurate treatment of both electronic structure and solvation. The perturbative GW approximation starting from density functional theory (DFT) calculations accurately models materials systems with varying dimensionality. Continuum solvation models enable efficient treatment of solvation effects in DFT calculations, but their applications with beyond-DFT electronic structure methods such as GW have been limited. Here, we introduce the frequency-dependent liquid polarizability from a nonlocal continuum solvation model in the screened Coulomb interaction of full-frequency GW calculations with a solvated DFT starting point. In this study, we show that the liquid screening contributions substantially reduce the HOMO–LUMO gap of molecules by 3–5 eV, while solvent effects on the DFT starting point negligibly impact the GW gap. The resulting framework facilitates the simultaneous electronic and solvation accuracy needed for first-principles electrochemistry.
Evaluation of the electron-nuclear dynamics and relaxation mechanisms of gold and silver nanoclusters and their alloys is important for future photocatalytic, light harvesting, and photoluminescence applications of these systems. In this work, the effect of silver doping on the nonradiative excited state relaxation dynamics of the atomically precise thiolate-protected gold nanocluster [Au 25- n Ag n (SH) 18 ] –1 ( n = 1, 12, 25) is studied theoretically. Time-dependent density functional theory is used to study excited states lying in the energy range 0.0–2.5 eV. The fewest switches surface hopping method with decoherence correction was used to investigate the dynamics of these states. The HOMO–LUMO gap increases significantly upon doping of 12 silver atoms but decreases for the pure silver nanocluster. Doped clusters show a different response for ground state population increase lifetimes and excited state population decay times in comparison to the undoped system. The ground state recovery times of the S 1 –S 6 states in the first excited peak were found to be longer for [Au 13 Ag 12 (SH) 18 ] –1 than the corresponding recovery times of other studied nanoclusters, suggesting that this partially doped nanocluster is best for preserving electrons in an excited state. The decay time constants were in the range of 2.0–20 ps for the six lowest energy excited states. Among the higher excited states, S 7 has the slowest decay time constant although it occurs more quickly than S 1 decay. Overall, these clusters follow common decay time constant trends and relaxation mechanisms due to the similarities in their electronic structures.
Lignin valorization is restricted by the stability of its methoxy groups, creating a critical need for efficient demethylation strategies. Here, in this study, density functional theory (DFT) is employed to dissect the mechanistic pathways of demethylation in lignin model compounds, guaiacol and syringol, using protic ionic liquids (PILs) that act as both solvent and catalyst. Conductor‐like Screening Model for Real Solvents (COSMO‐RS) analysis identifies monoethanolammonium acetate ([MEOA][Ace]) as the most promising medium, attributed to its strong hydrogen bonding network and solvation ability. By integrating implicit and explicit solvation models, it is revealed that an acid‐catalyzed hydrolytic mechanism governs demethylation, with PILs stabilizing crucial transition states and intermediates. Complementary electronic structure evaluations, including highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO‐LUMO) gap analysis, charge distribution, and electrostatic potential mapping, demonstrate how PILs lower energetic barriers and enhance reactivity. To simulate realistic environments, this study is extended to lignin dimer complexes with varying water content, uncovering how water‐bridged solvation reverses demethylation preference from guaiacyl to syringyl units. This mechanistic shift aligns with experimental observations showing faster S‐unit reactivity in hydrated systems. Together, these findings provide atomic‐level insight into lignin demethylation dynamics and highlight how tuning acid concentration in PILs can accelerate kinetics, enabling a rational pathway toward next‐generation biomass conversion technologies.
Ultrasmall metal nanoclusters (NCs) exhibit a quantized conduction band, hence, a distinct HOMO-LUMO gap (E g ). Such a quantized electronic structure gives rise to multiple discrete peaks in the optical absorption spectrum of the NCs. As the size grows to 130 gold atoms (Au 130 protected by ligands), electrical charging and optical behaviors seem to show certain metal-like features (hence, transition-sizes). To probe such behaviors, especially the potential ligand effect, we have devised the synthesis of Au 130 NCs protected by phenylethanethiolate and naphthalenethiolate, respectively, with the former having a nonconjugated separation between the aromatic molecular group and the metal core while the latter being in direct bonding. A careful comparison of these two Au 130 nanoclusters with the earlier reported analogues is carried out, including the structurally characterized Au 130 (pMBT) 50 and the aqueous counterpart. While all of these nanoclusters possess the same 80 free electron counts in the core, some notable differences in electrochemical and optical properties are found, which are attributed to the ligand effects. The obtained insights may stimulate further interest in the transition-sized nanoclusters and also promote their applications in optics, energy conversion, and biomedicine.
Developing an atomistic understanding of ionizing radiation induced changes to organic materials is necessary for intentional design of greener and more sustainable materials for radiation shielding and detection. Cocrystals are promising for these purposes, but a detailed understanding of how the specific intermolecular interactions within the lattice upon exposure to radiation affect the structural stability of the organic crystalline material is unknown. This study evaluates atomistic-level effects of γ radiation on both single- and multicomponent organic crystalline materials and how specific noncovalent interactions and packing within the crystalline lattice enhance structural stability. Dose studies were performed on all crystalline systems and evaluated via experimental and computational methods. Changes in crystallinity were evaluated by p-XRD and free radical formation was analyzed via EPR spectroscopy. Type of intermolecular interactions and packing within the crystal lattice was delineated and related to the specific free radical species formed and the structural integrity of each material. Periodic DFT and HOMO-LUMO surface mapping calculations provided atomistic-level identifications of the most probable sites for the radicals formed upon exposure to γ radiation and relate intermolecular interactions and molecular packing within the crystalline lattice to experimental results.
Electron transport in organic molecules and biomolecules is governed by electronic structure and molecular conformations. Despite recent progress, key challenges remain in understanding the role of intramolecular interactions and three-dimensional (3D) conformations on the electron transport behavior of organic molecules. In this work, the electronic properties of aromatic amide foldamers are characterized that organize into distinct 3D structures, including an extended secondary amide that adopts a trans-conformation and a folded N-methylated tertiary amide that adopts a cis-conformation. Results from single-molecule electronic experiments show that the extended secondary amide exhibits a fourfold enhancement in molecular conductance compared to the folded N-methylated tertiary amide, despite a longer contour length. The results show that extended amide molecules are governed by a through-bond electron transport mechanism, whereas folded amide molecules are dominated by through-space transport. Bulk spectroscopic characterization and density functional theory calculations further reveal that extended amides have a smaller HOMO–LUMO gap and larger transmission values compared to folded amides, consistent with single-molecule electronic experiments. Overall, this work shows that 3D molecular conformations significantly influence the electronic properties of single-molecule junctions.