Ligand dynamics in myoglobin: calculation of infrared spectra for photodissociated NO

We present molecular dynamics simulations of the photodissociated state of MbNO performed at 300 K using a fluctuating charge model for the nitric oxide (NO) ligand. After dissociation, NO is observed to remain mainly in the centre of the distal haem pocket, although some movement towards the primary docking site and the xenon-4 pocket can be seen. We calculate the NO infrared spectrum for the photodissociated ligand within the haem pocket and find a narrow peak in the range 1915-1922 cm(-1). The resulting blue shift of 1 to 8 cm(-1) compared to gas-phase NO is much smaller than the red shifts calculated and observed for carbon monoxide (CO) in Mb. A small splitting, due to NO in the xenon-4 pocket, is also observed. At lower temperatures, the spectra and conformational space explored by the ligand remain largely unchanged, but the electrostatic interactions with residue His64 become increasingly significant in determining the details of the ligand orientation within the distal haem pocket. The investigation of the effect of the L29F mutation reveals significant differences between the behaviour of NO and that of CO, and suggests a coupling between the ligand and the protein dynamics due to the different ligand dipole moments.

Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm−1 (compared to the experimental value of 10 cm−1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol−1 and agrees well with the experimentally estimated value of 0.29 kcal mol−1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.

A new interpretation of the bonding properties and UV–vis spectra of [M3(CO)12] clusters (M = Ru, Os): a TD-DFT study

DFT and TD-DFT calculations (ADF program) were performed in order to analyze the electronic structure of the [M-3(CO)(12)] clusters (M = Ru, Os) and interpret their electronic spectra. The highest occupied molecular orbitals are M-M bonding (sigma) involving different M-M bonds, both for Ru and Os. They participate in low-energy excitation processes and their depopulation should weaken M-M bonds in general. While the LUMO is M-NI and M-CO anti-bonding (sigma*), the next, higher-lying empty orbitals have a main contribution from CO (pi*) and either a small (Ru) or an almost negligible one (Os) from the metal atoms. The main difference between the two clusters comes from the different nature of these low-energy unoccupied orbitals that have a larger metal contribution in the case of ruthenium. The photochemical reactivity of the two clusters is reexamined and compared to earlier interpretations.

Infrared absorption spectra, radiative efficiencies, and global warming potentials of perfluorocarbons: Comparison between experiment and theory

Experimentally and theoretically determined infrared spectra are reported for a series of straight-chain perfluorocarbons: C2F6, C3F8, C4F10, C5F12, C6F14, and C8F18. Theoretical spectra were determined using both density functional (DFT) and ab initio methods. Radiative efficiencies (REs) were determined using the method of Pinnock et al. (1995) and combined with atmospheric lifetimes from the literature to determine global warming potentials (GWPs). Theoretically determined absorption cross sections were within 10% of experimentally determined values. Despite being much less computationally expensive, DFT calculations were generally found to perform better than ab initio methods. There is a strong wavenumber dependence of radiative forcing in the region of the fundamental C-F vibration, and small differences in wavelength between band positions determined by theory and experiment have a significant impact on the REs. We apply an empirical correction to the theoretical spectra and then test this correction on a number of branched chain and cyclic perfluoroalkanes. We then compute absorption cross sections, REs, and GWPs for an additional set of perfluoroalkenes.

Synthesis, X-ray crystal structures and properties of complex salts and sterically crowded heteroleptic complexes of group 10 metal ions with aromatic sulfonyl dithiocarbimates and triphenylphosphine ligand

Two new complex salts of the form (Bu4N)(2)[Ni(L)(2)] (1) and (Ph4P)(2)[Ni(L)(2)] (2) and four heteroleptic complexes cis-M(PPh3)(2)(L) [M = Ni(II) (3), Pd(II) (4), L = 4-CH3OC6H4SO2N=CS2] and cis-M(PPh3)(2)(L') [M = Pd(II) (5), Pt(II) (6), L' = C6H5SO2N=CS2] were prepared and characterized by elemental analyses, IR, H-1, C-13 and P-31 NMR and UV-Vis spectra, solution and solid phase conductivity measurements and X-ray crystallography. A minor product trans-Pd(PPh3)(2)(SH)(2), 4a was also obtained with the synthesis of 4. The NiS4 and MP2S2 core in the complex salts and heteroleptic complexes are in the distorted square-plane whereas in the trans complex, 4a the centrosymmetric PdS2P2 core is perforce square planar. X-ray crystallography revealed the proximity of the ortho phenyl proton of the PPh3 ligand to Pd(II) showing rare intramolecular C-H center dot center dot center dot Pd anagostic binding interactions in the palladium cis-5 and trans-4a complexes. The complex salts with sigma(rt) values similar to 10 (5) S cm (1) show semi-conductor behaviors. The palladium and platinum complexes show photoluminescence properties in solution at room temperature.

Origin of electronic absorption spectra of MLCT-excited and one-electron reduced 2,2′-bipyridine and 1,10-phenanthroline complexes

UV–Vis absorption spectra of one-electron reduction products and 3MLCT excited states of [ReICl(CO)3- (N,N)] (N,N = 2,20-bipyridine, bpy; 1,10-phenanthroline, phen) have been measured by low-temperature spectroelectrochemistry and UV–Vis transient absorption spectroscopy, respectively, and assigned by open-shell TD-DFT calculations. The characters of the electronic transitions are visualized and analyzed using electron density redistribution maps. It follows that reduced and excited states can be approximately formulated as [ReICl(CO)3(N,Nÿ)]ÿ and ⁄[ReIICl(CO)3(N,Nÿ)], respectively. UV–Vis spectra of the reduced complexes are dominated by IL transitions, plus weaker MLCT contributions. Excited-state spectra show an intense band in the UV region of 50% IL origin mixed with LMCT (bpy, 373 nm) or MLCT (phen, 307 nm) excitations. Because of the significant IL contribution, this spectral feature is akin to the principal IL band of the anions. In contrast, the excited-state visible spectral pattern arises from predominantly LMCT transitions, any resemblance with the reduced-state visible spectra being coincidental. The Re complexes studied herein are representatives of a broad class of metal a-diimines, for which similar spectroscopic behavior can be expected.

Structure and spectra of dichloro(hydroxy-methoxy-di(2-pyridyl)methane)copper(II)

The non-electrolyte dichloro(hydroxy-methoxy-di(2-pyridylmethane)copper(II), resulting from the reaction of di(2-pyridyl)ketone and copper(II) chloride in methanol solution, was isolated and characterized and its structure was determined by X-ray diffraction. The pyridyl nitrogens and the chloride anions virtually from a basal plane in which lies the copper atom, while the oxygen of the methoxy group is in an apical position at a distance of 2.497 (3)Å. The nitrogenous base adopts the boat conformation with the pyridyl rings forming a dihedral angle of 108.72 (14)°. The nearest interatomic copper distance of 3.940(3)Å precludes copper-copper interactions, while the proximity of copper to the out-of-plane chlorine atoms [3.109(3)Å] suggests weakly bound chloro-bridged dimers. Spectral changes indicate that protic molecules displace the methoxy group and water affords the corresponding 1,1-diol.

Influence of hydrogen bonding in the structure of mixed-ligand copper(II) chelates. Spectra and structure of [aquo((2,2'-dipyridylamine)-(3-chloro-2,4-pentanedionato)copper(II)] nitrate

The IR, the ligand field spectra and the crystal structure of the mixed-ligand compound [(aquo)2,2P1 , a = 8.718(5), b = 9.407(5), c = 13.484 (7) Å, = 94.17(4)°, = 105.12(5)°, = 119.75(5)°, Z = 2, R = 0.0332, R W = 0.0869).

Synthesis, characterisation, Mössbauer spectra, and structures of some trinuclear iron-tin clusters

Reactions of [Fe3(CO)12] with diaryltin species SnR2(R1= 2,4,6-triisopropylphenyl, R2= 2,6-diethylphenyl, R3= pentamethylphenyl) and with Sn[CH(PPh2)2]2 have been investigated. The tin reagents SnR2(R = R1 or R2) reacted under mild conditions to give in moderate yields the trinuclear species [Fe2(CO)8(µ-SnR12)]1 or [Fe2(CO)8(µ-SnR22)]2, as orange-red crystalline solids, which decompose in air on prolonged exposure. The compound [Fe2(CO)8(µ-SnR42)]3(R4= 2,4,6-triphenylphenyl) can be similarly obtained. Prolonged treatment of the carbonyl with the novel tin reagent SnR32, by contrast, afforded the known compound spiro-[(OC)8Fe2SnFe2(CO)8]4 for which data are briefly reported. Reactions with tin or lead reagents M[CH(PPh2)2]2(M = Sn or Pb) afforded [Fe2(CO)6(µ-CO)(µ-dppm)][dppm = 1,2-bis(diphenylphosphino)methane] rapidly and almost quantitatively. Full crystal and molecular structural data are reported for [Fe2(CO)8(µ-SnR12)] and [Fe2(CO)8(µ-SnR22)]. Mössbauer data are also presented for compounds 1–3, and interpreted in terms of the structural data for these and other systems.

The interplay of secondary Hg⋯S, Hg⋯N and Hg⋯π bonding interactions in supramolecular structures of phenylmercury(ii) dithiocarbamates

Three new phenylmercury(II) and one mercury(II) dithiocarbamate complexes viz. PhHg S2CN(PyCH2) Bz (1), PhHg S2CN(PyCH2)CH3 (2), PhHg S2CN(Bz)CH3 (3), and [Hg (NCS2(PyCH2)Bz)(2)] (4) (Py = pyridine; Bz = benzyl) have been synthesized and characterized by elemental analyses, IR, electronic absorption, H-1 and C-13 NMR spectroscopy. The crystal structures of 1, 2 and 3 showed a linear S-Hg-C core at the centre of the molecule, in which the metal atom is bound to the sulfur atom of the dithiocarbamate ligand and a carbon atom of the aromatic ring. In contrast the crystal structure of 4 showed a linear S-Hg-S core at the Hg(II) centre of the molecule. Weak intermolecular Hg center dot center dot center dot N (Py) interactions link molecules into a linear chain in the case of 1, whereas chains of dimers are formed in 2 through intermolecular Hg center dot center dot center dot N (Py) and Hg center dot center dot center dot S interactions. 3 forms a conventional face-to-edge dimeric structure through intermolecular Hg center dot center dot center dot S secondary bonding and 4 forms a linear chain of dimers through face-to-face Hg center dot center dot center dot S secondary bonding. In order to elucidate the nature of these secondary bonding interactions and the electronic absorption spectra of the complexes, ab initio quantum chemical calculations at the MP2 level and density functional theory calculations were carried out for 1-3. Complexes 1 and 2 exhibited photoluminescent properties in the solid state as well as in the solution phase. Studies indicate that Hg center dot center dot center dot S interactions decrease and Hg center dot center dot center dot N interactions increase the chances of photoluminescence in the solid phase

The troposphere-to-stratosphere transition in kinetic energy spectra and nonlinear spectral fluxes as seen in ECMWF analyses

Global horizontal wavenumber kinetic energy spectra and spectral fluxes of rotational kinetic energy and enstrophy are computed for a range of vertical levels using a T799 ECMWF operational analysis. Above 250 hPa, the kinetic energy spectra exhibit a distinct break between steep and shallow spectral ranges, reminiscent of dual power-law spectra seen in aircraft data and high-resolution general circulation models. The break separates a large-scale ‘‘balanced’’ regime in which rotational flow strongly dominates divergent flow and a mesoscale ‘‘unbalanced’’ regime where divergent energy is comparable to or larger than rotational energy. Between 230 and 100 hPa, the spectral break shifts to larger scales (from n 5 60 to n 5 20, where n is spherical harmonic index) as the balanced component of the flow preferentially decays. The location of the break remains fairly stable throughout the stratosphere. The spectral break in the analysis occurs at somewhat larger scales than the break seen in aircraft data. Nonlinear spectral fluxes defined for the rotational component of the flow maximize between about 300 and 200 hPa. Large-scale turbulence thus centers on the extratropical tropopause region, within which there are two distinct mechanisms of upscale energy transfer: eddy–eddy interactions sourcing the transient energy peak in synoptic scales, and zonal mean–eddy interactions forcing the zonal flow. A well-defined downscale enstrophy flux is clearly evident at these altitudes. In the stratosphere, the transient energy peak moves to planetary scales and zonal mean–eddy interactions become dominant.

Comments on some recent measurements of anomalously steep N2O and O3tracer spectra in the stratospheric surf zone

Recent aircraft measurements, primarily in the extratropics, of the horizontal variance of nitrous oxide (N2O) and ozone (O3) in the middle stratosphere indicate that horizontal spectra of the tracer variance scale nearly as k−2, where k is the spatial wavenumber along the aircraft flight track [Strahan and Mahlman, 1994; Bacmeister et al., 1996]. This spectral scaling has been regarded as inconsistent with the accepted picture of stratospheric tracer motion; large-scale quasi-two-dimensional tracer advection typically yields a k−1 scaling (i.e., the classical Batchelor spectrum). In this paper it is argued that the nearly k−2 scaling seen in the measurements is a natural outcome of quasi-two-dimensional filamentation of the polar vortex edge. The accepted picture of stratospheric tracer motion can thus be retained: no additional physical processes are needed to account for deviations from the Batchelor spectrum. Our argument is based on the finite lifetime of tracer filaments and on the “singularity spectrum” associated with a one-dimensional field composed of randomly spaced jumps in concentration.

Differentiation between Neotropical rainforest, dry forest and savannah ecosystems by their modern pollen spectra and implications for the fossil pollen record

Accurate differentiation between tropical forest and savannah ecosystems in the fossil pollen record is hampered by the combination of: i) poor taxonomic resolution in pollen identification, and ii) the high species diversity of many lowland tropical families, i.e. with many different growth forms living in numerous environmental settings. These barriers to interpreting the fossil record hinder our understanding of the past distributions of different Neotropical ecosystems and consequently cloud our knowledge of past climatic, biodiversity and carbon storage patterns. Modern pollen studies facilitate an improved understanding of how ecosystems are represented by the pollen their plants produce and therefore aid interpretation of fossil pollen records. To understand how to differentiate ecosystems palynologically, it is essential that a consistent sampling method is used across ecosystems. However, to date, modern pollen studies from tropical South America have employed a variety of methodologies (e.g. pollen traps, moss polsters, soil samples). In this paper, we present the first modern pollen study from the Neotropics to examine the modern pollen rain from moist evergreen tropical forest (METF), semi-deciduous dry tropical forest (SDTF) and wooded savannah (cerradão) using a consistent sampling methodology (pollen traps). Pollen rain was sampled annually in September for the years 1999–2001 from within permanent vegetation study plots in, or near, the Noel Kempff Mercado National Park (NKMNP), Bolivia. Comparison of the modern pollen rain within these plots with detailed floristic inventories allowed estimates of the relative pollen productivity and dispersal for individual taxa to be made (% pollen/% vegetation or ‘p/v’). The applicability of these data to interpreting fossil records from lake sediments was then explored by comparison with pollen assemblages obtained from five lake surface samples.

The global distribution of phytoplankton size spectrum and size classes from their light-absorption spectra derived from satellite data

The absorption spectra of phytoplankton in the visible domain hold implicit information on the phytoplankton community structure. Here we use this information to retrieve quantitative information on phytoplankton size structure by developing a novel method to compute the exponent of an assumed power-law for their particle-size spectrum. This quantity, in combination with total chlorophyll-a concentration, can be used to estimate the fractional concentration of chlorophyll in any arbitrarily-defined size class of phytoplankton. We further define and derive expressions for two distinct measures of cell size of mixed populations, namely, the average spherical diameter of a bio-optically equivalent homogeneous population of cells of equal size, and the average equivalent spherical diameter of a population of cells that follow a power-law particle-size distribution. The method relies on measurements of two quantities of a phytoplankton sample: the concentration of chlorophyll-a, which is an operational index of phytoplankton biomass, and the total absorption coefficient of phytoplankton in the red peak of visible spectrum at 676 nm. A sensitivity analysis confirms that the relative errors in the estimates of the exponent of particle size spectra are reasonably low. The exponents of phytoplankton size spectra, estimated for a large set of in situ data from a variety of oceanic environments (~ 2400 samples), are within a reasonable range; and the estimated fractions of chlorophyll in pico-, nano- and micro-phytoplankton are generally consistent with those obtained by an independent, indirect method based on diagnostic pigments determined using high-performance liquid chromatography. The estimates of cell size for in situ samples dominated by different phytoplankton types (diatoms, prymnesiophytes, Prochlorococcus, other cyanobacteria and green algae) yield nominal sizes consistent with the taxonomic classification. To estimate the same quantities from satellite-derived ocean-colour data, we combine our method with algorithms for obtaining inherent optical properties from remote sensing. The spatial distribution of the size-spectrum exponent and the chlorophyll fractions of pico-, nano- and micro-phytoplankton estimated from satellite remote sensing are in agreement with the current understanding of the biogeography of phytoplankton functional types in the global oceans. This study contributes to our understanding of the distribution and time evolution of phytoplankton size structure in the global oceans.

Modelling the time-evolution of phytoplankton size spectra from satellite remote sensing

A dynamic size-structured model is developed for phytoplankton and nutrients in the oceanic mixed layer and applied to extract phytoplankton biomass at discrete size fractions from remotely sensed, ocean-colour data. General relationships between cell size and biophysical processes (such as sinking, grazing, and primary production) of phytoplankton were included in the model through a bottom–up approach. Time-dependent, mixed-layer depth was used as a forcing variable, and a sequential data-assimilation scheme was implemented to derive model trajectories. From a given time-series, the method produces estimates of size-structured biomass at every observation, so estimates seasonal succession of individual phytoplankton size, derived here from remote sensing for the first time. From these estimates, normalized phytoplankton biomass size spectra over a period of 9 years were calculated for one location in the North Atlantic. Further analysis demonstrated that strong relationships exist between the seasonal trends of the estimated size spectra and the mixed-layer depth, nutrient biomass, and total chlorophyll. The results contain useful information on the time-dependent biomass flux in the pelagic ecosystem.