Developing a growth-transition matrix for the stock assessment of the green sea urchin (Strongylocentrotus droebachiensis) off Maine

The green sea urchin (Strongylocentrotus droebachiensis) is important to the economy of Maine. It is the state’s fourth largest fishery by value. The fishery has experienced a continuous decline in landings since 1992 because of decreasing stock abundance. Because determining the age of sea urchins is often difficult, a formal stock assessment demands the development of a size-structured population dynamic model. One of the most important components in a size-structured model is a growth-transition matrix. We developed an approach for estimating the growth-transition matrix using von Bertalanffy growth parameters estimated in previous studies of the green sea urchin off Maine. This approach explicitly considers size-specific variations associated with yearly growth increments for these urchins. The proposed growth-transition matrix can be updated readily with new information on growth, which is important because changes in stock abundance and the ecosystem will likely result in changes in sea urchin key life history parameters including growth. This growth-transition matrix can be readily incorporated into the size-structured stock assessment model that has been developed for assessing the green sea urchin stock off Maine.

Comparison between Two Methods to Calculate the Transition Matrix of Orbit Motion

Two methods to evaluate the state transition matrix are implemented and analyzed to verify the computational cost and the accuracy of both methods. This evaluation represents one of the highest computational costs on the artificial satellite orbit determination task. The first method is an approximation of the Keplerian motion, providing an analytical solution which is then calculated numerically by solving Kepler's equation. The second one is a local numerical approximation that includes the effect of J(2). The analysis is performed comparing these two methods with a reference generated by a numerical integrator. For small intervals of time (1 to 10s) and when one needs more accuracy, it is recommended to use the second method, since the CPU time does not excessively overload the computer during the orbit determination procedure. For larger intervals of time and when one expects more stability on the calculation, it is recommended to use the first method.

Inverse Sensitivity Analysis of Singular Solutions of FRF matrix in Structural System Identification

The problem of structural damage detection based on measured frequency response functions of the structure in its damaged and undamaged states is considered. A novel procedure that is based on inverse sensitivity of the singular solutions of the system FRF matrix is proposed. The treatment of possibly ill-conditioned set of equations via regularization scheme and questions on spatial incompleteness of measurements are considered. The application of the method in dealing with systems with repeated natural frequencies and (or) packets of closely spaced modes is demonstrated. The relationship between the proposed method and the methods based on inverse sensitivity of eigensolutions and frequency response functions is noted. The numerical examples on a 5-degree of freedom system, a one span free-free beam and a spatially periodic multi-span beam demonstrate the efficacy of the proposed method and its superior performance vis-a-vis methods based on inverse eigensensitivity.

Full characterization of the quantum linear-zigzag transition in atomic chains

A string of repulsively interacting particles exhibits a phase transition to a zigzag structure, by reducing the transverse trap potential or the interparticle distance. Based on the emergent symmetry Z2 it has been argued that this instability is a quantum phase transition, which can be mapped to an Ising model in transverse field. An extensive Density Matrix Renormalization Group analysis is performed, resulting in an high-precision evaluation of the critical exponents and of the central charge of the system, confirming that the quantum linear-zigzag transition belongs to the critical Ising model universality class. Quantum corrections to the classical phase diagram are computed, and the range of experimental parameters where quantum effects play a role is provided. These results show that structural instabilities of one-dimensional interacting atomic arrays can simulate quantum critical phenomena typical of ferromagnetic systems.

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.

Electronic structure and optical gain saturation of InAs1-xNx/GaAs quantum dots

The electronic band structures and optical gains of InAs1-xNx/GaAs pyramid quantum dots (QDs) are calculated using the ten-band k . p model and the valence force field method. The optical gains are calculated using the zero-dimensional optical gain formula with taking into consideration of both homogeneous and inhomogeneous broadenings due to the size fluctuation of quantum dots which follows a normal distribution. With the variation of QD sizes and nitrogen composition, it can be shown that the nitrogen composition and the strains can significantly affect the energy levels especially the conduction band which has repulsion interaction with nitrogen resonant state due to the band anticrossing interaction. It facilitates to achieve emission of longer wavelength (1.33 or 1.55 mu m) lasers for optical fiber communication system. For QD with higher nitrogen composition, it has longer emission wavelength and less detrimental effect of higher excited state transition, but nitrogen composition can affect the maximum gain depending on the factors of transition matrix element and the Fermi-Dirac distributions for electrons in the conduction bands and holes in the valence bands respectively. For larger QD, its maximum optical gain is greater at lower carrier density, but it is slowly surpassed by smaller QD as carrier concentration increases. Larger QD can reach its saturation gain faster, but this saturation gain is smaller than that of smaller QD. So the trade-off between longer wavelength, maximum optical, saturation gain, and differential gain must be considered to select the appropriate QD size according to the specific application requirement. (C) 2009 American Institute of Physics. [DOI 10.1063/1.3143025]

Electronic structure and optical gain saturation of InAs1-xNx/GaAs quantum dots

The electronic band structures and optical gains of InAs1-xNx/GaAs pyramid quantum dots (QDs) are calculated using the ten-band k . p model and the valence force field method. The optical gains are calculated using the zero-dimensional optical gain formula with taking into consideration of both homogeneous and inhomogeneous broadenings due to the size fluctuation of quantum dots which follows a normal distribution. With the variation of QD sizes and nitrogen composition, it can be shown that the nitrogen composition and the strains can significantly affect the energy levels especially the conduction band which has repulsion interaction with nitrogen resonant state due to the band anticrossing interaction. It facilitates to achieve emission of longer wavelength (1.33 or 1.55 mu m) lasers for optical fiber communication system. For QD with higher nitrogen composition, it has longer emission wavelength and less detrimental effect of higher excited state transition, but nitrogen composition can affect the maximum gain depending on the factors of transition matrix element and the Fermi-Dirac distributions for electrons in the conduction bands and holes in the valence bands respectively. For larger QD, its maximum optical gain is greater at lower carrier density, but it is slowly surpassed by smaller QD as carrier concentration increases. Larger QD can reach its saturation gain faster, but this saturation gain is smaller than that of smaller QD. So the trade-off between longer wavelength, maximum optical, saturation gain, and differential gain must be considered to select the appropriate QD size according to the specific application requirement. (C) 2009 American Institute of Physics. [DOI: 10.1063/1.3143025]

Custom coprocessor based matrix algorithms for image and signal processing

Matrix algorithms are important in many types of applications including image and signal processing. A close examination of the algorithms used in these, and related, applications reveals that many of the fundamental actions involve matrix algorithms such as matrix multiplication. This paper presents an investigation into the design and implementation of different matrix algorithms such as matrix operations, matrix transforms and matrix decompositions using a novel custom coprocessor system for MATrix algorithms based on Reconfigurable Computing (RCMAT). The proposed RCMAT architectures are scalable, modular and require less area and time complexity with reduced latency when compared with existing structures.

Discrete fractional order system vibrations

A theory of free vibrations of discrete fractional order (FO) systems with a finite number of degrees of freedom (dof) is developed. A FO system with a finite number of dof is defined by means of three matrices: mass inertia, system rigidity and FO elements. By adopting a matrix formulation, a mathematical description of FO discrete system free vibrations is determined in the form of coupled fractional order differential equations (FODE). The corresponding solutions in analytical form, for the special case of the matrix of FO properties elements, are determined and expressed as a polynomial series along time. For the eigen characteristic numbers, the system eigen main coordinates and the independent eigen FO modes are determined. A generalized function of visoelastic creep FO dissipation of energy and generalized forces of system with no ideal visoelastic creep FO dissipation of energy for generalized coordinates are formulated. Extended Lagrange FODE of second kind, for FO system dynamics, are also introduced. Two examples of FO chain systems are analyzed and the corresponding eigen characteristic numbers determined. It is shown that the oscillatory phenomena of a FO mechanical chain have analogies to electrical FO circuits. A FO electrical resistor is introduced and its constitutive voltage–current is formulated. Also a function of thermal energy FO dissipation of a FO electrical relation is discussed.

Dynamic game soundtrack generation in response to a continously varying emotional trajectory

Dynamic soundtracking presents various practical and aesthetic challenges to composers working with games. This paper presents an implementation of a system addressing some of these challenges with an affectively-driven music generation algorithm based on a second order Markov-model. The system can respond in real-time to emotional trajectories derived from 2-dimensions of affect on the circumplex model (arousal and valence), which are mapped to five musical parameters. A transition matrix is employed to vary the generated output in continuous response to the affective state intended by the gameplay.

Automated two-dimensional separation flow system with electrochemical preconcentration, stripping, capillary electrophoresis and contactless conductivity detection for trace metal ion analysis

This paper describes the automation of a fully electrochemical system for preconcentration, cleanup, separation and detection, comprising the hyphenation of a thin layer electrochemical flow cell with CE coupled with contactless conductivity detection (CE-C(4)D). Traces of heavy metal ions were extracted from the pulsed-flowing sample and accumulated on a glassy carbon working electrode by electroreduction for some minutes. Anodic stripping of the accumulated metals was synchronized with hydrodynamic injection into the capillary. The effect of the angle of the slant polished tip of the CE capillary and its orientation against the working electrode in the electrochemical preconcentration (EPC) flow cell and of the accumulation time were studied, aiming at maximum CE-C(4)D signal enhancement. After 6 min of EPC, enhancement factors close to 50 times were obtained for thallium, lead, cadmium and copper ions, and about 16 for zinc ions. Limits of detection below 25 nmol/L were estimated for all target analytes but zinc. A second separation dimension was added to the CE separation capabilities by staircase scanning of the potentiostatic deposition and/or stripping potentials of metal ions, as implemented with the EPC-CE-C(4)D flow system. A matrix exchange between the deposition and stripping steps, highly valuable for sample cleanup, can be straightforwardly programmed with the multi-pumping flow management system. The automated simultaneous determination of the traces of five accumulable heavy metals together with four non-accumulated alkaline and alkaline earth metals in a single run was demonstrated, to highlight the potentiality of the system.

Análise da estabilidade de sistemas dinâmicos periódicos usando Teoria de Sinha

Pós-graduação em Matemática - IBILCE

Extracellular matrix composition in COPD

Extracellular matrix (ECM) composition has an important role in determining airway structure. We postulated that ECM lung composition of chronic obstructive pulmonary disease (COPD) patients differs from that observed in smoking and nonsmoking subjects without airflow obstruction. We determined the fractional areas of elastic fibres, type-I, -III and -IV collagen, versican, decorin, biglycan, lumican, fibronectin and tenascin in different compartments of the large and small airways and lung parenchyma in 26 COPD patients, 26 smokers without COPD and 16 nonsmoking control subjects. The fractional area of elastic fibres was higher in non-obstructed smokers than in COPD and nonsmoking controls, in all lung compartments. Type-I collagen fractional area was lower in the large and small airways of COPD patients and in the small airways of non-obstructed smokers than in nonsmokers. Compared with nonsmokers, COPD patients had lower versican fractional area in the parenchyma, higher fibronectin fractional area in small airways and higher tenascin fractional area in large and small airways compartments. In COPD patients, significant correlations were found between elastic fibres and fibronectin and lung function parameters. Alterations of the major ECM components are widespread in all lung compartments of patients with COPD and may contribute to persistent airflow obstruction.

Multimodality and flexibility of stochastic gene expression

We consider a general class of mathematical models for stochastic gene expression where the transcription rate is allowed to depend on a promoter state variable that can take an arbitrary (finite) number of values. We provide the solution of the master equations in the stationary limit, based on a factorization of the stochastic transition matrix that separates timescales and relative interaction strengths, and we express its entries in terms of parameters that have a natural physical and/or biological interpretation. The solution illustrates the capacity of multiple states promoters to generate multimodal distributions of gene products, without the need for feedback. Furthermore, using the example of a three states promoter operating at low, high, and intermediate expression levels, we show that using multiple states operons will typically lead to a significant reduction of noise in the system. The underlying mechanism is that a three-states promoter can change its level of expression from low to high by passing through an intermediate state with a much smaller increase of fluctuations than by means of a direct transition.