Dephasing and Dissipation in a Source–Drain Model of Light‐Harvesting Systems

The energy transport process in natural‐light‐harvesting systems is investigated by solving the time‐dependent Schrödinger equation for a source–network–drain model incorporating the effects of dephasing and dissipation, owing to coupling with the environment. In this model, the network consists of electronically coupled chromophores, which can host energy excitations (excitons) and are connected to source channels, from which the excitons are generated, thereby simulating exciton creation from sunlight. After passing through the network, excitons are captured by the reaction centers and converted into chemical energy. In addition, excitons can reradiate in green plants as photoluminescent light or be destroyed by nonphotochemical quenching (NPQ). These annihilation processes are described in the model by outgoing channels, which allow the excitons to spread to infinity. Besides the photoluminescent reflection, the NPQ processes are the main outgoing channels accompanied by energy dissipation and dephasing. From the simulation of wave‐packet dynamics in a one‐dimensional chain, it is found that, without dephasing, the motion remains superdiffusive or ballistic, despite the strong energy dissipation. At an increased dephasing rate, the wave‐packet motion is found to switch from superdiffusive to diffusive in nature. When a steady energy flow is injected into a site of a linear chain, exciton dissipation along the chain, owing to photoluminescence and NPQ processes, is examined by using a model with coherent and incoherent outgoing channels. It is found that channel coherence leads to suppression of dissipation and multiexciton super‐radiance. With this method, the effects of NPQ and dephasing on energy transfer in the Fenna–Matthews–Olson complex are investigated. The NPQ process and the photochemical reflection are found to significantly reduce the energy‐transfer efficiency in the complex, whereas the dephasing process slightly enhances the efficiency. The calculated absorption spectrum reproduces the main features of the measured counterpart. As a comparison, the exciton dynamics are also studied in a linear chain of pigments and in a multiple‐ring system of light‐harvesting complexes II (LH2) from purple bacteria by using the Davydov D₁ ansatz. It is found that the exciton transport shows superdiffusion characteristics in both the chain and the LH2 rings.     

Fractional occupation numbers and spin density functional calculations of degenerate systems

We implemented ab initio self‐consistent field (SCF) fractional occupation numbers (FON) calculation with Dunlap's interpolation scheme for the twisted ethylene, which is a prototype molecule of a σ–π biradical system. The calculational results are compared with those of complete‐active‐space (CAS) SCF and spin‐unrestricted Kohn–Sham (UKS) calculations on potential surfaces, occupation numbers of natural orbitals, and correlation entropies. It was found that the UKS methods gave similar results to CASSCF, while the FON solutions appeared in only the nearly complete degenerate region. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 93: 317–323, 2003     

Methodological approach to study energetic and structural properties of nanostructured cadmium sulfide by using ab-initio molecular dynamics simulations

A single wurtzite phase of cadmium sulfide cluster is investigated by ab-initio molecular dynamics simulations at different temperatures, ranging from 100 K to 600 K. In this study we propose a possible procedure to characterize the CdS quantum dots system by means of molecular dynamics calculations using a standard Car-Parrinello scheme. In order to ensure the accuracy of the numerical approach, preliminary calculations to test pseudopotentials, cutoff and box size on both single atoms systems and Cd–Cd, S–S, Cd–S dimers have been performed. Calculated binding energies and bond lengths are obtained in good agreement with experimental data. Subsequently, an uncapped CdS cluster with size below 2 nm, 48 atoms of cadmium and 48 atoms of sulfur, in a wurtzite geometry was structurally optimized to minimize internal stresses. The CdS cluster has been carefully characterized structurally at several temperatures up to T = 600 K. At the temperature of 340 K atomic diffusion on the surface allows the onset of a new stable atomic configuration.     

Myoglobin-CO substate structures and dynamics: multidimensional vibrational echoes and molecular dynamics simulations

Spectrally resolved infrared stimulated vibrational echo data were obtained for sperm whale carbonmonoxymyoglobin (MbCO) at 300 K. The measured dephasing dynamics of the CO ligand are in agreement with dephasing dynamics calculated with molecular dynamics (MD) simulations for MbCO with the residue histidine-64 (His64) having its imidazole epsilon nitrogen protonated (N(epsilon)-H). The two conformational substate structures B(epsilon) and R(epsilon) observed in the MD simulations are assigned to the spectroscopic A(1) and A(3) conformational substates of MbCO, respectively, based on the agreement between the experimentally measured and calculated dephasing dynamics for these substates. In the A(1) substate, the N(epsilon)-H proton and N(delta) of His64 are approximately equidistant from the CO ligand, while in the A(3) substate, the N(epsilon)-H of His64 is oriented toward the CO, and the N(delta) is on the surface of the protein. The MD simulations show that dynamics of His64 represent the major source of vibrational dephasing of the CO ligand in the A(3) state on both femtosecond and picosecond time scales. Dephasing in the A(1) state is controlled by His64 on femtosecond time scales, and by the rest of the protein and the water solvent on longer time scales.     

Multidimensional umbrella sampling and replica‐exchange molecular dynamics simulations for structure prediction of transmembrane helix dimers

Structural information of a transmembrane (TM) helix dimer is useful in understanding molecular mechanisms of important biological phenomena such as signal transduction across the cell membrane. Here, we describe an umbrella sampling (US) scheme for predicting the structure of a TM helix dimer in implicit membrane using the interhelical crossing angle and the TM–TM relative rotation angles as the reaction coordinates. This scheme conducts an efficient conformational search on TM–TM contact interfaces, and its robustness is tested by predicting the structures of glycophorin A (GpA) and receptor tyrosine kinase EphA1 (EphA1) TM dimers. The nuclear magnetic resonance (NMR) structures of both proteins correspond to the global free‐energy minimum states in their free‐energy landscapes. In addition, using the landscape of GpA as a reference, we also examine the protocols of temperature replica‐exchange molecular dynamics (REMD) simulations for structure prediction of TM helix dimers in implicit membrane. A wide temperature range in REMD simulations, for example, 250–1000 K, is required to efficiently obtain a free‐energy landscape consistent with the US simulations. The interhelical crossing angle and the TM–TM relative rotation angles can be used as reaction coordinates in multidimensional US and be good measures for conformational sampling of REMD simulations. © 2013 Wiley Periodicals, Inc.     

Vibrational Dephasing in Ionic Liquids as a Signature of Hydrogen Bonding

Understanding both structure and dynamics is crucial for producing tailor‐made ionic liquids (ILs). We studied the vibrational and structural dynamics of medium versus weakly hydrogen‐bonded C?H groups of the imidazolium ring in ILs of the type [1‐alkyl‐3‐methylimidazolium][bis(trifluoromethanesulfonyl)imide] ([C_nmim][NTf₂]), with n=1, 2, and 8, by time‐resolved coherent anti‐Stokes Raman scattering (CARS) and quantum‐classical hybrid (QCH) simulations. From the time series of the CARS spectra, dephasing times were extracted by modeling the full nonlinear response. From the QCH calculations, pure dephasing times were obtained by analyzing the distribution of transition frequencies. Experiments and calculations reveal larger dephasing rates for the vibrational stretching modes of C(2)?H compared with the more weakly hydrogen‐bonded C(4,5)?H. This finding can be understood in terms of different H‐bonding motifs and the fast interconversion between them. Differences in population relaxation rates are attributed to Fermi resonance interactions.     

Characterization of tunneling systems in molecular versus polymer glasses by high-pressure photon echo spectroscopy

Intrinsic differences between tunneling two-level systems (TLSs) in molecular versus polymeric glasses are revealed by studying the effect of compression on TLS dynamics. Photon echo studies under variable low-temperature (1.1-2.3 K) and high-pressure (0-30 kbar) conditions have been performed to contrast the effect of compression on molecular [2-methyl-tetrahydrofuran (2MTHF)] versus polymer [Polymethylmethacrylate (PMMA)] glasses. The pressure-induced reduction in the magnitude of the optical dephasing rate of rhodamine 640 in a molecular glass (2MTHF) is found to be comparable to the volume decrease of the glass (e.g., approximately 20% at 30 kbar), indicating that TLSs in 2MTHF are associated with void space or low-density regions of the glass. In contrast, the relative pressure insensitivity observed for organic polymer glasses (PMMA) supports the idea that these TLSs are associated with side chain defects. The power-law exponent for the temperature-dependent dephasing in 2MTHF also decreased significantly at high pressure, suggesting a change in the form of the TLS density of states upon compression of the molecular glass.     

Microsecond simulations of DNA and ion transport in nanopores with novel ion–ion and ion–nucleotides effective potentials

We developed a novel scheme based on the grand‐canonical Monte Carlo/Brownian dynamics simulations and have extended it to studies of ion currents across three nanopores with the potential for single‐stranded DNA (ssDNA) sequencing: solid‐state nanopore Si₃N₄, α‐hemolysin, and E111N/M113Y/K147N mutant. To describe nucleotide‐specific ion dynamics compatible with ssDNA coarse‐grained model, we used the inverse Monte Carlo protocol, which maps the relevant ion–nucleotide distribution functions from all‐atom molecular dynamics (MD) simulations. Combined with the previously developed simulation platform for Brownian dynamics simulations of ion transport, it allows for microsecond‐ and millisecond‐long simulations of ssDNA dynamics in the nanopore with a conductance computation accuracy that equals or exceeds that of all‐atom MD simulations. In spite of the simplifications, the protocol produces results that agree with the results of previous studies on ion conductance across open channels and provide direct correlations with experimentally measured blockade currents and ion conductances that have been estimated from all‐atom MD simulations. © 2014 Wiley Periodicals, Inc.     

Atomistic models of three fluorinated polyimides in the amorphous state

Molecular models of three fluorinated polyimides based on the 4,4′‐(hexafluoroisopropylidene)diphthalic dianhydride (6FDA) have been studied using molecular dynamics (MD) simulations. The respective diamines were 4,4′‐hexafluoroisopropylidene dianiline (6FpDA), 3,3′‐hexafluoroisopropylidene dianiline (6FmDA), and 2,4,6‐trimethyl‐1,3phenylenediamine (DAM). Thirty independent samples were prepared using a hybrid pivot Monte Carlo‐MD generation technique and average densities were found to be in very good agreement with experiment. Model structures also agreed with available wide‐angle X‐ray scattering data. Cohesive energies, Hildebrand solubility parameters, fractional free volumes (FFV), void space distributions and intermolecular as well as intramolecular interactions were analyzed. The differences in bulk properties between both 6FDA‐6FpDA and 6FDA‐6FmDA isomers remain fairly small, although the configurations of the former are more extended. 6FDA‐DAM has a lower density, larger intermolecular distances, and higher free volume than the other two polyimides. Results are discussed with respect to their use as matrices for gas separation. © 2009 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 47: 1166–1180, 2009     

Temperature‐dependent changes in the hydrogen bonded hard segment network and microphase morphology in a model polyurethane: Experimental and simulation studies

Hydrogen bonding between hard segments has a critical effect on the morphology and properties of polyurethanes. Influence of temperature on hydrogen bonded urethane network and melting behavior of a model semicrystalline segmented polyurethane was investigated by experiments and simulations. Polyurethane was synthesized by the stoichiometric reaction between p‐phenylene diisocyanate and poly(tetramethylene oxide) (PTMO) with a molecular weight of 1000 g/mol. Simulations were carried out using dissipative particle dynamics (DPD) and molecular dynamics (MD) approaches. Experimental melting behavior obtained by various techniques was compared with simulations. DPD simulations showed a room temperature microphase morphology consisting of a three‐dimensional hydrogen‐bonded urethane hard segment network in a continuous and amorphous PTMO matrix. The first‐order melting transitions of crystalline urethane hard segments observed during the continuous isobaric heating in DPD and MD simulations (340–360 K) were in reasonably good agreement with those observed experimentally, such as AFM (320–340 K), WAXS (330–360 K), and FTIR (320–350 K) measurements. Quantitative verification of the melting of urethane hard segments was demonstrated by sharp discontinuities in energy versus temperature plots obtained by MD simulations due to substantial decrease in the number of hydrogen bonds above 340 K. © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2018 , 56, 182–192     

Dual‐level direct dynamics study on the reactions of SH (SD) with F2

The dual‐level direct dynamics study is carried out in the reactions of SH (SD) + F₂, and enthalpies of formation for the species FSH (FSD) are evaluated using isodesmic reactions. Furthermore, the dynamics calculations are performed using the variational transition‐state theory (VTST) with interpolated single‐point energies (ISPE) method at the approximate QCISD(T)/6‐311++G(3df, 2pd)//MPW1K/6‐311G(d, p) level. Over the temperature range of 200–2000 K, theoretical rate constants evaluated using canonical variational transition state theory (CVT) are in good agreement with the available experimental values. The kinetic isotope effects for the title reactions are not obvious. © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 710–716, 2005     

Pure optical dephasing dynamics in semiconducting single-walled carbon nanotubes

We report a detailed study of ultrafast exciton dephasing processes in semiconducting single-walled carbon nanotubes employing a sample highly enriched in a single tube species, the (6,5) tube. Systematic measurements of femtosecond pump-probe, two-pulse photon echo, and three-pulse photon echo peak shift over a broad range of excitation intensities and lattice temperature (from 4.4 to 292 K) enable us to quantify the timescales of pure optical dephasing (T(2)(*)), along with exciton-exciton and exciton-phonon scattering, environmental effects as well as spectral diffusion. While the exciton dephasing time (T(2)) increases from 205 fs at room temperature to 320 fs at 70 K, we found that further decrease of the lattice temperature leads to a shortening of the T(2) times. This complex temperature dependence was found to arise from an enhanced relaxation of exciton population at lattice temperatures below 80 K. By quantitatively accounting the contribution from the population relaxation, the corresponding pure optical dephasing times increase monotonically from 225 fs at room temperature to 508 fs at 4.4 K. We further found that below 180 K, the pure dephasing rate (1/T(2)(*)) scales linearly with temperature with a slope of 6.7 ± 0.6 μeV/K, which suggests dephasing arising from one-phonon scattering (i.e., acoustic phonons). In view of the large dynamic disorder of the surrounding environment, the origin of the long room temperature pure dephasing time is proposed to result from reduced strength of exciton-phonon coupling by motional narrowing over nuclear fluctuations. This consideration further suggests the occurrence of remarkable initial exciton delocalization and makes nanotubes ideal to study many-body effects in spatially confined systems.     

Midpoint cell method for hybrid (MPI+OpenMP) parallelization of molecular dynamics simulations

We have developed a new hybrid (MPI+OpenMP) parallelization scheme for molecular dynamics (MD) simulations by combining a cell‐wise version of the midpoint method with pair‐wise Verlet lists. In this scheme, which we call the midpoint cell method, simulation space is divided into subdomains, each of which is assigned to a MPI processor. Each subdomain is further divided into small cells. The interaction between two particles existing in different cells is computed in the subdomain containing the midpoint cell of the two cells where the particles reside. In each MPI processor, cell pairs are distributed over OpenMP threads for shared memory parallelization. The midpoint cell method keeps the advantages of the original midpoint method, while filtering out unnecessary calculations of midpoint checking for all the particle pairs by single midpoint cell determination prior to MD simulations. Distributing cell pairs over OpenMP threads allows for more efficient shared memory parallelization compared with distributing atom indices over threads. Furthermore, cell grouping of particle data makes better memory access, reducing the number of cache misses. The parallel performance of the midpoint cell method on the K computer showed scalability up to 512 and 32,768 cores for systems of 20,000 and 1 million atoms, respectively. One MD time step for long‐range interactions could be calculated within 4.5 ms even for a 1 million atoms system with particle‐mesh Ewald electrostatics. © 2014 Wiley Periodicals, Inc.     

Orthogonal sampling in free‐energy calculations of residue mutations in a tripeptide: TI versus λ‐LEUS

In a recent article (Bieler et al., J. Chem. Theory Comput. 2014, 10, 3006), we introduced a combination of λ‐dynamics and local‐elevation umbrella‐sampling termed λ‐LEUS to calculate free‐energy changes associated with alchemical processes using molecular dynamics simulations. This method was suggested to be more efficient than thermodynamic integration (TI), because the dynamical variation of the alchemical variable λ opens up pathways to circumvent barriers in the orthogonal space (defined by the N – 1 degrees of freedom that are not subjected to the sampling enhancement), a feature λ‐LEUS shares with Hamiltonian replica‐exchange (HR) approaches. However, the mutation considered, hydroquinone to benzene in water, was no real challenge in terms of orthogonal‐space properties, which were restricted to solvent‐relaxation processes. In the present article, we revisit the comparison between TI and λ‐LEUS considering non‐trivial mutations of the central residue X of a KXK tripeptide in water (with X = G, E, K, S, F, or Y). Side‐chain interactions that may include salt bridges, hydrogen bonds or steric clashes lead to slow relaxation in the orthogonal space, mainly in the two‐dimensional subspace spanned by the central and ψ dihedral angles of the peptide. The efficiency enhancement afforded by λ‐LEUS is confirmed in this more complex test system and can be attributed explicitly to the improved sampling of the orthogonal space. The sensitivity of the results to the nontrivial choices of a mass parameter and of a thermostat coupling time for the alchemical variable is also investigated, resulting in recommended ranges of 50 to 100 u nm² and 0.2 to 0.5 ps, respectively. © 2015 Wiley Periodicals, Inc.     

Ultrafast dynamics of myoglobin without the distal histidine: stimulated vibrational echo experiments and molecular dynamics simulations

Ultrafast protein dynamics of the CO adduct of a myoglobin mutant with the polar distal histidine replaced by a nonpolar valine (H64V) have been investigated by spectrally resolved infrared stimulated vibrational echo experiments and molecular dynamics (MD) simulations. In aqueous solution at room temperature, the vibrational dephasing rate of CO in the mutant is reduced by approximately 50% relative to the native protein. This finding confirms that the dephasing of the CO vibration in the native protein is sensitive to the interaction between the ligand and the distal histidine. The stimulated vibrational echo observable is calculated from MD simulations of H64V within a model in which vibrational dephasing is driven by electrostatic forces. In agreement with experiment, calculated vibrational echoes show slower dephasing for the mutant than for the native protein. However, vibrational echoes calculated for H64V do not show the quantitative agreement with measurements demonstrated previously for the native protein.     

Three pillars for achieving quantum mechanical molecular dynamics simulations of huge systems: Divide‐and‐conquer,density‐functional tight‐binding,and massively parallel computation

The linear‐scaling divide‐and‐conquer (DC) quantum chemical methodology is applied to the density‐functional tight‐binding (DFTB) theory to develop a massively parallel program that achieves on‐the‐fly molecular reaction dynamics simulations of huge systems from scratch. The functions to perform large scale geometry optimization and molecular dynamics with DC‐DFTB potential energy surface are implemented to the program called DC‐DFTB‐K. A novel interpolation‐based algorithm is developed for parallelizing the determination of the Fermi level in the DC method. The performance of the DC‐DFTB‐K program is assessed using a laboratory computer and the K computer. Numerical tests show the high efficiency of the DC‐DFTB‐K program, a single‐point energy gradient calculation of a one‐million‐atom system is completed within 60 s using 7290 nodes of the K computer. © 2016 Wiley Periodicals, Inc.     

Net charge changes in the calculation of relative ligand‐binding free energies via classical atomistic molecular dynamics simulation

The calculation of binding free energies of charged species to a target molecule is a frequently encountered problem in molecular dynamics studies of (bio‐)chemical thermodynamics. Many important endogenous receptor‐binding molecules, enzyme substrates, or drug molecules have a nonzero net charge. Absolute binding free energies, as well as binding free energies relative to another molecule with a different net charge will be affected by artifacts due to the used effective electrostatic interaction function and associated parameters (e.g., size of the computational box). In the present study, charging contributions to binding free energies of small oligoatomic ions to a series of model host cavities functionalized with different chemical groups are calculated with classical atomistic molecular dynamics simulation. Electrostatic interactions are treated using a lattice‐summation scheme or a cutoff‐truncation scheme with Barker–Watts reaction‐field correction, and the simulations are conducted in boxes of different edge lengths. It is illustrated that the charging free energies of the guest molecules in water and in the host strongly depend on the applied methodology and that neglect of correction terms for the artifacts introduced by the finite size of the simulated system and the use of an effective electrostatic interaction function considerably impairs the thermodynamic interpretation of guest‐host interactions. Application of correction terms for the various artifacts yields consistent results for the charging contribution to binding free energies and is thus a prerequisite for the valid interpretation or prediction of experimental data via molecular dynamics simulation. Analysis and correction of electrostatic artifacts according to the scheme proposed in the present study should therefore be considered an integral part of careful free‐energy calculation studies if changes in the net charge are involved. © 2013 The Authors Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Ab initio simulations of two‐dimensional electronic spectra: The SOS//QM/MM approach

Two‐dimensional electronic spectroscopy (2DES) is a cutting‐edge technique for investigating with high temporal resolution energy transfer, structure, and dynamics in a wide range of systems in physical chemistry, energy sciences, biophysics, and biocatalysis. However, the interpretation of 2DES is challenging and requires computational modeling. This perspective provides a roadmap for the development of computational tools that could be routinely applied to simulate 2DES spectra of multichromophoric systems active in the UV region (2DUV) using state‐of‐the‐art ab initio electronic structure methods within a quatum mechanics/molecular mechanics (QM/MM) scheme and the sum‐over‐states (SOS) approach (here called SOS//QM/MM). Multiconfigurational and multireference perturbative methods, such as the complete active space self‐consistent field and second‐order multireference perturbation theory (CASPT2) techniques, can be applied to reliably calculate the electronic properties of multichromophoric systems. Hybrid QM/MM method and molecular dynamics techniques can be used to assess environmental and conformational effects, respectively, that shape the 2D electronic spectra. DNA and proteins are important biological targets containing UV chromophores. We report ab initio simulation of 2DUV spectra of a cyclic tetrapeptide containing two interacting aromatic side chains, a model system for the study of protein structure and dynamics by means of 2DUV spectroscopy. © 2013 Wiley Periodicals, Inc.     

Deep desulfurization by oxidation using an active ionic liquid‐supported Zr metal–organic framework as catalyst

In this study, a Zr metal–organic framework (UIO‐66) was synthesized with zirconium tetrachloride and terephthalic acid using the solvent method. Then various masses of 1‐methylimidazolium‐3‐propylsulfonate hydrosulfate (PSMIMHSO₄) were supported on the UIO‐66 as catalysts, which were used for catalytic oxidative desulfurization. Sulfur removal using 400 mg of 40% PSMIMHSO₄ supported on the UIO‐66 of greater than 94% was obtained at 313 K for 20 min with an O/S molar ratio of 7:1. The results obtained in this work could provide useful information for the design of water‐stable metal–organic frameworks with permanent porosity in applications of catalytic oxidative desulfurization. Copyright © 2014 John Wiley & Sons, Ltd.