B3LYP Research Articles

Which Electronic Structure Method to Choose in Trajectory Surface Hopping Dynamics Simulations? Azomethane as a Case Study.

Nonadiabatic dynamics simulations have become a standard approach to explore photochemical reactions. Such simulations require underlying potential energy surfaces and couplings between them, calculated at a chosen level of theory, yet this aspect is rarely assessed. Here, in combination with the popular trajectory surface hopping dynamics method, we use a high-accuracy XMS-CASPT2 electronic structure level as a benchmark for assessing the performances of various post-Hartree-Fock methods (namely, CIS, ADC(2), CC2, and CASSCF) and exchange-correlation functionals (PBE, PBE0, and CAM-B3LYP) in a TD-DFT/TDA context, using the isomerization around a double bond as test case. Different relaxation pathways are identified, and the ability of the different methods to reproduce their relative importance and time scale is discussed. The results show that multireference electronic structure methods should be preferred, when studying nonadiabatic decay between excited and ground states. If not affordable, TD-DFT with TDA and hybrid functionals and ADC(2) are efficient alternatives but overestimate the nonradiative decay yield and thus may miss deexcitation pathways.

Inquiring 199Hg NMR Parameters by Combining Ab Initio Molecular Dynamics and Relativistic NMR Calculations.

Ab initio molecular dynamics (AIMD) sampling followed by relativistic density functional theory (DFT) 199Hg NMR calculations were performed for Hg organometallic complexes in water, dimethyl sulfoxide, and chloroform. The spin-orbit coupling, a relativistic effect, is a key factor for predicting δ(Hg) and 1J(Hg-C) accurately, in conjunction with a dynamic treatment of the systems. Good agreement between the theoretical and experimental results is reached by adopting implicit (based on a continuum model) and explicit (solvent molecules treated quantum mechanically) solvation models. Broader trends appearing in the experimental data available in the literature are reproduced by the calculations, and therefore, quantum chemistry is able to assist in the assignment and interpretation of 199Hg NMR data. Less pronounced trends, such as changes in the 199Hg chemical shift in different systems with the same atom types bound to Hg, are too weak to be predicted reliably by the current state-of-the-art theoretical methods based on AIMD sampling and relativistic DFT with hybrid functionals for NMR calculations.

A DFT study on photodissociation of the palladium tetrachloride anion PdCl42− in solution

Photodissociation of PdCl42− in water is investigated using the time-dependent DFT method with the polarized continuum model. First, to find out the exchange–correlation functional that can well describe the present system we make following two comparisons using many functionals: (i) the DFT absorption spectrum vs. experiment; (ii) the DFT one-dimensional (1D) potential energy curves (PECs) vs. the previous MRCI 1D PECs. As a result of comparison, it is found that the range-separated hybrid functionals well reproduce the observed spectrum and, among these, the M11 functional may be the most reasonable choice. Next, using M11, we calculate ground and excited-state 1D PECs for single Cl loss in water and the ground-state two-dimensional (2D) potential energy surface (PES) for double Cl loss in water. It is revealed in the 1D PECs in water that electronic excitations to the PECs that have large transition dipole moments with the ground electronic state may easily lead to the PdCl3− + Cl− dissociation channel, which is similar to single Cl loss in the gas phase. It is also found that the ground-state 2D PES in water correlates with the PdCl2 + 2Cl− dissociation channel with no transition state on the way to it. This is different from the prediction in the gas phase; the ground-state 2D PES correlates with the PdCl2 + Cl2− dissociation channel with a transition state on the minimum energy path leading to it.

Promoting electrocatalytic activity of single-phase semiconductor Y3Fe5O12 (YIG) dual-functional cathode/electrolyte for ceramic fuel cells

To construct hybrid functions of electrolyte and cathode materials has improved fuel cell power output and structural stability at low operational temperatures as a novel approach in ceramics fuel cells (CFC) research. We found that semiconductor material Y3Fe5O12 possesses dual-functions of electrolyte and cathode, which has been systematically studied on the ion transport, oxygen kinetics, redox reaction, and electrochemical performance. The sintering temperature strongly affects the YIG, resulting in weak metal oxygen (M-O) bonds, more active sites and enrich oxygen vacancies, which provide both an active surfaces towards oxygen reduction and high oxide ion conductivity. The hybrid functional YIG exhibits good peak power output of 423 mW/cm2 (as an electrolyte) and 531 mW/cm2 (as an cathode) at 540 °C, respectively. Moreover, YIG displays outstanding oxygen reduction reaction (ORR) and catalytic activity during the electrochemical measurements. This work may offer a novel strategy for developing hybrid functions for ceramic fuel cells to enable low operational temperatures.

TD-DFT analysis of excitation of the closed-form spironaphthopyran in methanol solution: The contribution of vibronic transitions and intermolecular hydrogen bonds

A set of forty hybrid functionals, combined with a 6–31++G(d,p) basis set and an IEFPCM solvent description, was applied to calculate the S0→S1 excitation of the spironaphthopyran (SNP) photochromic compound. None of the hybrid functionals resulted in Cspiro-O bond cleavage upon optimization of the S1 excited state of the SNP, in contrast to the spiropyran BIPS which we studied earlier. At the same time, the deformations of the SNP during excitation turned out to be very small (in contrast to those for BIPS) which made it possible to calculate its vibronic absorption spectrum. The best agreement between the calculated and experimental absorption bands was demonstrated by the SOGGA11X functional. The main result of this work is confirmation of the vibration participation in the breaking of the Cspiro-O bond in the S1 excited state of the SNP. This assumption based on experimental data was previously put forward by other authors. Strong H-bonds between the SNP and two methanol molecules enhance the photoinduced distortions of the latter and are slightly weakened upon excitation.

Paramagnetic Nuclear Magnetic Resonance Shifts for Triplet Systems and Beyond with Modern Relativistic Density Functional Methods.

An efficient framework for the calculation of paramagnetic NMR (pNMR) shifts within exact two-component (X2C) theory and (current-dependent) density functional theory (DFT) up to the class of local hybrid functionals (LHFs) is presented. Generally, pNMR shifts for systems with more than one unpaired electron depend on the orbital shielding contribution and a temperature-dependent term. The latter includes zero-field splitting (ZFS), hyperfine coupling (HFC), and the g-tensor. For consistency, we calculate these three tensors at the same level of theory, i.e., using scalar-relativistic X2C augmented with spin-orbit perturbation theory. Results for pNMR chemical shifts of transition-metal complexes reveal that this X2C-DFT framework can yield good results for both the shifts and the individual tensor contributions of metallocenes and related systems, especially if the HFC constant is large. For small HFC constants, the relative error is often large, and sometimes the sign may be off. 4d and 5d complexes with more complicated structures demonstrate the limitations of a fully DFT-based approach. Additionally, a Co-based complex with a very large ZFS and pronounced multireference character is not well described. Here, a hybrid DFT-multireference framework is necessary for accurate results. Our results show that X2C is sufficient to describe relativistic effects and computationally cheaper than a fully relativistic approach. Thus, it allows use of large basis sets for converged HFCs. Overall, current-dependent meta-generalized gradient approximations and LHFs show some potential; however, the currently available functionals leave a lot to be desired, and the predictive power is limited.

The interplay of density functional selection and crystal structure for accurate NMR chemical shift predictions.

Ab initio chemical shift prediction plays a central role in nuclear magnetic resonance (NMR) crystallography, and the accuracy with which chemical shifts can be predicted relative to experiment impacts the confidence with which structures can be assigned. For organic crystals, periodic density functional theory calculations with the gauge-including projector augmented wave (GIPAW) approximation and the PBE functional are widely used at present. Many previous studies have examined how using more advanced density functionals can increase the accuracy of predicted chemical shifts relative to experiment, but nearly all of those studies employed crystal structures that were optimized with generalized-gradient approximation (GGA) functionals. Here, we investigate how the accuracy of the predicted chemical shifts in organic crystals is affected by replacing GGA-level PBE-D3(BJ) crystal geometries with more accurate hybrid functional PBE0-D3(BJ) ones. Based on benchmark data sets containing 132 13C and 35 15N chemical shifts, plus case studies on testosterone, acetaminophen, and phenobarbital, we find that switching from GGA-level geometries and chemical shifts to hybrid-functional ones reduces 13C and 15N chemical shift errors by ∼40-60% versus experiment. However, most of the improvement stems from the use of the hybrid functional for the chemical shift calculations, rather than from the refined geometries. In addition, even with the improved geometries, we find that double-hybrid functionals still do not systematically increase chemical shift agreement with experiment beyond what hybrid functionals provide. In the end, these results suggest that the combination of GGA-level crystal structures and hybrid-functional chemical shifts represents a particularly cost-effective combination for NMR crystallography in organic systems.

Mechanisms of photoisomerization of the prenylated flavin mononucleotide cofactor: a theoretical study.

The enzymatic decarboxylation of α,β-unsaturated acids using the ferulic acid decarboxylase (Fdc1) enzyme and prenylated flavin mononucleotide (prFMN) cofactor is a potential, environmentally friendly reaction for the biosynthesis of styrene and its derivatives. However, experiments showed that the enzyme activity of Fdc1 depends on the ring structure of prFMN, namely, the iminium and ketimine forms, and the loss of enzyme activity results from prFMNim → prFMNket photoisomerization. To obtain insight into this photochemical process and to improve the enzyme efficiency of Fdc1, two proposed photoisomerization mechanisms with different proton sources for the acid-base reaction were studied herein using theoretical methods. The potential energy surfaces calculated using the density functional theory method with the Becke, 3-parameter, and Lee-Yang-Parr hybrid functionals and DZP basis set (DFT/B3LYP/DZP) and TD-DFT/B3LYP/DZP methods confirmed that the light-dependent reaction occurs in the rate-determining proton transfer process and that the mechanism involving intermolecular proton transfer between prFMNim and Glu282 (external base) is energetically more favorable than that involving intramolecular proton transfer in prFMNim (internal base). The thermodynamic results obtained from the transition state theory method suggested that the exothermic relaxation energy in the photo-to-thermal process can promote the spontaneous formation of a high-energy-barrier transition state, and an effective enzymatic decarboxylation could be achieved by slowing down the formation of the undesirable thermodynamically favorable product (prFMNket). Because the rate constant for formation of the high-energy-barrier transition state varies exponentially over the temperature range of 273-298 K, and experimental results have shown that incubating Fdc1 on ice results in a complete loss of enzyme activity, it is recommended to perform the decarboxylation reaction at 285 K to strike a balance between minimizing enzyme stability loss at 273 K and mitigating the effects of UV irradiation. The computational strategy and fundamental insights obtained in this study could serve as guidelines for future theoretical and experimental investigations on the same and similar photochemical systems.

Adsorption and vibrational spectroscopy of CO on the surface of MgO from periodic local coupled-cluster theory.

The adsorption of CO on the surface of MgO has long been a model problem in surface chemistry. Here, we report periodic Gaussian-based calculations for this problem using second-order perturbation theory (MP2) and coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)], with the latter two performed using a recently developed extension of the local natural orbital approximation to problems with periodic boundary conditions. The low cost of periodic local correlation calculations allows us to calculate the full CCSD(T) binding curve of CO approaching the surface of MgO (and thus the adsorption energy) and the two-dimensional potential energy surface (PES) as a function of the distance from the surface and the CO stretching coordinate. From the PES, we obtain the fundamental vibrational frequency of CO on MgO, whose shift from the gas phase value is a common experimental probe of surface adsorption. We find that CCSD(T) correctly predicts a positive frequency shift upon adsorption of +14.7 cm-1, in excellent agreement with the experimental shift of +14.3 cm-1. We use our CCSD(T) results to assess the accuracy of MP2, CCSD, and several density functional theory (DFT) approximations, including exchange correlation functionals and dispersion corrections. We find that MP2 and CCSD yield reasonable binding energies and frequency shifts, whereas many DFT calculations overestimate the magnitude of the adsorption energy by 5-15 kJ mol-1 and predict a negative frequency shift of about -20 cm-1, which we attribute to self-interaction-induced delocalization errors that are mildly ameliorated with hybrid functionals. Our findings highlight the accuracy and computational efficiency of the periodic local correlation for the simulation of surface chemistry with accurate wavefunction methods.

Calculation of electric field gradients with the exact two-component (X2C) quasi-relativistic method and its local approximations.

When calculating electric field gradients (EFGs), relativistic and electron correlation effects are crucial for obtaining accurate results, and the commonly used density functional methods produce unsatisfactory results, especially for heavy elements and/or strongly correlated systems. In this work, a stand-alone program is presented, which enables calculation of EFGs from the molecular orbitals supplied by an external high accuracy quantum chemical calculation and includes relativistic effects through the exact two-component (X2C) formalism and efficient local approximations to it. Application to BiN and BiP molecules shows that a high precision can be achieved in the calculation of nuclear quadrupole coupling constants of 209Bi by combining advanced ab initio methods with the X2C approach. For seventeen iron compounds, the Mössbauer nuclear quadrupole splittings (NQS) of 57Fe calculated using a double-hybrid functional method are in very good agreement with the experimental values. It is shown that, for strongly correlated molecules, the double-hybrid functionals are much more accurate than the commonly used hybrid functionals. The computer program developed in this study furnishes a useful utility for obtaining EFGs and related nuclear properties with high accuracy.

Predicting spin states of iron porphyrins with DFT methods including crystal packing effects and thermodynamic corrections.

Accurate computational treatment of spin states for transition metal complexes, exemplified by iron porphyrins, lies at the heart of quantum bioinorganic chemistry, but at the same time represents a great challenge for approximate density functional theory (DFT) methods, which are predominantly used. Here, the accuracy of DFT methods for spin-state splittings in iron porphyrin is assessed by probing the ability to correctly predict the ground states for six FeIII or FeII complexes experimentally characterized in solid state. For each case, molecular and periodic DFT calculations are employed to quantify the effect of porphyrin side substituents and the crystal packing effect (CPE) on the spin-state splitting. It is proposed to partition the total CPE into additive components, the direct and structural one, the importance of which is shown to significantly vary from case to case. By knowing the substituent effect, the CPE, and the Gibbs free energy thermodynamic correction from calculations, one can employ the experimental ground-state information in order to derive a quantitative constraint on the electronic energy difference for a simplified (porphin) model of the experimentally characterized metalloporphyrin. The constraints derived in such a way-in the form of single or double inequalities-are used to assess the accuracy of dispersion-corrected DFT methods for 6 spin-state splittings of [FeIII(P)(2-MeIm)2]+, [FeIII(P)(2-MeIm)]+, [FeII(P)(THF)2] and [FeII(P)] models (where P is porphin, 2-MeIm is 2-methylimidazole, THF is tetrahydrofuran). These data constitute the new benchmark set of spin states for crystalline iron porphyrins (SSCIP6). The highest accuracy is obtained in the case of double-hybrid functionals (B2PLYP-D3, DSD-PBEB95-D3), whereas hybrid functionals, especially those with reduced admixture of the exact exchange (B3LYP*-D3, TPSSh-D3), are found to considerably overstabilize the intermediate spin state, leading to incorrect ground-state prediction in FeIII porphyrins. The present approach, which can be generalized to other transition metal complexes, is not only useful in method benchmarking, but also sheds light on the interpretations of experimental data for metalloporphyrins, which are important models to understand the electronic properties of heme proteins.

A review of first-principles calculation methods for defects in semiconductors

Doping and defect control in semiconductors are essential prerequisites for their practical applications. First-principles calculations of defects based on density functional theory offer crucial guidance for doping and defect control. In this paper, the developments in the theoretical methods of first-principles semiconductor defect calculations are introduced. Firstly, we introduce the method of calculating the defect formation energy and finite-size errors to the formation energy caused by the supercell method. Then, we present corresponding image charge correction schemes, which include the widely used post-hoc corrections (such as Makov-Payne, Lany-Zunger, Freysoldt-Neugebauer-van de Walle schemes), the recently developed self-consistent potential correction which performs the image charge correction in the self-consistent loop for solving Kohn-Sham equations, and the self-consistent charge correction scheme which does not require an input of macroscopic dielectric constants. Further, we extend our discussion to charged defect calculations in low-dimensional semiconductors, elucidate the issue of charged defect formation energy divergence with the increase of vacuum thickness within the jellium model and introduce our theoretical model which solves this energy divergence issue by placing the ionized electrons or holes in the realistic host band-edge states instead of the virtual jellium state. Furthermore, we provide a brief overview of defect calculation correction methods due to the DFT band gap error, including the scissors operator, LDA+<i>U</i> and hybrid functionals. Finally, in order to describe the calculation of defect formation energy under illumination, we present our self-consistent two-Fermi-reservoir model, which can well predict the defect concentration and carrier concentration in the Mg doped GaN system under illumination. This work summarizes the recent developments regarding first-principles calculations of defects in semiconducting materials and low-dimensional semiconductors, under whether equilibrium conditions or non-equilibrium conditions, thus promoting further developments of doping and defect control within semiconductors.

Tailoring intersystem crossing in phosphorus corroles through axial chalcogenation: a detailed theoretical study.

Intersystem crossing (ISC) of visible-light absorbing metal-free corrole macrocycles can be greatly tuned by means of suitable chemical functionalization. Axially chalcogenated phosphorus corrole derivatives (XPCs; X = O, S, Se) are expected to show large spin-orbit coupling (SOC) via the heavy-atom effect and therefore a much improved ISC. Excited-state deactivation of XPCs including PC is studied using time-dependent optimally tuned range-separated hybrid functionals combined with a polarizable continuum model with toluene as a dielectric medium to account for polar solvent effects. PC and all XPCs are dynamically stable and also show favourable thermodynamic formation feasibility as confirmed by Gibbs free energy analysis. In spite of the relatively smaller contribution of P and X to the frontier molecular orbitals compared to the tetrapyrrolic ring, SOC is considerably improved due to the heavy-atom effect. While PC shows a one-order larger ISC rate of ∼107 s-1 than fluorescence, competitive fluorescence and ISC rates of ∼107 s-1 are found for OPC. In contrast, both SPC and SePC exhibit significantly larger ISC rates of ∼109 s-1 and ∼1013 s-1, respectively, with much smaller fluorescence rates of ∼107 s-1. Importantly, the first report of anti-Kasha's emission in metal-free corroles is predicted for OPC with a radiative rate of ∼109 s-1. Furthermore, calculated phosphorescence and ISC rates from the near-degenerate lowest excited triplets to the ground-state suggest millisecond to microsecond triplet lifetimes, signalling towards long-lived excited triplet formation. Overall, all three XPCs including PC could act as triplet photosensitizers and especially both SPC and SePC are predicted to be the highly efficient ones.

Quantum-to-classical modeling of monolayer Ge2Se2 and its application in photovoltaic devices.

Since the discovery of graphene in 2004, the unique properties of two-dimensional materials have sparked intense research interest regarding their use as alternative materials in various photonic applications. Transition metal dichalcogenide monolayers have been proposed as transport layers in photovoltaic cells, but the promising characteristics of group IV-VI dichalcogenides are yet to be thoroughly investigated. This manuscript reports on monolayer Ge2Se2 (a group IV-VI dichalcogenide), its optoelectronic behavior, and its potential application in photovoltaics. When employed as a hole transport layer, the material fosters an astonishing device performance. We use ab initio modeling for the material prediction, while classical drift-diffusion drives the device simulations. Hybrid functionals calculate electronic and optical properties to maintain high accuracy. The structural stability has been verified using phonon spectra. The E-k dispersion reveals the investigated material's key electronic properties. The calculations reveal a direct bandgap of 1.12 eV for monolayer Ge2Se2. We further extract critical optical parameters using the Kubo-Greenwood formalism and Kramers-Kronig relations. A significantly large absorption coefficient and a high dielectric constant inspired the design of a monolayer Ge2Se2-based solar cell, exhibiting a high open circuit voltage of V oc = 1.11 V, a fill factor of 87.66%, and more than 28% power conversion efficiency at room temperature. Our findings advocate monolayer Ge2Se2 for various optoelectronic devices, including next-generation solar cells. The hybrid quantum-to-macroscopic methodology presented here applies to broader classes of 2D and 3D materials and structures, showing a path to the computational design of future photovoltaic materials.

Metal-organic framework-based high-performance column chip for micro gas chromatography: hybrid function for simultaneous preconcentration and separation of volatile organic compounds.

Numerous attempts have been made to replace commercial bulky gas chromatography (GC) systems with compact GC systems for monitoring volatile organic compounds in indoor air. However, recently developed compact GC systems are still too bulky in terms of user convenience, portability, and on-site analysis. Hence, an advanced miniaturization of compact GC systems is needed. Importantly, the small and high-performance gas pretreatment chip devices should be developed for compact GC systems. This paper reports the development of a metal-organic framework (MOF)-coated hybrid micro gas chromatography column chip (hybrid GC chip) capable of preconcentration and separation on harmful volatile organic compounds at low-concentration in one single chip device. The hybrid GC chip, 2 cm × 2 cm in size, was fabricated using a microelectromechanical systems process. The synthesized MOF-5 particles were coated on the inner wall of the fabricated hybrid GC chip and acted as an adsorbent and a stationary phase. The developed hybrid GC chip afforded high preconcentration factors with 1033-1237 and high separation resolutions with 3.8-13.3. The developed column showed good performance as a gas preconcentrator and separation column, and is the first device to perform both roles in one single chip device.

Phase transitions and spectral shifts: a quantum mechanical exploration of vibrational frequency in magnesium ferrite.

Spinel ferrites represent an integral subset of magnetic materials, with their inherent properties largely influenced by cation occupancy and spin interaction. In this study, we present an in-depth theoretical exploration of the phase transition landscape of pure magnesium-ferrite, deploying hybrid functionals and a local Gaussian basis set to scrutinize the relaxed lattice structure, relative energy, magnetic properties, electronic characteristics, and vibrational frequencies. Our investigation reveals that the ground state of magnesium-ferrite is an open-shell system with an inverse structure. This is characterized by the complete occupancy of octahedral sites by magnesium atoms, with Iron atoms dispersed between both tetrahedral and octahedral sites. We found a relative energy difference of 0.41 eV between the antiferromagnetic ground configuration and the ferro arrangement within the inverse structure. Furthermore, our research also delved into the impact of spin rearrangement and inversion (X = 0.0, 0.5 and 1) on Raman and infrared spectra. Notably, the lattice distortion from the cubic form, apparent in the optimized structure, resonates in the IR and Raman spectra, resulting in significant splitting. The frequencies calculated in this study align well with experimental values, suggesting that the literature's current assignments warrant reevaluation in light of this new data. The results presented herein can be instrumental in detecting the phase of Mg ferrites from experimental spectra, thereby paving the way for a more profound comprehension of their properties and possible uses.

Benchmarking DFT-based excited-state methods for intermolecular charge-transfer excitations.

Intermolecular charge-transfer is a highly important process in biology and energy-conversion applications where generated charges need to be transported over several moieties. However, its theoretical description is challenging since the high accuracy required to describe these excited states must be accessible for calculations on large molecular systems. In this benchmark study, we identify reliable low-scaling computational methods for this task. Our reference results were obtained from highly accurate wavefunction calculations that restrict the size of the benchmark systems. However, the density-functional theory based methods that we identify as accurate can be applied to much larger systems. Since targeting charge-transfer states requires the unambiguous classification of an excited state, we first analyze several charge-transfer descriptors for their reliability concerning intermolecular charge-transfer and single out the charge-transfer distance calculated based on the variation of electron density upon excitation (DCT) as an optimal choice for our purposes. In general, best results are obtained for orbital-optimized methods and among those, the maximum overlap method proved to be the most numerically stable variant when using the initial MOs as reference orbitals. Favorable error cancellation with optimally-tuned range-separated hybrid functionals and a rather small basis set can provide an economical yet reasonable wavefunction when using time-dependent density functional theory, which provides relevant information about the excited-state character to be used in the orbital-optimized methods. The qualitative agreement makes these fast calculations attractive for high-throughput screening applications.

Nutrient removal and yield of different maize hybrids

Objective: to determine the macro and micronutrient removal values and potential yield of different hybrids, and, also to determine the relationship between grain nutrient removal and grain yield. Design/methodology/approach: to assess correlations and determine the association degree between the nutrient removal values and grain yield. Results: the total nutrient removal values were in N> K> Ca> Mg> P, and Mn> Fe> Zn> B> Cu order, which are higher values when compared to another research. Also, these provide the mineral content in grains, which is a nutritional quality-related parameter. Limitations on study/implications: increasing the number of hybrids, different fertilization rates, different soil conditions, and crop management practices should be evaluated to assess whether these influence/inhibit the final nutrient concentration and total removal in grain. Findings/conclusions: The total grain nutrient removal values varied as a function of hybrids, yield goal, and nutrient concentration in tissues. These values allow the adjustment of current fertilization rates. The same hybrids under different management practices (fertilization dose), or soil types, substantially influence the grain nutrient concentration and therefore total nutrient removal.

DFT investigation on the structural and vibrational behaviours of the non-protein amino acids in hybrid explicit/continuum solvent: a case of the zwitterions γ-aminobutyric and α - aminoisobutyric acids.

The influence of hybrid solvation models on the molecular structures and vibrational characteristics of g-aminobutyric acid (GABA) and a-aminoisobutyric acid (AIB) zwitterions was assessed by employing a variety of Density Functional Theory (DFT). The quantum chemical methods included the B3LYP and B3PW91 hybrid functionals and the 6‑311++G(d,p) basis set. The most stable conformation derived from the potential energy surface (PES) scans using the B3LYP/6-311++G(d,p) model chemistry for each studied molecule was predicted within a continuum environment represented by the COSMO and SMD solvation models. The stable structures were subsequently immersed in explicit/COSMO and explicit/SMD hybrid solvation models, where 10 and 8 water molecules were explicitly positioned around the functional groups of the GABA and AIB zwitterions, respectively. The number of water molecules chosen was sufficient to prevent proton transfer among the carboxylate group (COO-) and the ammonium group (NH3+) within each molecule under investigation. After optimizing the geometry of each hydrated complex, the normal vibrational modes were determined. The scaled theoretical frequencies obtained from the various model chemistries were then compared to available experimental data from infrared (IR) and Raman spectroscopy. In the case of GABA and AIB molecules, the comparisons revealed that the B3LYP/6-311++G(d,p) model chemistry yielded wavenumber values that closely matched the experimental IR and Raman data, particularly when the explicit/SMD solvent was employed. The computed results indicate deviations of less than 4% when compared to the experimental data for the two molecules under investigation.

On the choice of exchange–correlation functional for the correct band gap calculation in titanium oxide photocatalysts

In this work the band gap for a set of titanium oxide photocatalysts: TiO2 (rutile, brookite, anatase), SrTiO3, BaTiO3, CaTiO3, MgTiO3, KLaTiO4, K2La2Ti3O10, was theoretically investigated within density functional theory. For the PBE-based hybrid functional the band gap as a function of the fraction of the Hartree-Fock exchange energy (α) was calculated. A correlation between the Goldschmidt tolerance factor and the value of α parameter, that provides an agreement with the experimental band gap, was determined. A wide range of standard hybrid functionals was tested. It was found that the B1WC and WC1LYP hybrid functionals (with α = 0.16) provide the best agreement with experimental data for the set of studied compounds. However, for a more accurate description of the band structure of titanium oxide photocatalysts the α parameter should be chosen accounting for the structural distortion of TiO6 octahedra and their mutual arrangement. The results obtained can be used in predicting the properties of modified structures based on titanium oxides.