A study of the stability of subcycling algorithms in structural dynamics

Algorithms for explicit integration of structural dynamics problems with multiple time steps (subcycling) are investigated. Only one such algorithm, due to Smolinski and Sleith has proved to be stable in a classical sense. A simplified version of this algorithm that retains its stability is presented. However, as with the original version, it can be shown to sacrifice accuracy to achieve stability. Another algorithm in use is shown to be only statistically stable, in that a probability of stability can be assigned if appropriate time step limits are observed. This probability improves rapidly with the number of degrees of freedom in a finite element model. The stability problems are shown to be a property of the central difference method itself, which is modified to give the subcycling algorithm. A related problem is shown to arise when a constraint equation in time is introduced into a time-continuous space-time finite element model. (C) 1998 Elsevier Science S.A.

Anomalous fluorescence spectra in an extremely weak, broadband squeezed vacuum under monochromatic and bichromatic excitation

We show that a two-level atom interacting with an extremely weak squeezed vacuum can display resonance fluorescence spectra that are qualitatively different to those that can be obtained using fields with a classical analogue. We consider first the free space situation with monochromatic excitation, and then discuss a bichromatically driven two-level atom in a cavity as a practical scenario for experimentally detecting the anomalous features predicted. We show that in the bad cavity limit, the anomalous spectral features appear for a weak squeezed vacuum and large frequency differences of the bichromatic field, conditions which are easily accessible in laboratories. The advantage of bichromatic, as opposed to monochromatic, excitation is that there is no coherent scattering at line centre which could obscure the observations. A scaling law is derived, N similar to Omega(4) which relates the squeezed photon number to the Rabi frequency at which the anomalous features appear. (C) 1998 Elsevier Science B.V.

Relationships between blood pressure and heart rate in the saltwater crocodile Crocodylus porosus

The cardiac limb of the baroreflex loop was studied in the saltwater crocodile Crocodylus porosus, The classical pharmacological methodology using phenylephrine and sodium nitroprusside was used to trigger blood pressure changes, and the resulting alterations in heart rate were analysed quantitatively using a logistic function. Interindividual differences in resting heart rates and blood pressures were observed, but all seven animals displayed clear baroreflex responses. Atropine and sotalol greatly attenuated the response. A maximal baroreflex gain of 7.2 beats min(-1) kPa(-1) was found at a mean aortic pressure of 6.1 kPa, indicating the active role of the baroreflex in a wide pressure range encompassing hypotensive and hypertensive states. At the lowest mean aortic pressures (5.0 kPa), the synergistic role of the pulmonary-to-systemic shunt in buffering the blood pressure drop also contributes to blood pressure regulation, Pulse pressure showed a better correlation,vith heart rate and also a higher gain than mean aortic, systolic or diastolic pressures, and this is taken as an indicator of the existence of a differential control element working simultaneously with a linear proportional element.

Interference pattern with a dark center in resonance fluorescence from two atoms driven by a squeezed vacuum field

We study the resonance fluorescence from two interacting atoms driven by a squeezed vacuum field and show that this system produces an interference pattern with a dark center. We discuss the role of the interatomic interactions in this process and find that the interference pattern results from an unequal population of the symmetric and antisymmetric states of the two-atom system. We also identify intrinsically nonclassical effects versus classical squeezed field effects, (C) 1998 Elsevier Science B.V. All rights reserved.

Type 1 and Type 2 cytokines: from basic science to fungal infections

At the present time, it is clear that Th1 responses afford protection against the fungi; however, the development, maintenance and function of the protective immune responses are complex mechanisms and are influenced by multiple factors. The route of infection has been shown to affect initial cytokine production and, consequently, the induction of protective Th1 responses. The ability of different isolates of the same fungal agent to induce and sustain a protective response has also been emphasized. Protective immune responses have been shown to vary in genetically different mouse strains after infection. In addition, these protective responses, such as cellular influx and cytokine production, also vary within the same animal depending on the tissue infected. The functional dominance of certain cytokines over others in influencing development and maintenance of protective responses has been discussed. Certain cytokines may act differently in hosts lacking important components of their innate or immune repertoire. It is evident from these presentations that a more comprehensive understanding of the protective mechanisms against different fungal agents is emerging. However, there is still much to learn before cytokine modulatory therapy can be used effectively without risk in the human host.

Three-dimensional quantum solitons with parametric coupling

We consider the quantum field theory of two bosonic fields interacting via both parametric (cubic) and quartic couplings. In the case of photonic fields in a nonlinear optical medium, this corresponds to the process of second-harmonic generation (via chi((2)) nonlinearity) modified by the chi((3)) nonlinearity. The quantum solitons or energy eigenstates (bound-state solutions) are obtained exactly in the simplest case of two-particle binding, in one, two, and three space dimensions. We also investigate three-particle binding in one space dimension. The results indicate that the exact quantum solitons of this field theory have a singular, pointlike structure in two and three dimensions-even though the corresponding classical theory is nonsingular. To estimate the physically accessible radii and binding energies of the bound states, we impose a momentum cutoff on the nonlinear couplings. In the case of nonlinear optical interactions, the resulting radii and binding energies of these photonic particlelike excitations in highly nonlinear parametric media appear to be close to physically observable values.

Likelihood-based estimation of the regression model with scrambled responses

A significant problem in the collection of responses to potentially sensitive questions, such as relating to illegal, immoral or embarrassing activities, is non-sampling error due to refusal to respond or false responses. Eichhorn & Hayre (1983) suggested the use of scrambled responses to reduce this form of bias. This paper considers a linear regression model in which the dependent variable is unobserved but for which the sum or product with a scrambling random variable of known distribution, is known. The performance of two likelihood-based estimators is investigated, namely of a Bayesian estimator achieved through a Markov chain Monte Carlo (MCMC) sampling scheme, and a classical maximum-likelihood estimator. These two estimators and an estimator suggested by Singh, Joarder & King (1996) are compared. Monte Carlo results show that the Bayesian estimator outperforms the classical estimators in almost all cases, and the relative performance of the Bayesian estimator improves as the responses become more scrambled.

Split tolerance to a viral antigen expressed in thymic epithelium and keratinocytes

When expressed as a transgene from the keratin 14 (K14) promoter in an MHC class II-deficient mouse, I-Ab expressed in thymic cortical epithelium promotes positive but not negative selection of I-Ab-restricted CD4(+) T cells (Laufer, T. M. et al., Nature 1996. 383:81-85). Transgenic mice expressing the E7 protein of human papilloma virus 16 from the K14 promoter were studied to determine the consequence of expression of a cytoplasmic/nuclear protein from the K14 promoter. K14E7-transgenic mice express E7 in the thymus and skin without evidence for autoimmunity to E7. Repeated immunization of FVB(H-2(q)) or F1(C57BV6JxFVB) mice with E7 elicited similar antibody responses to the defined B cell epitopes of E7 in K14E7-transgenic and non-transgenic animals. In contrast, for each genetic background, a single immunization with E7 elicited demonstrable T cell proliferative responses to the major promiscuous T helper epitope of E7 in the transgenic but not the non-transgenic animals. Further,E7-immunized non-transgenic F1 (FVBxC57BL/6J) animals developed strong E7-specific cytotoxic T lymphocyte (CTL) responses and were protected against challenge with E7(+) tumors, whereas similarly immunized K14E7-transgenic animals had a markedly reduced CTL response to E7 and no E7-specific tumor protection was observed, although the antibody and CTL response to ovalbumin was normal. Expression of E7 protein as a transgene from the K14 promoter in the skin and thymus thus induces E7-specific tolerance in the cytotoxic T effector repertoire, together with expansion of the E7-specific T helper repertoire. These findings demonstrate that limited tissue distribution of an autoantigen may result in split tolerance to that autoantigen.

Generalized optimal lattice covering of finite-dimensional Euclidean space

A generalization of the classical problem of optimal lattice covering of R-n is considered. Solutions to this generalized problem are found in two specific classes of lattices. The global optimal solution of the generalization is found for R-2. (C) 1998 Elsevier Science Inc. All rights reserved.

Coherent molecular solitons in Bose-Einstein condensates

We analyze the coherent formation of molecular Bose-Einstein condensate (BEC) from an atomic BEG, using a parametric field theory approach. We point out the transition between a quantum soliton regime, where atoms couple in a local way to a classical soliton domain, where a stable coupled-condensate soliton can form in three dimensions. This gives the possibility of an intense, stable atom-laser output. [S0031-9007(98)07283-4].

Thermal diffusion in molecular dynamics simulations of thin film diamond deposition

Molecular dynamics simulations of carbon atom depositions are used to investigate energy diffusion from the impact zone. A modified Stillinger-Weber potential models the carbon interactions for both sp2 and sp3 bonding. Simulations were performed on 50 eV carbon atom depositions onto the (111) surface of a 3.8 x 3.4 x 1.0 nm diamond slab containing 2816 atoms in 11 layers of 256 atoms each. The bottom layer was thermostated to 300 K. At every 100th simulation time step (27 fs), the average local kinetic energy, and hence local temperature, is calculated. To do this the substrate is divided into a set of 15 concentric hemispherical zones, each of thickness one atomic diameter (0.14 nm) and centered on the impact point. A 50-eV incident atom heats the local impact zone above 10 000 K. After the initial large transient (200 fs) the impact zone has cooled below 3000 K, then near 1000 K by 1 ps. Thereafter the temperature profile decays approximately as described by diffusion theory, perturbed by atomic scale fluctuations. A continuum model of classical energy transfer is provided by the traditional thermal diffusion equation. The results show that continuum diffusion theory describes well energy diffusion in low energy atomic deposition processes, at distance and time scales larger than 1.5 nm and 1-2 ps, beyond which the energy decays essentially exponentially. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

Unsupervised cell nucleus segmentation with active contours

The task of segmenting cell nuclei from cytoplasm in conventional Papanicolaou (Pap) stained cervical cell images is a classical image analysis problem which may prove to be crucial to the development of successful systems which automate the analysis of Pap smears for detection of cancer of the cervix. Although simple thresholding techniques will extract the nucleus in some cases, accurate unsupervised segmentation of very large image databases is elusive. Conventional active contour models as introduced by Kass, Witkin and Terzopoulos (1988) offer a number of advantages in this application, but suffer from the well-known drawbacks of initialisation and minimisation. Here we show that a Viterbi search-based dual active contour algorithm is able to overcome many of these problems and achieve over 99% accurate segmentation on a database of 20 130 Pap stained cell images. (C) 1998 Elsevier Science B.V. All rights reserved.

Evolution of three-dimensional folds for a non-Newtonian plate in a viscous medium

We derive analytical solutions for the three-dimensional time-dependent buckling of a non-Newtonian viscous plate in a less viscous medium. For the plate we assume a power-law rheology. The principal, axes of the stretching D-ij in the homogeneously deformed ground state are parallel and orthogonal to the bounding surfaces of the plate in the flat state. In the model formulation the action of the less viscous medium is replaced by equivalent reaction forces. The reaction forces are assumed to be parallel to the normal vector of the deformed plate surfaces. As a consequence, the buckling process is driven by the differences between the in-plane stresses and out of plane stress, and not by the in-plane stresses alone as assumed in previous models. The governing differential equation is essentially an orthotropic plate equation for rate dependent material, under biaxial pre-stress, supported by a viscous medium. The differential problem is solved by means of Fourier transformation and largest growth coefficients and corresponding wavenumbers are evaluated. We discuss in detail fold evolutions for isotropic in-plane stretching (D-11 = D-22), uniaxial plane straining (D-22 = 0) and in-plane flattening (D-11 = -2D(22)). Three-dimensional plots illustrate the stages of fold evolution for random initial perturbations or initial embryonic folds with axes non-parallel to the maximum compression axis. For all situations, one dominant set of folds develops normal to D-11, although the dominant wavelength differs from the Biot dominant wavelength except when the plate has a purely Newtonian viscosity. However, in the direction parallel to D-22, there exist infinitely many modes in the vicinity of the dominant wavelength which grow only marginally slower than the one corresponding to the dominant wavelength. This means that, except for very special initial conditions, the appearance of a three-dimensional fold will always be governed by at least two wavelengths. The wavelength in the direction parallel to D-11 is the dominant wavelength, and the wavelength(s) in the direction parallel to D-22 is determined essentially by the statistics of the initial state. A comparable sensitivity to the initial geometry does not exist in the classic two-dimensional folding models. In conformity with tradition we have applied Kirchhoff's hypothesis to constrain the cross-sectional rotations of the plate. We investigate the validity of this hypothesis within the framework of Reissner's plate theory. We also include a discussion of the effects of adding elasticity into the constitutive relations and show that there exist critical ratios of the relaxation times of the plate and the embedding medium for which two dominant wavelengths develop, one at ca. 2.5 of the classical Biot dominant wavelength and the other at ca. 0.45 of this wavelength. We propose that herein lies the origin of parasitic folds well known in natural examples.

Integrability of a t - J model with impurities

A t - J model for correlated electrons with impurities is proposed. The impurities are introduced in such a way that integrability of the model in one dimension is not violated. The algebraic Bethe ansatz solution of the model is also given and it is shown that the Bethe states are highest weight states with respect to the supersymmetry algebra gl(2/1).

Atoms in squeezed light fields

Squeezed light is of interest as an example of a non-classical state of the electromagnetic field and because of its applications both in technology and in fundamental quantum physics. This review concentrates on one aspect of squeezed light, namely its application in atomic spectroscopy. The general properties, detection and application of squeezed light are first reviewed. The basic features of the main theoretical methods (master equations, quantum Langevin equations, coupled systems) used to treat squeezed light spectroscopy are then outlined. The physics of squeezed light interactions with atomic systems is dealt with first for the simpler case of two-level atoms and then for the more complex situation of multi-level atoms and multi-atom systems. Finally the specific applications of squeezed light spectroscopy are reviewed.