An official website of the United States government
ABSTRACT In this work, we present simulations of the optical rotation (OR) for five molecular crystals at density functional theory level with periodic boundary conditions (DFT‐PBC). Calculations are compared with experimental measurements and show semi‐quantitative agreement with experimental data for three of the crystals: tartatic acid, benzil, and pentaerythritol. For the other two crystals, aspartic acid and glutamic acid, the calculated data are in qualitative agreement with, but two orders of magnitude smaller than, the experimental data. We provide some arguments that support the theoretical predictions and suggest that the experiments should be revisited. We also find that the position of H centers provided in experimental X‐ray data is not sufficiently reliable for simulating OR, and better results are obtained when H atoms are allowed to relax while keeping heavier elements fixed at the experimental positions. Comparison with molecular cluster calculations with a better functional and a larger basis set indicate that the role of intermolecular interactions (reproduced with the PBC technique) is as or more important than the choice of model chemistry. Despite the current limitations in the level of theory that can be employed, these simulations provide a promising avenue to investigate the effect of intermolecular interactions on this sensitive electronic property of molecules and materials.
more »
« less
Surface Engineering of the Mechanical Properties of Molecular Crystals via an Atomistic Model for Computing the Facet Stress Response of Solids
Abstract Dynamic molecular crystals are an emerging class of crystalline materials that can respond to mechanical stress by dissipating internal strain in a number of ways. Given the serendipitous nature of the discovery of such crystals, progress in the field requires advances in computational methods for the accurate and high–throughput computation of the nanomechanical properties of crystals on specific facets which are exposed to mechanical stress. Here, we develop and apply a new atomistic model for computing the surface elastic moduli of crystals on any set of facets of interest using dispersion–corrected density functional theory (DFT−D) methods. The model was benchmarked against a total of 24 reported nanoindentation measurements from a diverse set of molecular crystals and was found to be generally reliable. Using only the experimental crystal structure of the dietary supplement, L–aspartic acid, the model was subsequently applied under blind test conditions, to correctly predict the growth morphology, facet and nanomechanical properties of L–aspartic acid to within the accuracy of the measured elastic stiffness of the crystal, 24.53±0.56 GPa. This work paves the way for the computational design and experimental realization of other functional molecular crystals with tailor–made mechanical properties.
more »
« less
- Award ID(s):
- 2137264
- PAR ID:
- 10642380
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Chemistry – A European Journal
- Volume:
- 30
- Issue:
- 38
- ISSN:
- 0947-6539
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Heng, Jerry (Ed.)The morphological evolution of organic crystals during crystallization depends on the face-specific growth rates. Classical growth rate models relate the face-specific growth rates to the crystal lattice, energy of stable facets, growth mechanism, and supersaturation. The complexities of these models have increased over time to account accurately for solution conditions, the structure of growth units, and their attachment rates. Such advanced growth rate models require several layers of computations to obtain attachment energies of facets, nucleation rates, kink density, and attachment rates. Among these, the most intensive and time-consuming computation is for attachment rates, which require molecular dynamic simulations. This substantially increases the overall computation time to predict the absolute growth rate for even one crystallization condition. Since it is nearly impossible to iterate such a growth rate model, optimization schemes cannot be implemented to identify solution conditions that favor specific crystal growth. To reduce the computational time for attachment rate calculations, we implement a group contribution method (GCM) that relates the properties of functional groups in a molecule to their attachment rates to the crystal lattice, thereby rapidly estimating the growth rates of organic crystals. The process of molecular attachment involves partial desolvation of a solvated molecule, referred to as a transition state, followed by total desolvation via spontaneous attachment to a crystal facet. The first step in GCM is to identify the equilibrium states of fully solvated and partially desolvated solute molecules. The degree of supersaturation dictates the extent of this equilibrium and, thereby, the activation barrier for the growth of crystals, according to transition state theory. Identifying this equilibrium phenomenon allows for capturing the functional-group-specific interactions that depend on molecular motion, which could be related to operating conditions such as temperature and pressure. The stochastic optimization technique with Monte-Carlo sampling allows an efficient optimization problem solution to obtain the group interaction parameters. The GCM approach is first validated for the estimation of growth rates of glutamic acid and L-histidine, and then extended to predict growth rates of alanine and glycine rapidly. The optimized parameters and GCM scheme can be used to estimate growth rates in other crystallization systems.more » « less
-
null (Ed.)Abstract Transition metal dichalcogenides (TMDs) offer superior properties over conventional materials in many areas such as in electronic devices. In recent years, TMDs have been shown to display a phase switching mechanism under the application of external mechanical strain, making them exciting candidates for phase change transistors. Molybdenum ditelluride (MoTe2) is one such material that has been engineered as a strain-based phase change transistor. In this work, we explore various aspects of the mechanical properties of this material by a suite of computational and experimental approaches. First, we present parameterization of an interatomic potential for modeling monolayer as well as multilayered MoTe2 films. For generating the empirical potential parameter set, we fit results from density functional theory calculations using a random search algorithm known as particle swarm optimization. The potential closely predicts structural properties, elastic constants, and vibrational frequencies of MoTe2 indicating a reliable fit. Our simulated mechanical response matches earlier larger scale experimental nanoindentation results with excellent prediction of fracture points. Simulation of uniaxial tensile deformation by molecular dynamics shows the complete non-linear stress-strain response up to failure. Mechanical behavior, including failure properties, exhibits directional anisotropy due to the variation of bond alignments with crystal orientation. Furthermore, we show the deterioration of mechanical properties with increasing temperature. Finally, we present computational and experimental evidence of an extended c-axis strain transfer length in MoTe2 compared to TMDs with smaller chalcogen atoms.more » « less
-
null (Ed.)9-Anthracene carboxylic acid derivatives comprise a family of thermally reversible photomechanical molecular crystals. The photomechanical response relies on a [4 + 4] photodimerization followed by dissociation that occurs on timescales of seconds to minutes. A combined theoretical and experimental investigation is undertaken to better understand how chemical modification of the anthracene core influences energetics of both the isolated molecule and the crystal lattice. We use both density functional theory and dispersion-corrected Moller–Plesset perturbation theory computational methods to establish orbital energies, photodimerization reaction energies, and lattice energies for a set of substituted 9-anthracene carboxylic acid molecules. The calculations reveal that steric interactions play a dominant role in the ability to form photodimers and indicate an energetic threshold of 80–90 kJ per mole for the dimerization reaction. Examination of intermolecular bonding in a subset of fluorinated 9ACs revealed the absence of H⋯F intermolecular bond formation and energy differences that can explain observed trends in the dissociation kinetics and mechanical reset times. Fluorescence recovery after photobleaching experiments shows that the photodimer dissociation kinetics depend on the amount of initial photodimer, preventing a straightforward correlation between halogen atom substitution and dissociation rates using the Bell–Evans–Polanyi principle. The results clarify how molecular structure affects intermolecular interactions and photoreactivity in this family of molecular crystals, but the origin of the complex photodimer dissociation dynamics remains an open question.more » « less
-
Abstract Like silicon, single crystals of organic semiconductors are pursued to attain intrinsic charge transport properties. However, they are intolerant to mechanical deformation, impeding their application in flexible electronic devices. Such contradictory properties, namely exceptional molecular ordering and mechanical flexibility, are unified in this work. We found that bis(triisopropylsilylethynyl)pentacene (TIPS‐P) crystals can undergo mechanically induced structural transitions to exhibit superelasticity and ferroelasticity. These properties arise from cooperative and correlated molecular displacements and rotations in response to mechanical stress. By utilizing a bending‐induced ferroelastic transition of TIPS‐P, flexible single‐crystal electronic devices were obtained that can tolerate strains (ϵ) of more than 13 % while maintaining the charge carrier mobility of unstrained crystals (μ>0.7 μ0). Our work will pave the way for high‐performance ultraflexible single‐crystal organic electronics for sensors, memories, and robotic applications.more » « less
