Numerical Study of Hypersonic Leeward Flow over a Blunt Nosed Delta Wing

A kinetic theory based Navier-Stokes solver has been implemented on a parallel supercomputer (Intel iPSC Touchstone Delta) to study the leeward flowfield of a blunt nosed delta wing at 30-deg incidence at hypersonic speeds (similar to the proposed HERMES aerospace plane). Computational results are presented for a series of grids for both inviscid and laminar viscous flows at Reynolds numbers of 225,000 and 2.25 million. In addition, comparisons are made between the present and two independent calculations of the some flows (by L. LeToullec and P. Guillen, and S. Menne) which were presented at the Workshop on Hypersonic Flows for Re-entry Problems, Antibes, France, 1991.

Structure-activity relationships of alpha-conotoxins targeting neuronal nicotinic acetylcholine receptors

alpha-Conotoxins that target the neuronal nicotinic acetylcholine receptor have a range of potential therapeutic applications and are valuable probes for examining receptor subtype selectivity. The three-dimensional structures of about half of the known neuronal specific alpha-conotoxins have now been determined and have a consensus fold containing a helical region braced by two conserved disulfide bonds. These disulfide bonds define the two-loop framework characteristic for alpha-conotoxins, CCXmCXnC, where loop 1 comprises four residues (m = 4) and loop 2 between three and seven residues (n = 3, 6 or 7). Structural studies, particularly using NMR spectroscopy have provided an insight into the role and spatial location of residues implicated in receptor binding and biological activity.

Molecular orbital calculations of two-electron states for P-donor solid-state spin qubits

We theoretically study the Hilbert space structure of two neighboring P-donor electrons in silicon-based quantum computer architectures. To use electron spins as qubits, a crucial condition is the isolation of the electron spins from their environment, including the electronic orbital degrees of freedom. We provide detailed electronic structure calculations of both the single donor electron wave function and the two-electron pair wave function. We adopted a molecular orbital method for the two-electron problem, forming a basis with the calculated single donor electron orbitals. Our two-electron basis contains many singlet and triplet orbital excited states, in addition to the two simple ground state singlet and triplet orbitals usually used in the Heitler-London approximation to describe the two-electron donor pair wave function. We determined the excitation spectrum of the two-donor system, and study its dependence on strain, lattice position, and interdonor separation. This allows us to determine how isolated the ground state singlet and triplet orbitals are from the rest of the excited state Hilbert space. In addition to calculating the energy spectrum, we are also able to evaluate the exchange coupling between the two donor electrons, and the double occupancy probability that both electrons will reside on the same P donor. These two quantities are very important for logical operations in solid-state quantum computing devices, as a large exchange coupling achieves faster gating times, while the magnitude of the double occupancy probability can affect the error rate.

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

Matrix elements of U(2n) generators in a multishell spin-orbit basis. II. The two-shell nonzero shift ACCs and U(2n) generator MEs

This is the second in a series of articles whose ultimate goal is the evaluation of the matrix elements (MEs) of the U(2n) generators in a multishell spin-orbit basis. This extends the existing unitary group approach to spin-dependent configuration interaction (CI) and many-body perturbation theory calculations on molecules to systems where there is a natural partitioning of the electronic orbital space. As a necessary preliminary to obtaining the U(2n) generator MEs in a multishell spin-orbit basis, we must obtain a complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. The zero-shift coefficients were obtained in the first article of the series. in this article, we evaluate the nonzero shift adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. We then demonstrate that the one-shell versions of these coefficients may be obtained by taking the Gelfand-Tsetlin limit of the two-shell formulas. These coefficients,together with the zero-shift types, then enable us to write down formulas for the U(2n) generator matrix elements in a two-shell spin-orbit basis. Ultimately, the results of the series may be used to determine the many-electron density matrices for a partitioned system. (C) 1998 John Wiley & Sons, Inc.

Quantum trajectory simulations of the fluorescence intensity from a two-level atom driven by a multichromatic field

The quantum trajectories method is illustrated for the resonance fluorescence of a two-level atom driven by a multichromatic field. We discuss the method for the time evolution of the fluorescence intensity in the presence of bichromatic and trichromatic driving fields. We consider the special case wherein one multichromatic field component is strong and resonant with the atomic transition whereas the other components are much weaker and arbitrarily detuned from the atomic resonance. We find that the phase-dependent modulations of the Rabi oscillations, recently observed experimentally [Q. Wu, D. J. Gauthier, and T. W. Mossberg, Phys. Rev. A 49, R1519 (1994)] for the special case when the weaker component of a bichromatic driving field is detuned from the atomic resonance by the strong-field Rabi frequency, appear also for detunings close to the subharmonics of the Rabi frequency. Furthermore, we show that for the atom initially prepared in one of the dressed states of the strong field component the modulations are not sensitive to the phase. We extend the calculations to the case of a trichromatic driving field and find that apart from the modulations of the amplitude there is a modulation of the frequency of the Rabi oscillations. Moreover, the time evolution of the fluorescence intensity depends on the phase regardless of the initial conditions and a phase-dependent suppression of the Rabi oscillations can be observed when the sideband fields are tuned to the subharmonics of the strong-field Rabi frequency. [S1050-2947(98)03501-X].

Structure-activity relationships of omega-conotoxins MVIIA, MVIIC and 14 loop splice hybrids at N and P/Q-type calcium channels

The omega-conotoxins are a set of structurally related, four-loop, six cysteine containing peptides, that have a range of selectivities for different subtypes of the voltage-sensitive calcium channel (VSCC). To investigate the basis of the selectivity displayed by these peptides, we have studied the binding affinities of two naturally occurring omega-conotoxins, MVIIA and MVIIC and a series of 14 MVIIA/MVIIC loop hybrids using radioligand binding assays for N and P/Q-type Ca2+ channels in rat brain tissue. A selectivity profile was developed from the ratio of relative potencies at N-type VSCCs (using [I-125]GVIA radioligand binding assays) and P/Q-type VSCCs (using [I-125]MVIIC radioligand binding assays). in these peptides, loops 2 and 4 make the greatest contribution to VSCC subtype selectivity, while the effects of loops 1 and 3 are negligible. Peptides with homogenous combinations of loop 2 and 4 display clear selectivity preferences, while those with heterogeneous combinations of loops 2 and 4 are less discriminatory. H-1 NMR spectroscopy revealed that the global folds of MVIIA, MVIIC and the 14 loop hybrid peptides were similar; however, several differences in local structure were identified. Based on the binding data and the 3D structures of MVIIA, GVIA and MVIIC, we have developed a preliminary pharmacophore based on the omega-conotoxin residues most Likely to interact with the N-type VSCC. (C) 1999 Academic Press.

Solution structure of BSTI: A new trypsin inhibitor from skin secretions of Bombina bombina

The three-dimensional solution structure of BSTI, a trypsin inhibitor from the European frog Bombina bombina, has been solved using H-1 NMR spectroscopy. The 60 amino acid protein contains five disulfide bonds, which were unambiguously determined to be Cvs (4-38), Cys (13-34), Cys (17-30), Cys (21-60), and Cys (40-54) by experimental restraints and subsequent structure calculations. The main elements of secondary structure are four beta -strands, arranged as two small antiparallel beta -sheets, The overall fold of BSTI is disk shaped and is characterized by the lack of a hydrophobic core. The presumed active site is located on a loop comprising residues 21-34, which is a relatively disordered region similar to that seen in many other protease inhibitors. However, the overall fold is different to other known protease inhibitors with the exception of a small family of inhibitors isolated from nematodes of the family Ascaris and recently also from the haemolymph of Apis mellifera. BSTI may thus be classified as a new member of this recently discovered family of protease inhibitors.

Comments on alternative calculations of the broadening of spectral lines of neutral sodium by H-atom collisions

With the exception of the sodium D-lines, recent calculations of line broadening cross sections for several multiplets of sodium by Leininger et al (Leininger T, Gadea F X and Dickinson A 2000 J. Phys. B: At. Mol. Opt. Phys. 33 1805) are in substantial disagreement with cross sections interpolated from the tables of Anstee and O'Mara (Anstee and O'Mara 1995 Mon. Not. R. Astron. Soc. 276 859) and Barklem and O'Mara (Barklem P S and O'Mara B J 1997 Mon. Not. R. Astron. Soc. 290 102). The discrepancy is as large as a factor of 3 for the 3p-4d multiplet. The two theories are tested by using the results of each to synthesize lines in the solar spectrum. It is found that generally the data from the theory of Anstee, Barklem and O'Mara produce the best match to the observed solar spectrum. It is found, using a simple model for reflection of the optical electron by the potential barrier between the two atoms, that the reflection coefficient is too large for avoided crossings with the upper states of subordinate lines to contribute to line broadening, supporting the neglect of avoided ionic crossings by Anstee, Barklem and O'Mara for these lines. The large discrepancies between the two sets of calculations is a result of an approximate treatment of avoided ionic crossings for these lines by Leininger et al (Leininger T, Gadea F X and Dickinson A 2000 J. Phys. B: At. Mol. Opt. Phys. 33 1805).

Modulation of ligand selectivity by mutation of the first extracellular loop of the human C5a receptor

The cyclic C5a receptor antagonist, phenylalanine [L-ornithine-proline-D-cyclohexylalanine-tryptophan-arginine] (F-[OPchaWR]), has similar to 1000-fold less affinity for the C5a receptor (C5aR) on murine polymorphonuclear leukocytes than on human. Analysis of C5aR from different species shows that a possible cause of this difference is the variation in the sequence of the first extracellular loop of the receptor. The mouse receptor contains Y at a position analogous to P-103 in the human receptor, and D at G(105). To test this hypothesis, we expressed human C5aR mutants ((PY)-Y-103, G(105)D and the double mutant, (PY)-Y-103/G(105)D) in RBL-2H3 cells and investigated the effects of these mutations on binding affinity and receptor activation. All three mutant receptors had a higher affinity for human C5a than the wild-type receptor, but showed no significant difference in the ability of F-[OPchaWR] to inhibit human C5a binding. However, all of the mutant receptors had substantially lower affinities for the weak agonist, C5a des Arg(74) (C5adR(74)), and two altered receptors (G(105)D and (PY)-Y-103/G(105)D) had much lower affinities for the C-terminal C5a agonist peptide analogue, L-tyrosine-serine-phenylalanine-lysine-proline-methionine-proline-leucine-D-alanine-arginine (YSFKPMPLaR). Although it is unlikely that differences at these residues are responsible for variations in the potency of F-[OPchaWR] across species, residues in the first extracellular loop are clearly involved in the recognition of both C5a and C5a agonists. The complex effects of mutating these residues on the affinity and response to C5a, C5adR(74), and the peptide analogues provide evidence of different binding modes for these ligands on the C5aR. (C) 2001 Elsevier Science Inc. All rights reserved.

Shear stress and sediment transport calculations for sheet flow under waves

A simple method is provided for calculating transport rates of not too fine (d(50) greater than or equal to 0.20 mm) sand under sheet flow conditions. The method consists of a Meyer-Peter-type transport formula operating on a time-varying Shields parameter, which accounts for both acceleration-asymmetry and boundary layer streaming. While velocity moment formulae, e.g.., = Constant x calibrated against U-tube measurements, fail spectacularly under some real waves (Ribberink, J.S., Dohmen-Janssen, C.M., Hanes, D.M., McLean, S.R., Vincent, C., 2000. Near-bed sand transport mechanisms under waves. Proc. 27th Int. Conf. Coastal Engineering, Sydney, ASCE, New York, pp. 3263-3276, Fig. 12), the new method predicts the real wave observations equally well. The reason that the velocity moment formulae fail under these waves is partly the presence of boundary layer streaming and partly the saw-tooth asymmetry, i.e., the front of the waves being steeper than the back. Waves with saw-tooth asymmetry may generate a net landward sediment transport even if = 0, because of the more abrupt acceleration under the steep front. More abrupt accelerations are associated with thinner boundary layers and greater pressure gradients for a given velocity magnitude. The two real wave effects are incorporated in a model of the form Q(s)(t) = Q(s)[theta(t)] rather than Q(S)(t) = Q(S)[u(infinity)(t)], i.e., by expressing the transport rate in terms of an instantaneous Shields parameter rather than in terms of the free stream velocity, and accounting for both streaming and accelerations in the 0(t) calculations. The instantaneous friction velocities u(*)(t) and subsequently theta(t) are calculated as follows. Firstly, a linear filter incorporating the grain roughness friction factor f(2.5) and a phase angle phi(tau) is applied to u(infinity)(t). This delivers u(*)(t) which is used to calculate an instantaneous grain roughness Shields parameter theta(2.5)(t). Secondly, a constant bed shear stress is added which corresponds to the streaming related bed shear stress -rho ($) over bar((u) over tilde(w) over tilde)(infinity) . The method can be applied to any u(infinity)(t) time series, but further experimental validation is recommended before application to conditions that differ strongly from the ones considered below. The method is not recommended for rippled beds or for sheet flow with typical prototype wave periods and d(50) < 0.20 turn. In such scenarios, time lags related to vertical sediment movement become important, and these are not considered by the present model. (C) 2002 Elsevier Science B.V. All rights reserved.

Chebyshev real wave packet propagation: H+O-2 (J=0) state-to-state reactive scattering calculations

In this paper we explore the relative performance of two recently developed wave packet methodologies for reactive scattering, namely the real wave packet Chebyshev domain propagation of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)] and the Lanczos subspace wave packet approach of Smith [J. Chem. Phys. 116, 2354 (2002); Chem. Phys. Lett. 336, 149 (2001)]. In the former method, a modified Schrodinger equation is employed to propagate the real part of the wave packet via the well-known Chebyshev iteration. While the time-dependent wave packet from the modified Schrodinger equation is different from that obtained using the standard Schrodinger equation, time-to-energy Fourier transformation yields wave functions which differ only trivially by normalization. In the Lanczos subspace approach the linear system of equations defining the action of the Green operator may be solved via either time-dependent or time-independent methods, both of which are extremely efficient due to the simple tridiagonal structure of the Hamiltonian in the Lanczos representation. The two different wave packet methods are applied to three dimensional reactive scattering of H+O-2 (total J=0). State-to-state reaction probabilities, product state distributions, as well as initial-state-resolved cumulative reaction probabilities are examined. (C) 2002 American Institute of Physics.

A comparative study of carbon gasification with O-2 and CO2 by density functional theory calculations

A comparative study of carbon gasification with O-2 and CO2 was conducted by using density functional theory calculations. It was found that the activation energy and the number of active sites in carbon gasification reactions are significantly affected by both the capacity and manner of gas chemisorption. O-2 has a strong adsorption capacity and the dissociative chemisorption of O-2 is thermodynamically favorable on either bare carbon surface or even isolated edge sites. As a result, a large number of semiquinone and o-quinone oxygen can be formed indicating a significant increase in the number of active sites. Moreover, the weaker o-quinone C-C bonds can also drive the reaction forward at (ca. 30%) lower activation energy. Epoxy oxygen forms under relatively high O-2 pressure, and it can only increase the number of active sites, not further reduce the activation energy. CO2 has a lower adsorption capacity. Dissociative chemisorption of CO2 can only occur on two consecutive edge sites and o-quinone oxygen formed from CO2 chemisorption is negligible, let alone epoxy oxygen. Therefore, CO2-carbon reaction needs (ca 30%) higher activation energy. Furthermore, the effective active sites are also reduced by the manner Of CO2 chemisorption. A combination of the higher activation energy and the fewer active sites leads to the much lower reaction rate Of CO2-carbon.

Secondary structure of the third extracellular loop responsible for ligand selectivity of a mammalian gonadotropin-releasing hormone receptor

The extracellular loop 3 (ECL3) of the mammalian gonadotropin-releasing hormone receptor (GnRH-R) contains an acidic amino acid (Glu(301) in the mouse GnRH-R,) that confers agonist selectivity for Are in mammalian GnRH. It is proposed that a specific conformation of ECL3 is necessary to orientate the carboxyl side chain of the acidic residue for interaction with Arg(8) of GnRH, which is supported by decreased affinity for Arg(8) GnRH but not Gln(8) GnRH when an adjacent Pro is mutated to Ala. To probe the structural contribution of the loop domain to the proposed presentation of the carboxyl side chain, we synthesized a model peptide (CGPEMLNRVSEPGC) representing residues 293-302 of mouse ECL3, where Cys and Gly residues are added symmetrically at the N and C termini, respectively, allowing the introduction of a disulfide bridge to simulate the distances at which the ECL3 is tethered to the transmembrane domains 6 and 7 of the receptor. The ability of the ECL3 peptide to bind GnRH with low affinity was demonstrated by its inhibition of GnRH stimulation of inositol phosphate production in cells expressing the GnRH-R. The CD bands of the ECL3 peptides exhibited a superposition of predominantly unordered structure and partial contributions from beta-sheet structure. Likewise, the analysis of the amide I and amide III bands from micro-Raman and FT Raman experiments revealed mainly unordered conformations of the cyclic and of the linear peptide. NMR data demonstrated the presence of a beta-hairpin among an ensemble of largely disordered structures in the cyclic peptide. The location of the turn linking the two strands of the hairpin was assigned to the three central residues L-296, N-297, and R-298. A small population of structured species among an ensemble of predominantly random coil conformation suggests that the unliganded receptor represents a variety of structural conformers, some of which have the potential to make contacts with the ligand. We propose a mechanism of receptor activation whereby binding of the agonist to the inactive receptor state induces and stabilizes a particular structural state of the loop domain, leading to further conformational rearrangements across the transmembrane domain and signal propagating interaction with G proteins. Interaction of the Glu(301) of the receptor with Arg(8) of GnRH induces a folded configuration of the ligand. Our proposal thus suggests that conformational changes of both ligand and receptor result from this interaction.

A comparison of two methods for the calculation of accumulated oxygen deficit

The aim of this study was to compare accumulated oxygen deficit data derived using two different exercise protocols with the aim of producing a less time-consuming test specifically for use with athletes. Six road and four track male endurance cyclists performed two series of cycle ergometer tests. The first series involved five 10 min sub-maximal cycle exercise bouts, a (V) over dotO(2peak) test and a 115% (V) over dotO(2peak) test. Data from these tests were used to estimate the accumulated oxygen deficit according to the calculations of Medbo et al. (1988). In the second series of tests, participants performed a 15 min incremental cycle ergometer test followed, 2 min later, by a 2 min variable resistance test in which they completed as much work as possible while pedalling at a constant rate. Analysis revealed that the accumulated oxygen deficit calculated from the first series of tests was higher (P< 0.02) than that calculated from the second series: 52.3 +/- 11.7 and 43.9 +/- 6.4 ml . kg(-1), respectively (mean +/- s). Other significant differences between the two protocols were observed for (V) over dot O-2peak, total work and maximal heart rate; all were higher during the modified protocol (P