Charge Density Analysis of Heterohalogen (Cl center dot center dot center dot F) and Homohalogen (F center dot center dot center dot F) Intermolecular Interactions

Experimental charge density distribution in 2-chloro-4-fluorobenzoic acid and 4-fluorobenzamide has been carried out using high resolution X-ray diffraction data collected at 100 K using Hansen-Coppens multipolar formalism of electron density. These compounds display short Cl center dot center dot center dot F and F center dot center dot center dot F interactions, respectively. The experimental results are compared with the theoretical charge densities using theoretical structure factors obtained from periodic quantum calculation at the B3LYP/6-31G** level. The topological features were derived from Bader's ``atoms in molecules'' (AIM) approach. Intermolecular Cl center dot center dot center dot F interaction in 2-chloro-4-fluorobenzoic acid is attractive in nature (type II interaction) while the nature of F center dot center dot center dot F interactions in 4-fluorobenzamide shows indication of a minor decrease in repulsion (type I interaction), though the extent of polarization on the fluorine atom is arguably small.

Scroll-Wave Dynamics in Human Cardiac Tissue: Lessons from a Mathematical Model with Inhomogeneities and Fiber Architecture

Cardiac arrhythmias, such as ventricular tachycardia (VT) and ventricular fibrillation (VF), are among the leading causes of death in the industrialized world. These are associated with the formation of spiral and scroll waves of electrical activation in cardiac tissue; single spiral and scroll waves are believed to be associated with VT whereas their turbulent analogs are associated with VF. Thus, the study of these waves is an important biophysical problem. We present a systematic study of the combined effects of muscle-fiber rotation and inhomogeneities on scroll-wave dynamics in the TNNP (ten Tusscher Noble Noble Panfilov) model for human cardiac tissue. In particular, we use the three-dimensional TNNP model with fiber rotation and consider both conduction and ionic inhomogeneities. We find that, in addition to displaying a sensitive dependence on the positions, sizes, and types of inhomogeneities, scroll-wave dynamics also depends delicately upon the degree of fiber rotation. We find that the tendency of scroll waves to anchor to cylindrical conduction inhomogeneities increases with the radius of the inhomogeneity. Furthermore, the filament of the scroll wave can exhibit drift or meandering, transmural bending, twisting, and break-up. If the scroll-wave filament exhibits weak meandering, then there is a fine balance between the anchoring of this wave at the inhomogeneity and a disruption of wave-pinning by fiber rotation. If this filament displays strong meandering, then again the anchoring is suppressed by fiber rotation; also, the scroll wave can be eliminated from most of the layers only to be regenerated by a seed wave. Ionic inhomogeneities can also lead to an anchoring of the scroll wave; scroll waves can now enter the region inside an ionic inhomogeneity and can display a coexistence of spatiotemporal chaos and quasi-periodic behavior in different parts of the simulation domain. We discuss the experimental implications of our study.

Flap-lag-torsion stability in hover and forward flight with three-dimensional wake

The contributions of full-wake dynamics in trim analysis are demonstrated for finding the control inputs and periodic responses simultaneously, as well as in Floquet eigenanalysis for finding the damping levels. The equations of flap bending, lag bending, and torsion are coupled with a three-dimensional, finite state wake, and low-frequency (<1/rev) to high frequency (>1/rev) multiblade modes are considered. Full blade-wake dynamics is used in trim analysis and Floquet eigenanalysis. A uniform cantilever blade in trimmed flight is investigated over a range of thrust levels, advance ratios, number of blades, and blade torsional frequencies. The investigation includes the convergence characteristics of control inputs, periodic responses, and damping levels with respect to the number of spatial azimuthal harmonics and radial shape functions in the wake representation. It also includes correlation with the measured lag damping of a three-bladed untrimmed rotor. The parametric study shows the dominant influence of wake dynamics on control inputs, periodic responses, and damping levels, and wake theory generally improves the correlation.

Neural network modeling of associative memory: Beyond the Hopfield model

A number of neural network models, in which fixed-point and limit-cycle attractors of the underlying dynamics are used to store and associatively recall information, are described. In the first class of models, a hierarchical structure is used to store an exponentially large number of strongly correlated memories. The second class of models uses limit cycles to store and retrieve individual memories. A neurobiologically plausible network that generates low-amplitude periodic variations of activity, similar to the oscillations observed in electroencephalographic recordings, is also described. Results obtained from analytic and numerical studies of the properties of these networks are discussed.

Parallel computing concepts and methods for floquet analysis of helicopter trim and stability

Floquet analysis is widely used for small-order systems (say, order M < 100) to find trim results of control inputs and periodic responses, and stability results of damping levels and frequencies, Presently, however, it is practical neither for design applications nor for comprehensive analysis models that lead to large systems (M > 100); the run time on a sequential computer is simply prohibitive, Accordingly, a massively parallel Floquet analysis is developed with emphasis on large systems, and it is implemented on two SIMD or single-instruction, multiple-data computers with 4096 and 8192 processors, The focus of this development is a parallel shooting method with damped Newton iteration to generate trim results; the Floquet transition matrix (FTM) comes out as a byproduct, The eigenvalues and eigenvectors of the FTM are computed by a parallel QR method, and thereby stability results are generated, For illustration, flap and flap-lag stability of isolated rotors are treated by the parallel analysis and by a corresponding sequential analysis with the conventional shooting and QR methods; linear quasisteady airfoil aerodynamics and a finite-state three-dimensional wake model are used, Computational reliability is quantified by the condition numbers of the Jacobian matrices in Newton iteration, the condition numbers of the eigenvalues and the residual errors of the eigenpairs, and reliability figures are comparable in both the parallel and sequential analyses, Compared to the sequential analysis, the parallel analysis reduces the run time of large systems dramatically, and the reduction increases with increasing system order; this finding offers considerable promise for design and comprehensive-analysis applications.

Singularity structure and chaotic behavior of the homopolar disk dynamo

We study in great detail a system of three first-order ordinary differential equations describing a homopolar disk dynamo (HDD). This system displays a large variety of behaviors, both regular and chaotic. Existence of periodic solutions is proved for certain ranges of parameters. Stability criteria for periodic solutions are given. The nonintegrability aspects of the HDD system are studied by investigating analytically the singularity structure of the system in the complex domain. Coexisting attractors (including period-doubling sequence) and coexisting strange attractors appear in some parametric regimes. The gluing of strange attractors and the ungluing of a strange attractor are also shown to occur. A period of bifurcation leading to chaos, not observed for other chaotic systems, is shown to characterize the chaotic behavior in some parametric ranges. The limiting case of the Lorenz system is also studied and is related to HDD.

Time-Resolved Photoluminescence and Microwave Conductivity at Semiconductor Electrodes: Depletion Layer Effects

Analytical expressions which include depletion layer effects on low-injection carrier relaxation are being presented for the first time here. Starting from the continuity equation for the minority carriers, we derive expressions for the output signal pertinent to time-resolved microwave and luminescence experiments. These are valid for the time domain that usually overlaps with the time scales of surface processes, such as charge transfer and trapping. Apart from the usual pulse form of illumination, theoretical expressions pertaining to other forms of illumination such as switch-on and switch-off transient modes, a periodic mode, and a steady state and their various inter-relationships are derived here. The expressions obtained are seen to be generalizations of existing flat-band low-injection results in the Limit of early or initial band bendings. The importance of the depletion layer as an experimental parameter is clearly seen in the limit of larger band bendings wherein it is shown, unlike the flat-band case, to exhibit pure exponential forms of carrier relaxation. Our results are consistent with the main conclusions of the numerical and experimental work published recently. Furthermore, this work provides the actual functional relationships between the applied potential and observed carrier decay. This should enable one to extract the surface kinetic parameters, after deciding on the dominant mode of carrier relaxation at the interface, whether charge transfer or trapping, by studying the potential dependence of the fate of relaxation.

Plate-shaped transformation products in zirconium-base alloys

Plate-shaped products resulting from martensitic, diffusional, and mixed mode transformations in zirconium-base alloys are compared. in the present study. These alloys are particularly suitable for the comparison in view of the fact that the lattice correspondence between the parent beta (bcc) and the product alpha (hcp) or gamma-hydride (fct) phases are remarkably similar for different types of transformations. Crystallographic features such as orientation relations, habit planes, and interface structures associated with these transformations have been compared:, with a view toward examining whether the transformation mechanisms have characteristic imprints on these experimental observables. Martensites exhibiting dislocated lath, internally twinned plate, and self-accommodating three-plate cluster morphologies have been encountered in Zr-2.5Nb alloy. Habit planes corresponding to all these morphologies have been found to be consistent with the predictions based on the invariant plane strain (IFS) criterion. Different morphologies have been found to reflect the manner in which the neighboring martensite variants are assembled. Lattice-invariant shears (LISs) for all these cases have been identified to be either {10 (1) over bar 1}(alpha) ((1) over bar 123)(alpha) slip or twinning on (10 (1) over bar 1)(alpha) planes. Widmanstatten alpha precipitates, forming in a step-quenching treatment, have been shown to have a lath morphology, the alpha/beta interface being decorated with a periodic array of (c + a) dislocations at a spacing of 8 to 10 nm. The line vectors of these dislocations are nearly parallel to the invariant lines. The alpha precipitates, forming in the retained beta phase on aging, exhibit an internally twinned structure with a zigzag habit plane. Average habit planes for the morphologies have been found to lie near the {103}(beta) - {113}(beta) poles, which are close to the specific variant of the {112}(beta) plane, which transforms into a prismatic plane of the type {1 (1) over bar 00}(alpha). The crystallography of the formation of the gamma-hydride phase (fct) from both the alpha and beta phases is seen to match the IFS predictions. While the beta-gamma transformation can be treated approximately as a simple shear on the basal plane involving a change in the stacking sequence, the alpha-gamma transformation call be conceptually broken into a alpha --> beta transformation following the Burgers correspondence and the simple beta-gamma shear process. The active eutectoid decomposition in the Zr-Cu system, beta --> alpha + beta', has been described in terms of cooperative growth of the alpha phase from the beta phase through the Burgers correspondence and of the partially ordered beta' (structurally similar to the equilibrium Zr2Cu phase) through an ordering process. Similarities and differences in crystallographic features of these transformations have been discussed. and the importance of the invariant line vector in deciding the geometry of the corresponding habit planes has been pointed out.

Halogen Bonding in 2,5-Dichloro-1,4-benzoquinone: Insights from Experimental and Theoretical Charge Density Analysis

Experimental charge density distribution in 2, 5-dichloro-1, 4-benzoquinone has been carried out using high resolution X-ray diffraction data at 90 K to quantitatively evaluate the nature of C-Cl center dot center dot center dot O=C halogen bond in molecular crystals. Additionally, the halogen bond is studied from geometrical point of view and the same has been visualized using Hirshfeld surface analysis. The obtained results from experimental charge density analysis are compared with periodic quantum calculations using B3LYP 6-31G(d,p) level of theory. The topological values at bond critical point, three-dimensional static deformation density features and electrostatic potential isosurfaces unequivocally establish the attractive nature of C-Cl center dot center dot center dot O=C halogen bond in crystalline lattice.

Spectral asymptotics of the Helmholtz model in fluid-solid structures

A model representing the vibrations of a coupled fluid-solid structure is considered. This structure consists of a tube bundle immersed in a slightly compressible fluid. Assuming periodic distribution of tubes, this article describes the asymptotic nature of the vibration frequencies when the number of tubes is large. Our investigation shows that classical homogenization of the problem is not sufficient for this purpose. Indeed, our end result proves that the limit spectrum consists of three parts: the macro-part which comes from homogenization, the micro-part and the boundary layer part. The last two components are new. We describe in detail both macro- and micro-parts using the so-called Bloch wave homogenization method. Copyright (C) 1999 John Wiley & Sons, Ltd.

The 2D Coulomb gas on a 1D lattice

The statistical mechanics of a two-dimensional Coulomb gas confined to one dimension is studied, wherein hard core particles move on a ring. Exact self-duality is shown for a version of the sine-Gordon model arising in this context, thereby locating the transition temperature exactly. We present asymptotically exact results for the correlations in the model and characterize the low- and high-temperature phases. Numerical simulations provide support to these renormalization group calculations. Connections with other interesting problems, such as the quantum Brownian motion of a panicle in a periodic potential and impurity problems, are pointed out.

Zero-point entropy in 'spin ice'

Common water ice (ice I-h) is an unusual solid-the oxygen atoms form a periodic structure but the hydrogen atoms are highly disordered due to there being two inequivalent O-H bond lengths'. Pauling showed that the presence of these two bond lengths leads to a macroscopic degeneracy of possible ground states(2,3), such that the system has finite entropy as the temperature tends towards zero. The dynamics associated with this degeneracy are experimentally inaccessible, however, as ice melts and the hydrogen dynamics cannot be studied independently of oxygen motion(4). An analogous system(5) in which this degeneracy can be studied is a magnet with the pyrochlore structure-termed 'spin ice'-where spin orientation plays a similar role to that of the hydrogen position in ice I-h. Here we present specific-heat data for one such system, Dy2Ti2O7, from which we infer a total spin entropy of 0.67Rln2. This is similar to the value, 0.71Rln2, determined for ice I-h, SO confirming the validity of the correspondence. We also find, through application of a magnetic field, behaviour not accessible in water ice-restoration of much of the ground-state entropy and new transitions involving transverse spin degrees of freedom.

Driven Heisenberg magnets: Nonequilibrium criticality, spatiotemporal chaos and control

We drive a d-dimensional Heisenberg magnet using an anisotropic current. The continuum Langevin equation is analysed using a dynamical renormalization group and numerical simulations. We discover a rich steady-state phase diagram, including a critical point in a new nonequilibrium universality class, and a spatiotemporally chaotic phase. The latter may be controlled in a robust manner to target spatially periodic steady states with helical order.

Minimizing buffer requirements under rate-optimal schedule in regular dataflow networks

Large-grain synchronous dataflow graphs or multi-rate graphs have the distinct feature that the nodes of the dataflow graph fire at different rates. Such multi-rate large-grain dataflow graphs have been widely regarded as a powerful programming model for DSP applications. In this paper we propose a method to minimize buffer storage requirement in constructing rate-optimal compile-time (MBRO) schedules for multi-rate dataflow graphs. We demonstrate that the constraints to minimize buffer storage while executing at the optimal computation rate (i.e. the maximum possible computation rate without storage constraints) can be formulated as a unified linear programming problem in our framework. A novel feature of our method is that in constructing the rate-optimal schedule, it directly minimizes the memory requirement by choosing the schedule time of nodes appropriately. Lastly, a new circular-arc interval graph coloring algorithm has been proposed to further reduce the memory requirement by allowing buffer sharing among the arcs of the multi-rate dataflow graph. We have constructed an experimental testbed which implements our MBRO scheduling algorithm as well as (i) the widely used periodic admissible parallel schedules (also known as block schedules) proposed by Lee and Messerschmitt (IEEE Transactions on Computers, vol. 36, no. 1, 1987, pp. 24-35), (ii) the optimal scheduling buffer allocation (OSBA) algorithm of Ning and Gao (Conference Record of the Twentieth Annual ACM SIGPLAN-SIGACT Symposium on Principles of Programming Languages, Charleston, SC, Jan. 10-13, 1993, pp. 29-42), and (iii) the multi-rate software pipelining (MRSP) algorithm (Govindarajan and Gao, in Proceedings of the 1993 International Conference on Application Specific Array Processors, Venice, Italy, Oct. 25-27, 1993, pp. 77-88). Schedules generated for a number of random dataflow graphs and for a set of DSP application programs using the different scheduling methods are compared. The experimental results have demonstrated a significant improvement (10-20%) in buffer requirements for the MBRO schedules compared to the schedules generated by the other three methods, without sacrificing the computation rate. The MBRO method also gives a 20% average improvement in computation rate compared to Lee's Block scheduling method.

Persistent passive hopping and juggling is possible even with plastic collisions

We describe simple one-dimensional models of passive (no energy input, no control), generally dissipative, vertical hopping and one-ball juggling. The central observation is that internal passive system motions can conspire to eliminate collisions in these systems. For hopping, two point masses are connected by a spring and the lower mass has inelastic collisions with the ground. For juggling, a lower point-mass hand is connected by a spring to the ground and an upper point-mass ball is caught with an inelastic collision and then re-thrown into gravitational free flight. The two systems have identical dynamics. Despite inelastic collisions between non-zero masses, these systems have special symmetric energy-conserving periodic motions where the collision is at zero relative velocity. Additionally, these special periodic motions have a non-zero sized, one-sided region of attraction on the higher-energy side. For either very large or very small mass ratios, the one-sided region of attraction is large. These results persist for mildly non-linear springs and non-constant gravity. Although non-collisional damping destroys the periodic motions, small energy injection makes the periodic motions stable, with a two-sided region of attraction. The existence of such special energy conserving solutions for hopping and juggling points to possibly useful strategies for both animals and robots. The lossless motions are demonstrated with a table-top experiment.