In Silico Ultrafast Nonlinear Spectroscopy Meets Experiments: The Case of Perylene Bisimide Dye.

Spectroscopy simulations are of paramount importance for the interpretation of experimental electronic spectra, the disentangling of overlapping spectral features, and the tracing of the microscopic origin of the observed signals. Linear and nonlinear simulations are based on the results drawn from electronic structure calculations that provide the necessary parameterization of the molecular systems probed by light. Here, we investigate the applicability of excited-state properties obtained from linear-response time-dependent density functional theory (TDDFT) in the description of nonlinear spectra by employing the pseudowavefunction approach and compare them with benchmarks from highly accurate RASSCF/RASPT2 calculations and with high temporal resolution experimental results. As a test case, we consider the prediction of femtosecond transient absorption and two-dimensional electronic spectroscopy of a perylene bisimide dye in solution. We find that experimental signals are well reproduced by both theoretical approaches, showing that the computationally cheaper TDDFT can be a suitable option for the simulation of nonlinear spectroscopy of molecular systems that are too large to be treated with higher-level RASSCF/RASPT2 methods.

Chemiresistive Detection of Gaseous Hydrocarbons and Interrogation of Charge Transport in Cu[Ni(2,3-pyrazinedithiolate)2] by Gas Adsorption.

The development of new chemiresistive materials for use in chemical sensors that operate near ambient conditions could potentially reduce the costs of implementation, encouraging their use in new areas. Conductive metal-organic frameworks represent one intriguing class of materials for further investigation in this area, given their vast structural diversity and the specificity of adsorbate interactions afforded by their crystallinity. Here, we re-examine the electronic conductivity of the desolvated and acetonitrile-solvated microporous framework Cu[Ni(pdt)2] (pdt2- = 2,3-pyrazinedithiolate), and find that the conductivity in the pristine material is 200-fold greater than in the solvated state, highlighting the sensitivity of sample conductivity to guest inclusion. Additionally, the desolvated material is demonstrated to selectively adsorb the gaseous hydrocarbons ethane, ethylene, acetylene, propane, propylene, and cis-2-butene at ambient temperature. Investigation of the effect of gas adsorption on conductivity using an in situ measurement cell reveals a chemiresistive response for each adsorbate, and the change in conductivity with adsorbate pressure closely follows an empirical model identical in form to the Langmuir-Freundlich equation. The relative sensitivity of the framework to each adsorbate is, surprisingly, not correlated with binding strength. Instead, the differences in chemiresistive response between adsorbates are found to correlate strongly with gas phase specific heat capacity of the adsorbate. Nanoconfinement effects, manifesting as a relative deviation from the expected chemiresistive response, may influence charge transport in the case of the largest adsorbate considered, cis-2-butene. Time-resolved conductance and adsorption measurements additionally show that the chemiresistive response of the sensor equilibrates on a shorter time scale than gas adsorption, suggesting that interparticle contacts limit conduction through the bulk material and that conductivity at the crystallite surfaces is most responsive to gas adsorption.

A Modular Implementation for the Simulation of 1D and 2D Solid-State NMR Spectra of Quadrupolar Nuclei in the Virtual Multifrequency Spectrometer-Draw Graphical Interface.

We present the implementation of the solid state (SoS)NMR module for the simulation of several 1D and 2D NMR spectra of all the elements in the periodic table in the virtual multifrequency spectrometer (VMS). This module is fully integrated with the graphical user interface of VMS (VMS-Draw) [Licari et al., J. Comput. Chem. 36, 2015, 321-334], a freeware tool which allows a user-friendly handling of structures and analyses of advanced spectroscopical properties of chemical compounds-from model systems to real-world applications. Besides the numerous modules already available in VMS for the study of electronic, optical, vibrational, vibronic, and EPR properties, here the simulation of NMR spectra is presented with a particular emphasis on those techniques usually employed to investigate solid state systems. The SoSNMR module benefits from its ability to work under both periodic and nonperiodic conditions, such that small molecules/molecular clusters can be treated, as well as extended three-dimensional systems enforcing (or not) translational periodicity. These features allow VMS to simulate spectra resulting from NMR calculations by some popular quantum chemistry codes, namely Gaussian09/16, Castep, and Quantum Espresso. The effectiveness of the SoSNMR module of VMS is examined throughout the manuscript, and applied to simulate 1D static, MAS, and VAS NMR spectra as well as 2D correlation (90°, MAS) and MQMAS spectra of active NMR nuclei embedded in different amorphous and crystalline systems of actual interest in chemistry and material science. Finally, the program is able to simulate the spectra of both the total ensemble of spin-active nuclei present in the system and of subensembles differentiated depending on the chemical environment of the first and second coordination sphere in a very general way applicable to any kind of systems.

Rolling resistance contribution to a road pavement life cycle carbon footprint analysis.

PURPOSE: Although the impact of road pavement surface condition on rolling resistance has been included in the life cycle assessment (LCA) framework of several studies in the last years, there is still a high level of uncertainty concerning the methodological assumptions and the parameters that can affect the results. In order to adopt pavement carbon footprint/LCA as a decision-making tool, it is necessary to explore the impact of the chosen methods and assumptions on the LCA results. METHODS: This paper provides a review of the main models describing the impact of the pavement surface properties on vehicle fuel consumption and analyses the influence of the methodological assumptions related to the rolling resistance on the LCA results. It compares the CO2 emissions, calculated with two different rolling resistance models existing in literature, and performs a sensitivity test on some specific input variables (pavement deterioration rate, traffic growth, and emission factors/fuel efficiency improvement). RESULTS AND DISCUSSION: The model used to calculate the impact of the pavement surface condition on fuel consumption significantly affects the LCA results. The pavement deterioration rate influences the calculation in both models, while traffic growth and fuel efficiency improvement have a limited impact on the vehicle CO2 emissions resulting from the pavement condition contribution to rolling resistance. CONCLUSIONS AND RECOMMENDATIONS: Existing models linking pavement condition to rolling resistance and hence vehicle emissions are not broadly applicable to the use phase of road pavement LCA and further research is necessary before a widely-used methodology can be defined. The methods of modelling and the methodological assumptions need to be transparent in the analysis of the impact of the pavement surface condition on fuel consumption, in order to be interpreted by decision makers and implemented in an LCA framework. This will be necessary before product category rules (PCR) for pavement LCA can be extended to include the use phase.

Modeling emission features of salicylidene aniline molecular crystals: A QM/QM' approach.

A new computational protocol relying on the use of electrostatic embedding, derived from QM/QM' ONIOM calculations, to simulate the effect of the crystalline environment on the emission spectra of molecular crystals is here applied to the ß-form of salicylidene aniline (SA). The first singlet excited states (S1 ) of the SA cis-keto and trans-keto conformers, surrounded by a cluster of other molecules representing the crystalline structure, were optimized by using a QM/QM' ONIOM approach with and without electronic embedding. The model system consisting of the central salicylidene aniline molecule was treated at the DFT level by using either the B3LYP, PBE0, or the CAM-B3LYP functional, whereas the real system was treated at the HF level. The CAM-B3LYP/HF level of theory provides emission energies in good agreement with experiment with differences of -20/-32 nm (cis-keto form) and -8/-14 nm (trans-keto form), respectively, whereas notably larger differences are obtained using global hybrids. Though such differences on the optical properties arise from the density functional choice, the contribution of the electronic embedding is rather independent of the functional used. This plays in favor of a more general applicability of the present protocol to other crystalline molecular systems.

Insights into structural and dynamical features of water at halloysite interfaces probed by DFT and classical molecular dynamics simulations.

Density functional theory calculations and classical molecular dynamics simulations have been used to investigate the structure and dynamics of water molecules on kaolinite surfaces and confined in the interlayer of a halloysite model of nanometric dimension. The first technique allowed us to accurately describe the structure of the tetrahedral-octahedral slab of kaolinite in vacuum and in interaction with water molecules and to assess the performance of two widely employed empirical force fields to model water/clay interfaces. Classical molecular dynamics simulations were used to study the hydrogen bond network structure and dynamics of water adsorbed on kaolinite surfaces and confined in the halloysite interlayer. The results are in nice agreement with the few experimental data available in the literature, showing a pronounced ordering and reduced mobility of water molecules at the hydrophilic octahedral surfaces of kaolinite and confined in the halloysite interlayer, with respect to water interacting with the hydrophobic tetrahedral surfaces and in the bulk. Finally, this investigation provides new atomistic insights into the structural and dynamical properties of water-clay interfaces, which are of fundamental importance for both natural processes and industrial applications.

Prediction uncertainty of density functional approximations for properties of crystals with cubic symmetry.

The performance of a method is generally measured by an assessment of the errors between the method's results and a set of reference data. The prediction uncertainty is a measure of the confidence that can be attached to a method's prediction. Its estimation is based on the random part of the errors not explained by reference data uncertainty, which implies an evaluation of the systematic component(s) of the errors. As the predictions of most density functional approximations (DFA) present systematic errors, the standard performance statistics, such as the mean of the absolute errors (MAE or MUE), cannot be directly used to infer prediction uncertainty. We investigate here an a posteriori calibration method to estimate the prediction uncertainty of DFAs for properties of solids. A linear model is shown to be adequate to address the systematic trend in the errors. The applicability of this approach to modest-size reference sets (28 systems) is evaluated for the prediction of band gaps, bulk moduli, and lattice constants with a wide panel of DFAs.

Unraveling the polymorphism of [(p-cymene)Ru(κN-INA)Cl2] through dispersion-corrected DFT and NMR GIPAW calculations.

The structural and (13)C/(1)H NMR parameters of the four crystal forms (1α, 1·H2O, 1ß, and 1γ) of the solid wheel-and-axle (WAA) metal-organic compound [(p-cymene)Ru(κN-INA)Cl2] have been studied by means of periodic DFT calculations. The quality of the results obtained strongly depends on a correct description of long-range interactions; thus, in the geometry refinement protocol used, the pure DFT functionals need to be coupled with a dispersion-correction term (B3LYP-D2, B3LYP-D*). The solid-state (13)C/(1)H NMR δ(iso) parameters and (13)C MAS NMR spectra, calculated by means of the PBE-GIPAW method, agree well with the experimental data for the four crystal forms (mean absolute deviations of the (13)C and (1)H δ(iso) data values lie in the ranges 1.3-2.9 and 0.3-1.0 ppm, respectively). In this context, some revisions in the experimental assignment of the (13)C/(1)H NMR δ(iso) parameters of the 1·H2O, 1ß, and 1γ crystal forms can be suggested. The mismatch in the assignment seems to be due to the rotation of the -COOH moiety, which occurs at the 1α-1·H2O transition and was not considered in the experiments. Finally, the results obtained suggest the presence of two COOH···Cl hydrogen bonds of comparable strength established by the two molecules in the asymmetric unit of the 1γ polymorph, in partial disagreement with previous findings.

Computational Protocol for Modeling Thermochromic Molecular Crystals: Salicylidene Aniline As a Case Study.

A computational protocol that combines periodic and QM/QM' calculations has been applied to investigate the structural (geometrical and electronic) and photophysical absorption properties of the salicylidene aniline (SA) thermochromic molecular crystal. The protocol consists of three different steps, namely (i) the description of the molecular crystal using a periodic approach taking into account dispersion interactions, (ii) the identification of reliable finite models (clusters), and (iii) the calculation of vertical transition energies including environmental effects through the use of an electronic embedding model (QM/QM' ONIOM approach). The encouraging results obtained in this work for the ß polymorph of SA, both in terms of accuracy and computational cost, open the way to the simulation and the prediction of the photophysical behavior of other molecular crystals, especially those much less well characterized experimentally.