THz to far-infrared spectra of the known crystal polymorphs of phenylalanine.

There is renewed interest in the structure of the essential amino acid phenylalanine in the solid state. Three new polymorphs were found in the years 2012 to 2014. Here, we investigate the structure, stability, and energetical ordering of these phases using first-principles simulations at the level of density functional theory incorporating van der Waals interactions. Two of the distinct crystal forms are found to be structurally similar and energetically very close after vibrational free energy corrections have been taken into account. Infrared absorption spectra are likewise calculated and compared to experimental measurements. By combining measurements obtained with a commercial Fourier transform infra-red spectrometer and a homemade air-photonics-based THz time domain spectrometer, we could carry out this comparison in the vibrational frequency region from 1 to 40 THz. The excellent agreement of the line positions and the established energy ranking allow us to identify the most stable polymorph of phenylalanine.

Exact non-adiabatic coupling vectors for the time-dependent density functional based tight-binding method.

We report on non-adiabatic coupling vectors between electronic excited states for the time-dependent-density functional theory based tight-binding (TD-DFTB) method. The implementation includes orbital relaxation effects that have been previously neglected and covers also the case of range-separated exchange-correlation functionals. Benchmark calculations with respect to first principles TD-DFT highlight the large dependence of non-adiabatic couplings on the functional. Closer investigations of the topology around a conical intersection between excited states show that TD-DFTB delivers near-exact values of the Berry phase, which paves the way for consistent non-adiabatic molecular dynamics simulations for large systems.

Molecular electronic refrigeration against parallel phonon heat leakage channels.

Due to their structured density of states, molecular junctions provide rich resources to filter and control the flow of electrons and phonons. Here we compute the out of equilibrium current-voltage characteristics and dissipated heat of some recently synthesized oligophenylenes (OPE3) using the Density Functional based Tight-Binding (DFTB) method within Non-Equilibrium Green's Function Theory (NEGF). We analyze the Peltier cooling power for these molecular junctions as function of a bias voltage and investigate the parameters that lead to optimal cooling performance. In order to quantify the attainable temperature reduction, an electro-thermal circuit model is presented, in which the key electronic and thermal transport parameters enter. Overall, our results demonstrate that the studied OPE3 devices are compatible with temperature reductions of several K. Based on the results, some strategies to enable high performance devices for cooling applications are briefly discussed.

Pathways to fluorescence via restriction of intramolecular motion in substituted tetraphenylethylenes.

The design of materials with enhanced luminescence properties is a fast-developing field due to the potential applicability of these materials as light-emitting diodes or for bioimaging. A transparent way to enhance the emission properties of interesting molecular candidates is blocking competing and unproductive non-radiative relaxation pathways by the restriction of intramolecular motions. Rationalized functionalization is an important possibility to achieve such restrictions. Using time-dependent density functional theory (TD-DFT) based on the ωB97XD functional and the semiempirical tight-binding method including long-range corrections (TD-LC-DFTB), this work investigates the effect of functionalization of the paradigmatic tetraphenylethylene (TPE) on achieving restricted access to conical intersections (RACI). Photodynamical surface hopping simulations have been performed on a larger set of compounds including TPE and ten functionalized TPE compounds. Functionalization has been achieved by means of electron-withdrawing groups, bulky groups which block the relaxation channels via steric hindrance and groups capable of forming strong hydrogen bonds, which restrict the motion via the formation of hydrogen bond channels. Most of the investigated functionalized TPE candidates show ultrafast deactivation to the ground state due to their still existing structural flexibility, but two examples, one containing -CN and -CF3 groups and a second characterized by a network of hydrogen bonds, have been identified as interesting candidates for creating efficient luminescence properties in solution.

P∩N Bridged Cu(I) Dimers Featuring Both TADF and Phosphorescence. From Overview towards Detailed Case Study of the Excited Singlet and Triplet States.

We present an overview over eight brightly luminescent Cu(I) dimers of the type Cu2X2(P∩N)3 with X = Cl, Br, I and P∩N = 2-diphenylphosphino-pyridine (Ph2Ppy), 2-diphenylphosphino-pyrimidine (Ph2Ppym), 1-diphenylphosphino-isoquinoline (Ph2Piqn) including three new crystal structures (Cu2Br2(Ph2Ppy)3 1-Br, Cu2I2(Ph2Ppym)3 2-I and Cu2I2(Ph2Piqn)3 3-I). However, we mainly focus on their photo-luminescence properties. All compounds exhibit combined thermally activated delayed fluorescence (TADF) and phosphorescence at ambient temperature. Emission color, decay time and quantum yield vary over large ranges. For deeper characterization, we select Cu2I2(Ph2Ppy)3, 1-I, showing a quantum yield of 81%. DFT and SOC-TDDFT calculations provide insight into the electronic structures of the singlet S1 and triplet T1 states. Both stem from metal+iodide-to-ligand charge transfer transitions. Evaluation of the emission decay dynamics, measured from 1.2 ≤ T ≤ 300 K, gives ∆E(S1-T1) = 380 cm-1 (47 meV), a transition rate of k(S1→S0) = 2.25 × 106 s-1 (445 ns), T1 zero-field splittings, transition rates from the triplet substates and spin-lattice relaxation times. We also discuss the interplay of S1-TADF and T1-phosphorescence. The combined emission paths shorten the overall decay time. For OLED applications, utilization of both singlet and triplet harvesting can be highly favorable for improvement of the device performance.

Analytical Time-Dependent Long-Range Corrected Density Functional Tight Binding (TD-LC-DFTB) Gradients in DFTB+: Implementation and Benchmark for Excited-State Geometries and Transition Energies.

The absorption and emission of light is a ubiquitous process in chemical and biological processes, making a theoretical description inevitable for understanding and predicting such properties. Although ab initio and DFT methods are capable of describing excited states with good accuracy in many cases, the investigation of dynamical processes and the need to sample the phase space in complex systems often requires methods with reduced computational costs but still sufficient accuracy. In the present work, we report the derivation and implementation of analytical nuclear gradients for time-dependent long-range corrected density functional tight binding (TD-LC-DFTB) in the DFTB+ program. The accuracy of the TD-LC-DFTB potential-energy surfaces is benchmarked for excited-state geometries and adiabatic as well as vertical transition energies. The benchmark set consists of more than 100 organic molecules taken as subsets from available benchmark sets. The reported method yields a mean deviation of 0.31 eV for adiabatic excitation energies with respect to CC2. In order to study more subtle effects, seminumerical second derivatives based on the analytical gradients are employed to simulate vibrationally resolved UV/vis spectra. This extensive test exhibits few problematic cases, which can be traced back to the parametrization of the repulsive potential.

Excitonic and charge transfer interactions in tetracene stacked and T-shaped dimers.

Extended quantum chemical calculations were performed for the tetracene dimer to provide benchmark results, analyze the excimer survival process, and explore the possibility of using long-range-corrected (LC) time-dependent second-order density functional tight-biding (DFTB2) for this system. Ground- and first-excited-state optimized geometries, vertical excitations at relevant minima, and intermonomer displacement potential energy curves (PECs) were calculated for these purposes. Ground-state geometries were optimized with the scaled-opposite-spin (SOS) second-order Møller-Plesset perturbation (MP2) theory and LC-DFT (density functional theory) and LC-DFTB2 levels. Excited-state geometries were optimized with SOS-ADC(2) (algebraic diagrammatic construction to second-order) and the time-dependent approaches for the latter two methods. Vertical excitations and PECs were compared to multireference configuration interaction DFT (DFT/MRCI). All methods predict the lowest-energy S0 conformer to have monomers parallel and rotated relative to each other and the lowest S1 conformer to be of a displaced-stacked type. LC-DFTB2, however, presents some relevant differences regarding other conformers for S0. Despite some state-order inversions, overall good agreement between methods was observed in the spectral shape, state character, and PECs. Nevertheless, DFT/MRCI predicts that the S1 state should acquire a doubly excited-state character relevant to the excimer survival process and, therefore, cannot be completely described by the single reference methods used in this work. PECs also revealed an interesting relation between dissociation energies and the intermonomer charge-transfer interactions for some states.

Simulation of Structural Evolution Using Time-Dependent Density-Functional Based Tight-Binding Method.

A faster and more efficient quantum mechanical simulation method for application to complicated issues of real systems beyond model cases has long been sought after. The density-functional based tight-binding (DFTB) method has successfully explained the atomistic and electronic properties of semiconductors, surfaces, and nanostructures. In addition, the time-dependent formalism implemented in DFTB showed high efficiency in terms of computational cost. In this study, we demonstrated the structural and electronic evolution of small molecules induced by a laser pulse using the time-dependent DFTB (TD-DFTB) method. We identified the critical fluence of the input laser for structural dissociations in carbon chains and fullerenes, which related to the structural stability. The excitation energies of several molecules calculated by TD-DFTB agreed with the experimental values.

Towards a simplified description of thermoelectric materials: accuracy of approximate density functional theory for phonon dispersions.

We calculate the phonon-dispersion relations of several two-dimensional materials and diamond using the density-functional based tight-binding approach (DFTB). Our goal is to verify if this numerically efficient method provides sufficiently accurate phonon frequencies and group velocities to compute reliable thermoelectric properties. To this end, the results are compared to available DFT results and experimental data. To quantify the accuracy for a given band, a descriptor is introduced that summarizes contributions to the lattice conductivity that are available already in the harmonic approximation. We find that the DFTB predictions depend strongly on the employed repulsive pair-potentials, which are an important prerequisite of this method. For carbon-based materials, accurate pair-potentials are identified and lead to errors of the descriptor that are of the same order as differences between different local and semi-local DFT approaches.

The Fragment Molecular Orbital Method Based on Long-Range Corrected Density-Functional Tight-Binding.

The presently available linear scaling approaches to density-functional tight-binding (DFTB) based on the fragment molecular orbital (FMO) method are severely impacted by the problem of artificial charge transfer due to the self-interaction error (SIE), which hampers the simulation of zwitterionic systems such as biopolymers or ionic liquids. Here we report an extension of FMO-DFTB where we included a long-range corrected (LC) functional designed to mitigate the DFTB SIE, called the FMO-LC-DFTB method, resulting in a robust method which succeeds in simulating zwitterionic systems. Both energy and analytic gradient are developed for the gas phase and the polarizable continuum model of solvation. The scaling of FMO-LC-DFTB with system size N is shown to be almost linear, O( N1.13-1.28), and its numerical accuracy is established for a variety of representative systems including neutral and charged polypeptides. It is shown that pair interaction energies between fragments for two mini-proteins are in excellent agreement with results from long-range corrected density functional theory. The new method was employed in long time scale (1 ns) molecular dynamics simulations of the tryptophan cage protein (PDB: 1L2Y ) in the gas phase for four different protonation states and in stochastic global minimum structure searches for 1-ethyl-3-methylimidazolium nitrate ionic liquid clusters containing up to 2300 atoms.

Ultrafast Nonadiabatic Cascade and Subsequent Photofragmentation of Extreme Ultraviolet Excited Caffeine Molecule.

Ultrafast XUV chemistry is offering new opportunities to decipher the complex dynamics taking place in highly excited molecular states and thus better understand fundamental natural phenomena as molecule formation in interstellar media. We used ultrashort XUV light pulses to perform XUV pump-IR probe experiments in caffeine as a model of prebiotic molecule. We observed a 40 fs decay of excited cationic states. Guided by quantum calculations, this time scale is interpreted in terms of a nonadiabatic cascade through a large number of highly correlated states. This shows that the correlation driven nonadiabatic relaxation seems to be a general process for highly excited states, which might impact our understanding of molecular processing in interstellar media.

New candidates for the global minimum of medium-sized silicon clusters: A hybrid DFTB/DFT genetic algorithm applied to Si n , n = 8-80.

In this study, we perform a systematic search to find the possible lowest energy structure of silicon nanoclusters Sin (n = 8-80) by means of an evolutionary algorithm. The fitness function for this search is the total energy of density functional tight binding (DFTB). To be on firm ground, we take several low energy structures of DFTB and perform further geometrical optimization by density functional theory (DFT). Then we choose structures with the lowest DFT total energy and compare them with the reported lowest energy structures in the literature. In our search, we found several lowest energy structures that were previously unreported. We further observe a geometrical transition at n = 27 from elongated to globular structures. In addition, the optical gap of the lowest energy structures is investigated by time-dependent DFTB (TD-DFTB) and time-dependent DFT (TD-DFT). The results show the same trend in TD-DFTB and TD-DFT for the optical gap. We also find a sudden drop in the optical gap at n = 27, precisely where the geometrical transition occurs.

An efficient approximate algorithm for nonadiabatic molecular dynamics.

We propose a modification to the nonadiabatic surface hopping calculation method formulated in a paper by Yu et al. [Phys. Chem. Chem. Phys. 16, 25883 (2014)], which is a multidimensional extension of the Zhu-Nakamura theory with a practical diabatic gradient estimation algorithm. In our modification, their diabatic gradient estimation algorithm, which is based on a simple interpolation of the adiabatic potential energy surfaces, is replaced by an algorithm using the numerical derivatives of the adiabatic gradients. We then apply the algorithm to several models of nonadiabatic dynamics, both analytic and ab initio models, to numerically demonstrate that our method indeed widens the applicability and robustness of their method. We also discuss the validity and limitations of our new nonadiabatic surface hopping method while considering in mind potential applications to excited-state dynamics of biomolecules or unconventional nonadiabatic dynamics such as radiation decay processes in ultraintense X-ray fields.

Parametrization and Benchmark of Long-Range Corrected DFTB2 for Organic Molecules.

We present the parametrization and benchmark of long-range corrected second-order density functional tight binding (DFTB), LC-DFTB2, for organic and biological molecules. The LC-DFTB2 model not only improves fundamental orbital energy gaps but also ameliorates the DFT self-interaction error and overpolarization problem, and further improves charge-transfer excited states significantly. Electronic parameters for the construction of the DFTB2 Hamiltonian as well as repulsive potentials were optimized for molecules containing C, H, N, and O chemical elements. We use a semiautomatic parametrization scheme based on a genetic algorithm. With the new parameters, LC-DFTB2 describes geometries and vibrational frequencies of organic molecules similarly well as third-order DFTB3/3OB, the de facto standard parametrization based on a GGA functional. LC-DFTB2 performs well also for atomization and reaction energies, however, slightly less satisfactorily than DFTB3/3OB.

Nonadiabatic Dynamics of Cycloparaphenylenes with TD-DFTB Surface Hopping.

We implemented a version of the decoherence-corrected fewest switches surface hopping based on linear-response time-dependent density functional tight binding (TD-DFTB), enhanced by transition density analysis. The method has been tested for the gas-phase relaxation dynamics of two cycloparaphenylene molecules, [8]CPP and [10]CPP, explaining some important features of their nonadiabatic dynamics, such as the origin of their long fluorescence lifetimes (related to the slow radiative emission from the S1 state) and the trend of increasing the fluorescence rate with the molecular size (related to an increase in the S1-S0 energy gaps and oscillator strengths in the larger molecule). The quality of the TD-DFTB electronic structure information was assessed through four quantities: excitation energies; charge-transfer (CT) numbers, which estimate the charge transfer character of states; participation ratio (PR), which describes delocalization of electronic density; and participation ratio of natural transition orbitals (PRNTO), which describes the multiconfigurational character of states. These quantities were computed during dynamics and recomputed for the same geometries with the higher-level long-range-corrected TD-LC-DFTB and a lower-level single-determinant approximation for the excited states, SD-(LC)-DFTB. Taking TD-LC-DFTB as the standard, TD-DFTB underestimates the excitation energies by ∼0.5 eV and overestimates CT and PR. SD-DFTB underestimates excitation energies and overestimates CT to the same extent that TD-DFTB does, but it predicts reasonable PR distributions. SD-LC-DFTB leads to an extreme overestimation of the excitation energies by ∼3 eV, overestimates the charge transfer character of the state, but predicts the PR values very close to those obtained with TD-LC-DFTB.

Charge transfer excitations from particle-particle random phase approximation-Opportunities and challenges arising from two-electron deficient systems.

The particle-particle random phase approximation (pp-RPA) is a promising method for studying charge transfer(CT) excitations. Through a detailed analysis on two-electron deficient systems, we show that the pp-RPA is always able to recover the long-distance asymptotic -1/R trend for CT excitations as a result of the concerted effect between orbital energies and the pp-RPA kernel. We also provide quantitative results for systems with relatively short donor-acceptor distances. With conventional hybrid or range-separated functionals, the pp-RPA performs much better than time-dependent density functional theory (TDDFT), although it still gives underestimated results which are not as good as TDDFT with system-dependent tuned functionals. For pp-RPA, there remain three great challenges in dealing with CT excitations. First, the delocalized frontier orbitals in strongly correlated systems often lead to difficulty with self-consistent field convergence as well as an incorrect picture with about half an electron transferred. Second, the commonly used density functionals often underestimate the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital (LUMO) for the two-electron deficient species, resulting in systems with delocalized orbitals. Third, the performance of pp-RPA greatly depends on the energy difference between the LUMO and a higher virtual orbital. However, the meaning of the orbital energies for higher virtual orbitals is still not clear. We also discuss the performance of an approximate pp-RPA scheme that uses density functional tight binding (pp-DFTB) as reference and demonstrate that the aforementioned challenges can be overcome by adopting suitable range-separated hybrid functionals. The pp-RPA and pp-DFTB are thus promising general approaches for describing charge transfer excitations.

Time-Dependent Extension of the Long-Range Corrected Density Functional Based Tight-Binding Method.

We present a consistent linear response formulation of the density functional based tight-binding method for long-range corrected exchange-correlation functionals (LC-DFTB). Besides a detailed account of derivation and implementation of the method, we also test the new scheme on a variety of systems considered to be problematic for conventional local/semilocal time-dependent density functional theory (TD-DFT). To this class belong the optical properties of polyacenes and nucleobases, as well as charge transfer excited states in molecular dimers. We find that the approximate LC-DFTB method exhibits the same general trends and similar accuracy as range-separated DFT methods at significantly reduced computational cost. The scheme should be especially useful in the determination of the electronic excited states of very large molecules, for which conventional TD-DFT is supposed to fail due to a multitude of artificial low energy states.

Dynamic Characteristics of Aggregation Effects of Organic Dyes in Dye-Sensitized Solar Cells.

Two organic dyes (LS-1 and IQ4) containing identical electron donor and acceptor units but distinct π units result in significantly different power conversion efficiency of the corresponding dye-sensitized solar cells (DSSCs): LS-1, 4.4%, and IQ4, 9.2%. Herein, we combine first-principle calculations and molecular dynamics to explore the aggregation effects of LS-1 and IQ4 by comparing their optical properties and intermolecular electronic couplings. The calculated absorption spectra are in good agreement with the experimental observations and reveal them to be evidently affected by the dimerization. Furthermore, molecular dynamics simulations show that steric hindrance induced by the diphenylquinoxaline unit in IQ4 can elongate the distances between intermolecular π units or electron donors, which are responsible for the fact that the intermolecular electronic coupling of LS-1 is about 10 times larger than that of IQ4. More importantly, the aggregated IQ4 remains almost perpendicular to the TiO2 surface, whereas LS-1 gradually tilts during the dynamic simulation, impacting electron injection and recombination in several ways, which clarifies why IQ4 leads to larger photocurrent and higher conversion efficiency. The deep understanding of the dye aggregation effects sheds new light on the complex factors determining DSSC function and paves the way for rational design of high-efficiency self-anti-aggregation sensitizers.

Efficient emission facilitated by multiple energy level transitions in uniform graphitic carbon nitride films deposited by thermal vapor condensation.

Graphitic carbon nitride (g-CN) films are important components of optoelectronic devices, but current techniques for their production, such as drop casting and spin coating, fail to deliver uniform and pinhole-free g-CN films on solid substrates. Here, versatile, cost-effective, and large-area growth of uniform and pinhole-free g-CN films is achieved by using a thermal vapor condensation method under atmospheric pressure. A comparison of the X-ray diffraction and Fourier transform infrared data with the calculated infrared spectrum confirmed the graphitic build-up of films composed of tri-s-triazine units. These g-CN films possess multiple active energy states including π*, π, and lone-pair states, which facilitate their efficient (6% quantum yield in the solid state) photoluminescence, as confirmed by both experimental measurements and theoretical calculations.

Signatures in vibrational and UV-visible absorption spectra for identifying cyclic hydrocarbons by graphene fragments.

To promote possible applications of graphene in molecular identification based on stacking effects, in particular in recognizing aromatic amino acids and even sequencing nucleobases in life sciences, we comprehensively study the interaction between graphene segments and different cyclic organic hydrocarbons including benzene (C6H6), cyclohexane (C6H12), benzyne (C6H4), cyclohexene (C6H10), 1,3-cyclohexadiene (C6H8(1)) and 1,4-cyclohexadiene (C6H8(2)), using the density-functional tight-binding (DFTB) method. Interestingly, we find obviously different characteristics in Raman vibrational and ultraviolet visible absorption spectra of the small molecules adsorbed on the graphene sheet. Specifically, we find that both spectra involve clearly different characteristic peaks, belonging to the different small molecules upon adsorption, with the ones of ionized molecules being more substantial. Further analysis shows that the adsorptions are almost all due to the presence of dispersion energy in neutral cases and involve charge transfer from the graphene to the small molecules. In contrast, the main binding force in the ionic adsorption systems is the electronic interaction. The results present clear signatures that can be used to recognize different kinds of aromatic hydrocarbon rings on graphene sheets. We expect that our findings will be helpful for designing molecular recognition devices using graphene.