Ab Initio Fragment Method for Calculating Molecular X-ray Diffraction

Extracting the electronic structure signal from X-ray and electron scattering in the gas phase

X-ray and electron scattering from free gas-phase molecules is examined using the independent atom model (IAM) and ab initio electronic structure calculations. The IAM describes the effect of the molecular geometry on the scattering, but does not account for the redistribution of valence electrons due to, for instance, chemical bonding. By examining the total, i.e. energy-integrated, scattering from three molecules, fluoroform (CHF3), 1,3-cyclohexadiene (C6H8) and naphthalene (C10H8), the effect of electron redistribution is found to predominantly reside at small-to-medium values of the momentum transfer (q ≤ 8 Å-1) in the scattering signal, with a maximum percent difference contribution at 2 ≤ q ≤ 3 Å-1. A procedure to determine the molecular geometry from the large-q scattering is demonstrated, making it possible to more clearly identify the deviation of the scattering from the IAM approximation at small and intermediate q and to provide a measure of the effect of valence electronic structure on the scattering signal.

Automatic Clustering of Excited-State Trajectories: Application to Photoexcited Dynamics

We introduce automatic clustering as a computationally efficient tool for classifying and interpreting trajectories from simulations of photo-excited dynamics. Trajectories are treated as time-series data, with the features for clustering selected by variance mapping of normalized data. The L2-norm and dynamic time warping are proposed as suitable similarity measures for calculating the distance matrices, and these are clustered using the unsupervised density-based DBSCAN algorithm. The silhouette coefficient and the number of trajectories classified as noise are used as quality measures for the clustering. The ability of clustering to provide rapid overview of large and complex trajectory data sets, and its utility for extracting chemical and physical insight, is demonstrated on trajectories corresponding to the photochemical ring-opening reaction of 1,3-cyclohexadiene, noting that the clustering can be used to generate reduced dimensionality representations in an unbiased manner.

Mapping the photochemistry of cyclopentadiene: from theory to ultrafast X-ray scattering

The photoinduced ring-conversion reaction when cyclopentadiene (CP) is excited at 5.10 eV is simulated using surface-hopping semiclassical trajectories with XMS(3)-CASPT2(4,4)/cc-pVDZ electronic structure theory. In addition, PBE0/def2-SV(P) is employed for ground state propagation of the trajectories. The dynamics is propagated for 10 ps, mapping both the nonadiabatic short-time dynamics (<300 fs) and the increasingly statistical dynamics on the electronic ground state. The short-time dynamics yields a mixture of hot CP and bicyclo[2.1.0]pentene (BP), with the two products reached via different regions of the same conical intersection seam. On the ground state, we observe slow conversion from BP to CP which is modelled by RRKM theory with a transition state determined using PBE0/def2-TZVP. The CP products are furthermore associated with ground state hydrogen shifts and some H-atom dissociation. Finally, the prospects for detailed experimental mapping using novel ultrafast X-ray scattering experiments are discussed and observables for such experiments are predicted. In particular, we assess the possibility of retrieving electronic states and their populations alongside the structural dynamics.

Fragmentation and transferability in Hirshfeld atom refinement

Hirshfeld atom refinement (HAR) is one of the most effective methods for obtaining accurate structural parameters for hydrogen atoms from X-ray diffraction data. Unfortunately, it is also relatively computationally expensive, especially for larger molecules due to wavefunction calculations. Here, a fragmentation approach has been tested as a remedy for this problem. It gives an order of magnitude improvement in computation time for larger organic systems and is a few times faster for metal-organic systems at the cost of only minor differences in the calculated structural parameters when compared with the original HAR calculations. Fragmentation was also applied to polymeric and disordered systems where it provides a natural solution to problems that arise when HAR is applied. The concept of fragmentation is closely related to the transferable aspherical atom model (TAAM) and allows insight into possible ways to improve TAAM. Hybrid approaches combining fragmentation with the transfer of atomic densities between chemically similar atoms have been tested. An efficient handling of intermolecular interactions was also introduced for calculations involving fragmentation. When applied in fragHAR (a fragmentation approach for polypeptides) as a replacement for the original approach, it allowed for more efficient calculations. All of the calculations were performed with a locally modified version of Olex2 combined with a development version of discamb2tsc and ORCA. Care was taken to efficiently use the power of multicore processors by simple implementation of load-balancing, which was found to be very important for lowering computational time.

Towards high-resolution X-ray scattering as a probe of electron correlation

We demonstrate that X-ray scattering can be used as a probe of electron–electron correlation.

Imaging of femtosecond bond breaking and charge dynamics in ultracharged peptides

X-ray free-electrons lasers have revolutionized the method of imaging biological macromolecules such as proteins, viruses and cells by opening the door to structural determination of both single particles and crystals at room temperature. By utilizing high intensity X-ray pulses on femtosecond timescales, the effects of radiation damage can be reduced. Achieving high resolution structures will likely require knowledge of how radiation damage affects the structure on an atomic scale, since the experimentally obtained electron densities will be reconstructed in the presence of radiation damage. Detailed understanding of the expected damage scenarios provides further information, in addition to guiding possible corrections that may need to be made to obtain a damage free reconstruction. In this work, we have quantified the effects of ionizing photon-matter interactions using first principles molecular dynamics. We utilize density functional theory to calculate bond breaking and charge dynamics in three ultracharged molecules and two different structural conformations that are important to the structural integrity of biological macromolecules, comparing to our previous studies on amino acids. The effects of the ultracharged states and subsequent bond breaking in real space are studied in reciprocal space using coherent diffractive imaging of an ensemble of aligned biomolecules in the gas phase.

Quantum-Mechanical Structure Optimization of Protein Crystals and Analysis of Interactions in Periodic Systems

A fast quantum-mechanical approach, density-functional tight-binding combined with the fragment molecular orbital method and periodic boundary conditions, is used to optimize atomic coordinates and cell parameters for a set of protein crystals: 1ETL, 5OQZ, 3Q8J, 1CBN, and 2VB1. Good agreement between experimental and calculated structures is obtained for both atomic coordinates and cell parameters. Sterical clashes present in the experimental structures are corrected by simulations. The partition analysis is extended to treat periodic boundary conditions and applied to analyze protein-solvent interactions in crystals.

fragHAR: towards ab initio quantum-crystallographic X-ray structure refinement for polypeptides and proteins

The first ab initio aspherical structure refinement against experimental X-ray structure factors for polypeptides and proteins using a fragmentation approach to break up the protein into residues and solvent, thereby speeding up quantum-crystallographic Hirshfeld atom refinement (HAR) calculations, is described. It it found that the geometric and atomic displacement parameters from the new fragHAR method are essentially unchanged from a HAR on the complete unfragmented system when tested on dipeptides, tripeptides and hexapeptides. The largest changes are for the parameters describing H atoms involved in hydrogen-bond interactions, but it is shown that these discrepancies can be removed by including the interacting fragments as a single larger fragment in the fragmentation scheme. Significant speed-ups are observed for the larger systems. Using this approach, it is possible to perform a highly parallelized HAR in reasonable times for large systems. The method has been implemented in the TONTO software.

On the use of the Obara–Saika recurrence relations for the calculation of structure factors in quantum crystallography

Modern methods of quantum crystallography are techniques firmly rooted in quantum chemistry and, as in many quantum chemical strategies, electron densities are expressed as two-centre expansions that involve basis functions centred on atomic nuclei. Therefore, the computation of the necessary structure factors requires the evaluation of Fourier transform integrals of basis function products. Since these functions are usually Cartesian Gaussians, in this communication it is shown that the Fourier integrals can be efficiently calculated by exploiting an extension of the Obara–Saika recurrence formulas, which are successfully used by quantum chemists in the computation of molecular integrals. Implementation and future perspectives of the technique are also discussed.