Johannes Neugebauer

Quantum chemical calculation of vibrational spectra of large molecules—Raman and IR spectra for Buckminsterfullerene

In this work we demonstrate how different modern quantum chemical methods can be efficiently combined and applied for the calculation of the vibrational modes and spectra of large molecules. We are aiming at harmonic force fields, and infrared as well as Raman intensities within the double harmonic approximation, because consideration of higher order terms is only feasible for small molecules. In particular, density functional methods have evolved to a powerful quantum chemical tool for the determination of the electronic structure of molecules in the last decade. Underlying theoretical concepts for the calculation of intensities are reviewed, emphasizing necessary approximations and formal aspects of the introduced quantities, which are often not explicated in detail in elementary treatments of this topic. It is shown how complex quantum chemistry program packages can be interfaced to new programs in order to calculate IR and Raman spectra. The advantages of numerical differentiation of analytical gradients, dipole moments, and static, as well as dynamic polarizabilities, are pointed out. We carefully investigate the influence of the basis set size on polarizabilities and their spatial derivatives. This leads us to the construction of a hybrid basis set, which is equally well suited for the calculation of vibrational frequencies and Raman intensities. The efficiency is demonstrated for the highly symmetric C60, for which we present the first all‐electron density functional calculation of its Raman spectrum.

Fundamental vibrational frequencies of small polyatomic molecules from density-functional calculations and vibrational perturbation theory

An extensive study of fundamental frequencies and anharmonic vibrational constants for polyatomic molecules obtained from Becke three parameter Lee–Yang–Parr (B3LYP) and Becke–Perdew (BP86) density functional calculations is presented. These calculations are based on standard perturbation theory, and are compared to correlation-corrected vibrational self-consistent field (CC-VSCF) calculations for the water dimer. The anharmonic corrections obtained from density-functional calculations compare well with experimental values and with results from correlated ab initio methods. While fundamental frequencies from B3LYP calculations are reliable, they are considerably too small for BP86 calculations. Consequently, the good agreement of unscaled harmonic frequencies from BP86 calculations with experimental frequencies is due to an error cancellation effect. This is of importance for the prediction of vibrational spectra for large molecules, because the perturbation theory approach naturally becomes unreliable fo...

The merits of the frozen-density embedding scheme to model solvatochromic shifts

We investigate the usefulness of a frozen-density embedding scheme within density-functional theory [J. Phys. Chem. 97, 8050 (1993)] for the calculation of solvatochromic shifts. The frozen-density calculations, particularly of excitation energies have two clear advantages over the standard supermolecule calculations: (i) calculations for much larger systems are feasible, since the time-consuming time-dependent density functional theory (TDDFT) part is carried out in a limited molecular orbital space, while the effect of the surroundings is still included at a quantum mechanical level. This allows a large number of solvent molecules to be included and thus affords both specific and nonspecific solvent effects to be modeled. (ii) Only excitations of the system of interest, i.e., the selected embedded system, are calculated. This allows an easy analysis and interpretation of the results. In TDDFT calculations, it avoids unphysical results introduced by spurious mixings with the artificially too low charge-transfer excitations which are an artifact of the adiabatic local-density approximation or generalized gradient approximation exchange-correlation kernels currently used. The performance of the frozen-density embedding method is tested for the well-studied solvatochromic properties of the n-->pi(*) excitation of acetone. Further enhancement of the efficiency is studied by constructing approximate solvent densities, e.g., from a superposition of densities of individual solvent molecules. This is demonstrated for systems with up to 802 atoms. To obtain a realistic modeling of the absorption spectra of solvated molecules, including the effect of the solvent motions, we combine the embedding scheme with classical molecular dynamics (MD) and Car-Parrinello MD simulations to obtain snapshots of the solute and its solvent environment, for which then excitation energies are calculated. The frozen-density embedding yields estimated solvent shifts in the range of 0.20-0.26 eV, in good agreement with experimental values of between 0.19 and 0.21 eV.

Subsystem density‐functional theory

Subsystem density‐functional theory (subsystem DFT) has developed into a powerful alternative to Kohn–Sham DFT for quantum chemical calculations of complex systems. It exploits the idea of representing the total electron density as a sum of subsystem densities. The optimum total density is found by minimizing the total energy with respect to each of the subsystem densities, which breaks down the electronic‐structure problem into effective subsystem problems. This enables calculations on large molecular aggregates and even (bio‐)polymers without system‐specific parameterizations. We provide a concise review of the underlying theory, typical approximations, and embedding approaches related to subsystem DFT such as frozen‐density embedding (FDE). Moreover, we discuss extensions and applications of subsystem DFT and FDE to molecular property calculations, excited states, and wave function in DFT embedding methods. Furthermore, we outline recent developments for reconstruction techniques of embedding potentials arising in subsystem DFT, and for using subsystem DFT to incorporate constraints into DFT calculations.

Accurate frozen-density embedding potentials as a first step towards a subsystem description of covalent bonds

The frozen-density embedding (FDE) scheme [Wesolowski and Warshel, J. Phys. Chem. 97, 8050 (1993)] relies on the use of approximations for the kinetic-energy component v(T)[rho(1),rho(2)] of the embedding potential. While with approximations derived from generalized-gradient approximation kinetic-energy density functional weak interactions between subsystems such as hydrogen bonds can be described rather accurately, these approximations break down for bonds with a covalent character. Thus, to be able to directly apply the FDE scheme to subsystems connected by covalent bonds, improved approximations to v(T) are needed. As a first step toward this goal, we have implemented a method for the numerical calculation of accurate references for v(T). We present accurate embedding potentials for a selected set of model systems, in which the subsystems are connected by hydrogen bonds of various strength (water dimer and F-H-F(-)), a coordination bond (ammonia borane), and a prototypical covalent bond (ethane). These accurate potentials are analyzed and compared to those obtained from popular kinetic-energy density functionals.

A flexible implementation of frozen‐density embedding for use in multilevel simulations

A new implementation of frozen‐density embedding (FDE) in the Amsterdam Density Functional (ADF) program package is presented. FDE is based on a subsystem formulation of density‐functional theory (DFT), in which a large system is assembled from an arbitrary number of subsystems, which are coupled by an effective embedding potential. The new implementation allows both an optimization of all subsystems as a linear‐scaling alternative to a conventional DFT treatment, the calculation of one active fragment in the presence of a frozen environment, and intermediate setups, in which individual subsystems are fully optimized, partially optimized, or completely frozen. It is shown how this flexible setup can facilitate the application of FDE in multilevel simulations.

A mode-selective quantum chemical method for tracking molecular vibrations applied to functionalized carbon nanotubes

The theoretical investigation of mesoscopic objects requires new techniques which are particularly suited for the study of selected aspects of these systems. Vibrational spectroscopy is a main source for structural information on heterogeneous systems. We present an efficient quantum chemical method, which relies on a modified Davidson algorithm for targeting selected vibrations in infrared and Raman spectra. This approach is applied to the characteristic breathing modes of single-walled carbon nanotubes.

Gas-Phase C−H and N−H Bond Activation by a High Valent Nitrido-Iron Dication and 〈NH〉-Transfer to Activated Olefins

A tetrapodal pentadentate nitrogen ligand (2,6-bis(1,1-di(aminomethyl)ethyl)pyridine, 1) is used for the synthesis of the azido-iron(III) complex [(1)Fe(N3)]X2 where X is either Br or PF6. By means of electrospray ionization mass spectrometry, the dication [(1)Fe(N3)]2+ can be transferred into the gas phase as an intact entity. Upon collisional activation, [(1)Fe(N3)]2+ undergoes an expulsion of molecular nitrogen to afford the dicationic nitrido-iron species [(1)FeN]2+ as an intermediate, which upon further activation can intramolecularly activate C-H- and N-H bonds of the chelating ligand 1 or can transfer an NH unit in bimolecular reactions with activated olefins. The precursor dication [(1)Fe(N3)]2+, the resulting nitrido species [(1)FeN]2+, and its possible isomers are investigated by mass spectrometric experiments, isotopic labeling, and complementary computational studies using density functional theory.

Quantum Chemical Description of Absorption Properties and Excited‐State Processes in Photosynthetic Systems

The theoretical description of the initial steps in photosynthesis has gained increasing importance over the past few years. This is caused by more and more structural data becoming available for light-harvesting complexes and reaction centers which form the basis for atomistic calculations and by the progress made in the development of first-principles methods for excited electronic states of large molecules. In this Review, we discuss the advantages and pitfalls of theoretical methods applicable to photosynthetic pigments. Besides methodological aspects of excited-state electronic-structure methods, studies on chlorophyll-type and carotenoid-like molecules are discussed. We also address the concepts of exciton coupling and excitation-energy transfer (EET) and compare the different theoretical methods for the calculation of EET coupling constants. Applications to photosynthetic light-harvesting complexes and reaction centers based on such models are also analyzed.

Subsystem-Based Theoretical Spectroscopy of Biomolecules and Biomolecular Assemblies

The absorption properties of chromophores in biomolecular systems are subject to several fine-tuning mechanisms. Specific interactions with the surrounding protein environment often lead to significant changes in the excitation energies, but bulk dielectric effects can also play an important role. Moreover, strong excitonic interactions can occur in systems with several chromophores at close distances. For interpretation purposes, it is often desirable to distinguish different types of environmental effects, such as geometrical, electrostatic, polarization, and response (or differential polarization) effects. Methods that can be applied for theoretical analyses of such effects are reviewed herein, ranging from continuum and point-charge models to explicit quantum chemical subsystem methods for environmental effects. Connections to physical model theories are also outlined. Prototypical applications to optical spectra and excited states of fluorescent proteins, biomolecular photoreceptors, and photosynthetic protein complexes are discussed.