Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

The hardness profile as a tool to detect spurious stationary points in the potential energy surface

The energy and hardness profile for a series of inter and intramolecular conformational changes at several levels of calculation were computed. The hardness profiles were found to be calculated as the difference between the vertical ionization potential and electron affinity. The hardness profile shows the correct number of stationary points independently of the basis set and methodology used. It was found that the hardness profiles can be used to check the reliability of the energy profiles for those chemical system

Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential VA on a one-atom density ρ A or ρ B. To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale

Distributed all-atom simulations of intrinsically unstructured proteins

The Computational Biophysics Group at the Universitat Pompeu Fabra (GRIB-UPF) hosts two unique computational resources dedicated to the execution of large scale molecular dynamics (MD) simulations: (a) the ACMD molecular-dynamics software, used on standard personal computers with graphical processing units (GPUs); and (b) the GPUGRID. net computing network, supported by users distributed worldwide that volunteer GPUs for biomedical research. We leveraged these resources and developed studies, protocols and open-source software to elucidate energetics and pathways of a number of biomolecular systems, with a special focus on flexible proteins with many degrees of freedom. First, we characterized ion permeation through the bactericidal model protein Gramicidin A conducting one of the largest studies to date with the steered MD biasing methodology. Next, we addressed an open problem in structural biology, the determination of drug-protein association kinetics; we reconstructed the binding free energy, association, and dissaciociation rates of a drug like model system through a spatial decomposition and a Makov-chain analysis. The work was published in the Proceedings of the National Academy of Sciences and become one of the few landmark papers elucidating a ligand-binding pathway. Furthermore, we investigated the unstructured Kinase Inducible Domain (KID), a 28-peptide central to signalling and transcriptional response; the kinetics of this challenging system was modelled with a Markovian approach in collaboration with Frank Noe’s group at the Freie University of Berlin. The impact of the funding includes three peer-reviewed publication on high-impact journals; three more papers under review; four MD analysis components, released as open-source software; MD protocols; didactic material, and code for the hosting group.

Non-invasive monitoring of diaphragmatic timing by means of surface contact sensors: an experimental study in dogs

Background: Non-invasive monitoring of respiratory muscle function is an area of increasing research interest, resulting in the appearance of new monitoring devices, one of these being piezoelectric contact sensors. The present study was designed to test whether the use of piezoelectric contact (non-invasive) sensors could be useful in respiratory monitoring, in particular in measuring the timing of diaphragmatic contraction.Methods: Experiments were performed in an animal model: three pentobarbital anesthetized mongrel dogs. The motion of the thoracic cage was acquired by means of a piezoelectric contact sensor placed on the costal wall. This signal is compared with direct measurements of the diaphragmatic muscle length, made by sonomicrometry. Furthermore, to assess the diaphragmatic function other respiratory signals were acquired: respiratory airflow and transdiaphragmatic pressure. Diaphragm contraction time was estimated with these four signals. Using diaphragm length signal as reference, contraction times estimated with the other three signals were compared with the contraction time estimated with diaphragm length signal.Results: The contraction time estimated with the TM signal tends to give a reading 0.06 seconds lower than the measure made with the DL signal (-0.21 and 0.00 for FL and DP signals, respectively), with a standard deviation of 0.05 seconds (0.08 and 0.06 for FL and DP signals, respectively). Correlation coefficients indicated a close link between time contraction estimated with TM signal and contraction time estimated with DL signal (a Pearson correlation coefficient of 0.98, a reliability coefficient of 0.95, a slope of 1.01 and a Spearman's rank-order coefficient of 0.98). In general, correlation coefficients and mean and standard deviation of the difference were better in the inspiratory load respiratory test than in spontaneous ventilation tests.Conclusion: The technique presented in this work provides a non-invasive method to assess the timing of diaphragmatic contraction in canines, using a piezoelectric contact sensor placed on the costal wall.

Determination of sea surface salinity and wind speed by L-band microwave radiometry from a fixed platform

The European Space Agency Soil Moisture andOcean Salinity (SMOS) mission aims at obtaining global maps ofsoil moisture and sea surface salinity from space for large-scale andclimatic studies. It uses an L-band (1400–1427 MHz) MicrowaveInterferometric Radiometer by Aperture Synthesis to measurebrightness temperature of the earth’s surface at horizontal andvertical polarizations ( h and v). These two parameters will beused together to retrieve the geophysical parameters. The retrievalof salinity is a complex process that requires the knowledge ofother environmental information and an accurate processing ofthe radiometer measurements. Here, we present recent resultsobtained from several studies and field experiments that were partof the SMOS mission, and highlight the issues still to be solved.

Surface topography and mixed-pixel effects on the simulated L-band brightness temperatures.

The impact of topography and mixed pixels on L-band radiometric observations over land needs to be quantified to improve the accuracy of soil moisture retrievals. For this purpose, a series of simulations has been performed with an improved version of the soil moisture and ocean salinity (SMOS) end-to-end performance simulator (SEPS). The brightness temperature generator of SEPS has been modified to include a 100-m-resolution land cover map and a 30-m-resolution digital elevation map of Catalonia (northeast of Spain). This high-resolution generator allows the assessment of the errors in soil moisture retrieval algorithms due to limited spatial resolution and provides a basis for the development of pixel disaggregation techniques. Variation of the local incidence angle, shadowing, and atmospheric effects (up- and downwelling radiation) due to surface topography has been analyzed. Results are compared to brightness temperatures that are computed under the assumption of an ellipsoidal Earth.

Cholesterol and Fatty Acids Regulate Dynamic Caveolin Trafficking through the Golgi Complex and between the Cell Surface and Lipid Bodies

Caveolins are a crucial component of plasma membrane (PM) caveolae but have also been localized to intracellular compartments, including the Golgi complex and lipid bodies. Mutant caveolins associated with human disease show aberrant trafficking to the PM and Golgi accumulation. We now show that the Golgi pool of mainly newly synthesized protein is detergent-soluble and predominantly in a monomeric state, in contrast to the surface pool. Caveolin at the PM is not recognized by specific caveolin antibodies unless PM cholesterol is depleted. Exit from the Golgi complex of wild-type caveolin-1 or -3, but not vesicular stomatitis virus-G protein, is modulated by changing cellular cholesterol levels. In contrast, a muscular dystrophy-associated mutant of caveolin-3, Cav3P104L, showed increased accumulation in the Golgi complex upon cholesterol treatment. In addition, we demonstrate that in response to fatty acid treatment caveolin can follow a previously undescribed pathway from the PM to lipid bodies and can move from lipid bodies to the PM in response to removal of fatty acids. The results suggest that cholesterol is a rate-limiting component for caveolin trafficking. Changes in caveolin flux through the exocytic pathway can therefore be an indicator of cellular cholesterol and fatty acid levels.

Ab initio study of NOx compounds adsorption on SnO2 surface

An ab initio study of the adsorption processes on NOx compounds on (1 1 0) SnO2 surface is presented with the aim of providing theoretical hints for the development of improved NOx gas sensors. From first principles calculations (DFT¿GGA approximation), the most relevant NO and NO2 adsorption processes are analyzed by means of the estimation of their adsorption energies. The resulting values and the developed model are also corroborated with experimental desorption temperatures for NO and NO2, allowing us to explain the temperature-programmed desorption experiments. The interference of the SO2 poisoning agent on the studied processes is discussed and the adsorption site blocking consequences on sensing response are analyzed.

Surface states in template synthesized tin oxide nanoparticles

Tin-oxide nanoparticles with controlled narrow size distributions are synthesized while physically encapsulated inside silica mesoporous templates. By means of ultraviolet-visible spectroscopy, a redshift of the optical absorbance edge is observed. Photoluminescence measurements corroborate the existence of an optical transition at 3.2 eV. The associated band of states in the semiconductor gap is present even on template-synthesized nanopowders calcined at 800°C, which contrasts with the evolution of the gap states measured on materials obtained by other methods. The gap states are thus considered to be surface localized, disappearing with surface faceting or being hidden by the surface-to-bulk ratio decrease.

Gas collisions and pressure quenching of the photoluminescence of silicon nanopowder grown by plasma-enhanced chemical vapor deposition

The quenching of the photoluminescence of Si nanopowder grown by plasma-enhanced chemical vapor deposition due to pressure was measured for various gases ( H2, O2, N2, He, Ne, Ar, and Kr) and at different temperatures. The characteristic pressure, P0, of the general dependence I(P) = I0¿exp(¿P/P0) is gas and temperature dependent. However, when the number of gas collisions is taken as the variable instead of pressure, then the quenching is the same within a gas family (mono- or diatomic) and it is temperature independent. So it is concluded that the effect depends on the number of gas collisions irrespective of the nature of the gas or its temperature.

Influence of premetallization surface treatment on the formation of Schottky Au-nGaN contacts

The influence of premetallization surface preparation on the structural, chemical, and electrical properties of Au-nGaN interfaces has been investigated by x-ray photoemission spectroscopy (XPS), current-voltage measurement (I-V) and cross-section transmission electron microscopy (TEM). XPS analysis showed that the three GaN substrate treatments investigated i.e., ex situ hydrofluoric acid etch, in situ anneal in ultrahigh-vacuum (UHV), and in situ Ga reflux cleaning in UHV result in surfaces increasingly free of oxygen contamination. XPS and TEM characterization of Au-nGaN formed after the three premetallization surface treatments show that HF etching and UHV annealing produce abrupt, well-defined interfaces. Conversely, GaN substrate cleaning in a Ga flux results in Au/GaN intermixing. I-V characterization of Au¿nGaN contacts yields a Schottky barrier height of 1.25 eV with a very low-ideality factor and very good contact uniformity for the premetallization UHV anneal, while the Ga reflux cleaning results in a much lower barrier (0.85 eV), with poor ideality and uniformity. I-V and XPS results suggest a high density of acceptor states at the surface, which is further enhanced by UHV annealing. These results are discussed in the context of current models of Schottky barrier formation.