a Photosynthesis includes capturing sunlight by an assembly of molecules, called chlorophylls, and directing the harvested energy in the form of electronic excitations to the reaction center. Here we report, using realspace density functional theory and time-dependent density functional theory together with GW calculations, the optical and electronic properties of the two main chlorophylls in green plants, namely, chlorophylls a and b. Furthermore, we estimate the dipole and primitive quadrupole electric moments of these molecules. We employ Casida's assignment ansatz to study the absorption spectra of the chlorophylls in the two main red and blue regions at various environments with different exchangecorrelation functionals. In addition, we obtain the band gap of chlorophylls a and b, which are all in remarkable agreement with experimental observations.
Density functional theory (DFT) including van der Waals (vdW) interactions and accounting for zero-point energy (ZPE) is believed to provide a good description of crystalline ice phases [B. Pamuk et al., Phys. Rev. Lett. 108, 193003 (2012)]. Given the computational cost of DFT, it is not surprising that extensive phonon calculations, which yield the ZPE, have only been done for a limited amount of ice structures. Computationally convenient force fields on the other hand are the method of choice for large systems and/or dynamical simulations, e.g., of supercooled water. Here, we present a systematic comparison for seven hydrogen-ordered crystalline ice phases (Ih, IX, II, XIII, XIV, XV, and VIII) between many commonly used nonpolarizable force fields and density functionals, including some recently developed meta-GGA functionals and accounting for vdW interactions. Starting from the experimentally determined crystal structures, we perform space-group-constrained structural relaxations. These provide the starting point for highly accurate phonon calculations that yield effectively volume-dependent ZPEs within the quasiharmonic approximation. In particular, when including ZPE, the force fields show a remarkably good performance for equilibrium volumes and cohesive energies superior to many density functionals. A decomposition of the cohesive energies into intramolecular deformation, electrostatic, and vdW contributions quantifies the differences between force fields and DFT. Results for the equilibrium volumes and phase transition pressures for all studied force fields are much more strongly affected by ZPE than all studied density functionals. We track this down to significantly smaller shifts of the O-H-stretch modes and compare with experimental data from Raman spectroscopy.
A potential function is presented for describing a system of flexible H 2 O molecules based on the single-center multipole expansion (SCME) of the electrostatic interaction. The model, referred to as SCME/f, includes the variation of the molecular quadrupole moment as well as the dipole moment with changes in bond length and angle so as to reproduce results of high-level electronic structure calculations. The multipole expansion also includes fixed octupole and hexadecapole moments, as well as anisotropic dipole−dipole, dipole−quadrupole, and quadrupole−quadrupole polarizability tensors. The model contains five adjustable parameters related to the repulsive interaction and damping functions in the electrostatic and dispersion interactions. Their values are adjusted to reproduce the lowest energy isomers of small clusters, (H 2 O) n with n = 2−6, as well as measured properties of the ice Ih crystal. Subsequent calculations of the energy difference between the various isomer configurations of the clusters show that SCME/f gives good agreement with results of electronic structure calculations and represents a significant improvement over the previously presented rigid SCME potential function. Analysis of the vibrational frequencies of the clusters and structural properties of ice Ih crystal show the importance of accurately describing the variation of the quadrupole moment with molecular structures.
Calorimetric
studies on ice II reveal a surprising H
2
O/D
2
O isotope effect. While the ice II to ice Ic transition
is endothermic for H
2
O, it is exothermic for D
2
O samples. The transition enthalpies are +40 and −140 J/mol,
respectively, where such a sign change upon isotope substitution is
unprecedented in ice research. To understand the observations we employ
force field calculations using two water models known to perform well
for H
2
O ice phases and their vibrational properties. These
simulations reveal that the isotope effect can be traced back to zero-point
energy. q-TIP4P/F fares better and is able to account for approximately
three-fourths of the isotope effect, while MB-pol only catches approximately
one-third. Phonon and configurational entropy contributions are necessary
to predict reasonable transition enthalpies, but they do not have
an impact on the isotope effect. We suggest to use these calorimetric
isotope data as a benchmark for water models.
The anomalous volume isotope effect (VIE) of ice Ih is
calculated
and analyzed based on the quasi-harmonic approximation to account
for nuclear quantum effects in the Helmholtz free energy. While a
lot of recently developed polarizable many-body potential functions
give a normal VIE contrary to experimental results, we find that one
of them, MB-pol, yields the anomalous VIE in good agreement with the
most recent high-resolution neutron diffraction measurementsbetter
than DFT calculations. The short-range three-body terms in the MB-pol
function, which are fitted to CCSD(T) calculations, are found to have
a surprisingly large influence. A vibrational mode group decomposition
of the zero-point pressure together with a hitherto unconsidered benchmark
value for the intramolecular stretching modes of H2O ice
Ih obtained from Raman spectroscopy data unveils the reason for the
VIE: a delicate competition between the latter and the librations.
A potential function describing a system of flexible water molecules based on a single center multipole expansion of the electrostatic interactions is described, denoted f-SCME. The potential function includes a quadrupole moment surface (QMS) that reproduces results of high level configuration interaction calculations in addition to the commonly used dipole moment surface (DMS) developed by Partridge and Schwenke.The use of the so-called M-site models based on the DMS atomic charges to represent the QMS is explored, and some improvements presented. The potential function also includes the static octupole and hexadecapole moments and anisotropic dipole-dipole, dipole-quadrupole and quadrupole-quadrupole polarizability tensors as well as dispersion interaction of the original rigid SCME potential [SCME, Wikfeldt et al, PCCP
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.