We simulate oil spills of 1500 and 4500m3/day lasting 14, 45, and 90days in the spawning grounds of the commercial fish species, Northeast Arctic cod. Modeling the life history of individual fish eggs and larvae, we predict deviations from the historical pattern of recruitment to the adult population due to toxic oil exposures. Reductions in survival for pelagic stages of cod were 0-10%, up to a maximum of 43%. These reductions resulted in a decrease in adult cod biomass of <3% for most scenarios, up to a maximum of 12%. In all simulations, the adult population remained at full reproductive potential with a sufficient number of juveniles surviving to replenish the population. The diverse age distribution helps protect the adult cod population from reductions in a single year's recruitment after a major oil spill. These results provide insights to assist in managing oil spill impacts on fisheries.
The authors present a method for calculating the electrostatic potential directly in a straightforward manner. While traditional methods for calculating the electrostatic potential usually involve solving the Poisson equation iteratively, the authors obtain the electrostatic interaction potential by performing direct numerical integration of the Coulomb-law expression using finite-element functions defined on a grid. The singularity of the Coulomb operator is circumvented by an integral transformation and the resulting auxiliary integral is obtained using Gaussian quadrature. The three-dimensional finite-element basis is constructed as a tensor (outer) product of one-dimensional functions, yielding a partial factorization of the expressions. The resulting algorithm has, without using any prescreening or other computational tricks, a formal computational scaling of O(N4∕3), where N is the size of the grid. The authors show here how to implement the method for efficiently running on parallel computers. The matrix multiplications of the innermost loops are completely independent, yielding a parallel algorithm with the computational costs scaling practically linearly with the number of processors.
Global currents: Quantum chemical methodology has been used to investigate the magnetically induced currents of both neutral Buckminsterfullerene and its spherically aromatic +10 ion and has allowed quantitative measurement of the current strengths. Global sphere currents (see picture), instead of local ring currents, are found to be the defining feature of fullerene aromaticity. Supporting information for this article is available on the WWW under http://www.wiley-vch.de/contents/jc_2002/2005/z462348_s.pdf or from the author. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.
The nuclear magnetic shieldings and magnetically induced ring currents have been calculated for the planar ring-shaped hydrogen fluoride trimer (HF)3 at correlated ab initio and density functional theory levels. Calculations of the magnetically induced current densities using the gauge-including magnetically induced current (GIMIC) method show that, contrary to a recent suggestion, (HF)3 has, at the MP2/TZVPP level, a very small ring-current susceptibility of 0.37 nA/T. Thus, only a weak net current is passing across the H···F hydrogen bond. An external magnetic field perpendicular to the ring plane induces strong edge currents circling around each HF molecule giving rise to a nonvanishing magnetic shielding at the center of the ring. The GIMIC results are supported by calculations of the long-range magnetic shielding function; the long-range magnetic shielding is very small, indicating that the magnetically induced ring-current is very weak. The surprisingly large nucleus-independent chemical shift (NICS) value for (HF)3 was recently taken as an indication of "H-bonded aromaticity". The NICS value calculated at the CCSD/QZ2P level is 2.77 ppm. The present GIMIC and aromatic ring-current shielding study shows that some care has to be taken when using NICS values as aromaticity indices.
An entry from the Cambridge Structural Database, the world’s repository for small molecule crystal structures. The entry contains experimental data from a crystal diffraction study. The deposited dataset for this entry is freely available from the CCDC and typically includes 3D coordinates, cell parameters, space group, experimental conditions and quality measures.
Gas-phase NMR spectra of (11)B, (10)B, and (19)F in BF(3) are reported, and high-level ab initio calculations of the corresponding NMR shielding constants are described. Extrapolation of the measured resonance frequencies to the zero-density limit ensures that the results correspond to the ab initio values for an isolated molecule. Simultaneous measurements of (3)He resonance frequencies and application of the calculated shielding constants allow us to determine improved values of the nuclear magnetic dipole moments of (11)B and (10)B. The magnetic moments of both isotopes are also determined independently by comparing with the (19)F spectral parameters (frequencies and shielding constants). The separately derived nuclear magnetic moments are in good agreement, whereas the literature moments of both (11)B and (10)B are noticeably less accurate.
Magnetically induced current densities, nuclear magnetic shieldings, and electric polarizabilities of planar ring-shaped hydrocarbons have been studied at the density-functional theory level using the Becke−Perdew (BP86) functional. The current densities were calculated using the Gauge-Including Magnetically Induced Current (GIMIC) method employing gauge-including atomic orbitals. The GIMIC calculations yield rules to estimate the global and local ring−current strengths as well as the current pathways for the hydrocarbon nanorings. For the overall antiaromatic molecules, aromatic groups such as benzene, naphthalene, anthracene, and pyrene moieties localize the ring current making the global ring currents vanish. The ability of the edge groups to localize the currents is related to the aromatic character of the molecule as a whole. The local ring current prefers to follow the edges of the group. Phenalenyl corner moieties are found to introduce strong global ring currents, whereas with fused benzene and pyrene corner groups the global ring current vanishes. Fused benzene rings in the corner or along the edge of overall antiaromatic molecules sustain local ring currents of about the same size as for a free benzene molecule. For the overall aromatic molecules, the global ring current is split along the bonds of the edge moieties, but the detailed division fulfilling Kirchhoff's current law is not easily predictable and must be calculated for each individual bond. At the phenalenyl corner moieties, the global ring current follows the innermost route isolating the rest of the group from the main delocalization pathway. A hydrocarbon nanoring sustaining strong ring currents should be large and formally aromatic with many and large aromatic moieties along the edges. A clear correlation between the strength of the global ring currents and the size of the electric polarizabilities is obtained. The calculated 1H NMR shieldings of a proton in immediate contact to the global ring current vary between 22 ppm and 67 ppm in the studied molecules. The trend correlates well with the global ring−current strengths, which are in the range of 0−88 nA/T. The 13C NMR shieldings are also sensitive to the strength of the global ring current, but they vary less systematically and are not as good an indicator of the current strength as the hydrogen shieldings.
Coinage and alkali metal mixed clusters, M4Na- (M = Cu, Au) have been investigated experimentally using photoelectron spectroscopy and computationally at correlated ab initio levels. The related Cu4Li-, Ag4Li-, Ag4Na-, and Au4Li- clusters as well as the neutral Cu4Li2 and Cu4Na2 clusters have also been studied computationally. The calculations show that the two lowest isomers of the negatively charged clusters include a pyramidal C4v structure and a planar C2v species. For Cu4Li- and Cu4Na-, the C4v structure is calculated at correlated ab initio level to be 30.9 and 16.9 kJ/mol below the planar C2v isomer, whereas the planar isomers of Au4Li- and Au4Na- are found to be 29.7 and 49.4 kJ/mol below the pyramidal ones. For Ag4Li- and Ag4Na-, the pyramidal isomers are the lowest ones. Comparison of the calculated and measured photoelectron spectra of Cu4Na- and Au4Na- shows that the pyramidal Cu4Na- cluster of C4v symmetry and the planar Au4Na- of C2v symmetry are detected experimentally. Calculations of the magnetically induced current density in Cu4Li- and Cu4Li2 using the Gauge-Including Magnetically Induced Current (GIMIC) method show that strong ring currents are sustained mainly by the highest-occupied molecular orbital primarily derived from the Cu 4s. The GIMIC calculations thus show that the Cu4(2-) ring is -aromatic and that the d orbitals do not play any significant role for the electron delocalization effects. The present study does not support the notion that the square-planar Cu4(2-) is the first example of d-orbital aromatic molecules.
