2012-01-01
A natural bond orbital (NBO) analysis of unpaired electron spin density in metalloproteins is presented, which allows a fast and robust calculation of paramagnetic NMR parameters. Approximately 90% of the unpaired electron spin density occupies metal–ligand NBOs, allowing the majority of the density to be modeled by only a few NBOs that reflect the chemical bonding environment. We show that the paramagnetic relaxation rate of protons can be calculated accurately using only the metal–ligand NBOs and that these rates are in good agreement with corresponding rates measured experimentally. This holds, in particular, for protons of ligand residues where the point-dipole approximation breaks down. To describe the paramagnetic relaxation of heavy nuclei, also the electron spin density in the local orbitals must be taken into account. Geometric distance restraints for 15N can be derived from the paramagnetic relaxation enhancement and the Fermi contact shift when local NBOs are included in the analysis. Thus, the NBO approach allows us to include experimental paramagnetic NMR parameters of 15N nuclei as restraints in a structure optimization protocol. We performed a molecular dynamics simulation and structure determination of oxidized rubredoxin using the experimentally obtained paramagnetic NMR parameters of 15N. The corresponding structures obtained are in good agreement with the crystal structure of rubredoxin. Thus, the NBO approach allows an accurate description of the geometric structure and the dynamics of metalloproteins, when NMR parameters are available of nuclei in the immediate vicinity of the metal-site. PMID:22329704
Accurate parameter estimation for unbalanced three-phase system.
Chen, Yuan; So, Hing Cheung
2014-01-01
Smart grid is an intelligent power generation and control console in modern electricity networks, where the unbalanced three-phase power system is the commonly used model. Here, parameter estimation for this system is addressed. After converting the three-phase waveforms into a pair of orthogonal signals via the α β-transformation, the nonlinear least squares (NLS) estimator is developed for accurately finding the frequency, phase, and voltage parameters. The estimator is realized by the Newton-Raphson scheme, whose global convergence is studied in this paper. Computer simulations show that the mean square error performance of NLS method can attain the Cramér-Rao lower bound. Moreover, our proposal provides more accurate frequency estimation when compared with the complex least mean square (CLMS) and augmented CLMS. PMID:25162056
Accurate and robust estimation of camera parameters using RANSAC
NASA Astrophysics Data System (ADS)
Zhou, Fuqiang; Cui, Yi; Wang, Yexin; Liu, Liu; Gao, He
2013-03-01
Camera calibration plays an important role in the field of machine vision applications. The popularly used calibration approach based on 2D planar target sometimes fails to give reliable and accurate results due to the inaccurate or incorrect localization of feature points. To solve this problem, an accurate and robust estimation method for camera parameters based on RANSAC algorithm is proposed to detect the unreliability and provide the corresponding solutions. Through this method, most of the outliers are removed and the calibration errors that are the main factors influencing measurement accuracy are reduced. Both simulative and real experiments have been carried out to evaluate the performance of the proposed method and the results show that the proposed method is robust under large noise condition and quite efficient to improve the calibration accuracy compared with the original state.
Machine learning of parameters for accurate semiempirical quantum chemical calculations
Dral, Pavlo O.; von Lilienfeld, O. Anatole; Thiel, Walter
2015-04-14
We investigate possible improvements in the accuracy of semiempirical quantum chemistry (SQC) methods through the use of machine learning (ML) models for the parameters. For a given class of compounds, ML techniques require sufficiently large training sets to develop ML models that can be used for adapting SQC parameters to reflect changes in molecular composition and geometry. The ML-SQC approach allows the automatic tuning of SQC parameters for individual molecules, thereby improving the accuracy without deteriorating transferability to molecules with molecular descriptors very different from those in the training set. The performance of this approach is demonstrated for the semiempiricalmore » OM2 method using a set of 6095 constitutional isomers C7H10O2, for which accurate ab initio atomization enthalpies are available. The ML-OM2 results show improved average accuracy and a much reduced error range compared with those of standard OM2 results, with mean absolute errors in atomization enthalpies dropping from 6.3 to 1.7 kcal/mol. They are also found to be superior to the results from specific OM2 reparameterizations (rOM2) for the same set of isomers. The ML-SQC approach thus holds promise for fast and reasonably accurate high-throughput screening of materials and molecules.« less
Machine learning of parameters for accurate semiempirical quantum chemical calculations
Dral, Pavlo O.; von Lilienfeld, O. Anatole; Thiel, Walter
2015-04-14
We investigate possible improvements in the accuracy of semiempirical quantum chemistry (SQC) methods through the use of machine learning (ML) models for the parameters. For a given class of compounds, ML techniques require sufficiently large training sets to develop ML models that can be used for adapting SQC parameters to reflect changes in molecular composition and geometry. The ML-SQC approach allows the automatic tuning of SQC parameters for individual molecules, thereby improving the accuracy without deteriorating transferability to molecules with molecular descriptors very different from those in the training set. The performance of this approach is demonstrated for the semiempirical OM2 method using a set of 6095 constitutional isomers C_{7}H_{10}O_{2}, for which accurate ab initio atomization enthalpies are available. The ML-OM2 results show improved average accuracy and a much reduced error range compared with those of standard OM2 results, with mean absolute errors in atomization enthalpies dropping from 6.3 to 1.7 kcal/mol. They are also found to be superior to the results from specific OM2 reparameterizations (rOM2) for the same set of isomers. The ML-SQC approach thus holds promise for fast and reasonably accurate high-throughput screening of materials and molecules.
Accurate and transferable extended Hückel-type tight-binding parameters
NASA Astrophysics Data System (ADS)
Cerdá, J.; Soria, F.
2000-03-01
We show how the simple extended Hückel theory can be easily parametrized in order to yield accurate band structures for bulk materials, while the resulting optimized atomic orbital basis sets present good transferability properties. The number of parameters involved is exceedingly small, typically ten or eleven per structural phase. We apply the method to almost fifty elemental and compound bulk phases.
Direct computation of parameters for accurate polarizable force fields
Verstraelen, Toon Vandenbrande, Steven; Ayers, Paul W.
2014-11-21
We present an improved electronic linear response model to incorporate polarization and charge-transfer effects in polarizable force fields. This model is a generalization of the Atom-Condensed Kohn-Sham Density Functional Theory (DFT), approximated to second order (ACKS2): it can now be defined with any underlying variational theory (next to KS-DFT) and it can include atomic multipoles and off-center basis functions. Parameters in this model are computed efficiently as expectation values of an electronic wavefunction, obviating the need for their calibration, regularization, and manual tuning. In the limit of a complete density and potential basis set in the ACKS2 model, the linear response properties of the underlying theory for a given molecular geometry are reproduced exactly. A numerical validation with a test set of 110 molecules shows that very accurate models can already be obtained with fluctuating charges and dipoles. These features greatly facilitate the development of polarizable force fields.
Accurate 3D quantification of the bronchial parameters in MDCT
NASA Astrophysics Data System (ADS)
Saragaglia, A.; Fetita, C.; Preteux, F.; Brillet, P. Y.; Grenier, P. A.
2005-08-01
The assessment of bronchial reactivity and wall remodeling in asthma plays a crucial role in better understanding such a disease and evaluating therapeutic responses. Today, multi-detector computed tomography (MDCT) makes it possible to perform an accurate estimation of bronchial parameters (lumen and wall areas) by allowing a quantitative analysis in a cross-section plane orthogonal to the bronchus axis. This paper provides the tools for such an analysis by developing a 3D investigation method which relies on 3D reconstruction of bronchial lumen and central axis computation. Cross-section images at bronchial locations interactively selected along the central axis are generated at appropriate spatial resolution. An automated approach is then developed for accurately segmenting the inner and outer bronchi contours on the cross-section images. It combines mathematical morphology operators, such as "connection cost", and energy-controlled propagation in order to overcome the difficulties raised by vessel adjacencies and wall irregularities. The segmentation accuracy was validated with respect to a 3D mathematically-modeled phantom of a pair bronchus-vessel which mimics the characteristics of real data in terms of gray-level distribution, caliber and orientation. When applying the developed quantification approach to such a model with calibers ranging from 3 to 10 mm diameter, the lumen area relative errors varied from 3.7% to 0.15%, while the bronchus area was estimated with a relative error less than 5.1%.
Accurate Prediction of Docked Protein Structure Similarity.
Akbal-Delibas, Bahar; Pomplun, Marc; Haspel, Nurit
2015-09-01
One of the major challenges for protein-protein docking methods is to accurately discriminate nativelike structures. The protein docking community agrees on the existence of a relationship between various favorable intermolecular interactions (e.g. Van der Waals, electrostatic, desolvation forces, etc.) and the similarity of a conformation to its native structure. Different docking algorithms often formulate this relationship as a weighted sum of selected terms and calibrate their weights against specific training data to evaluate and rank candidate structures. However, the exact form of this relationship is unknown and the accuracy of such methods is impaired by the pervasiveness of false positives. Unlike the conventional scoring functions, we propose a novel machine learning approach that not only ranks the candidate structures relative to each other but also indicates how similar each candidate is to the native conformation. We trained the AccuRMSD neural network with an extensive dataset using the back-propagation learning algorithm. Our method achieved predicting RMSDs of unbound docked complexes with 0.4Å error margin. PMID:26335807
Identification of accurate nonlinear rainfall-runoff models with unique parameters
NASA Astrophysics Data System (ADS)
Schoups, G.; Vrugt, J. A.; Fenicia, F.; van de Giesen, N.
2009-04-01
We propose a strategy to identify models with unique parameters that yield accurate streamflow predictions, given a time-series of rainfall inputs. The procedure consists of five general steps. First, an a priori range of model structures is specified based on prior general and site-specific hydrologic knowledge. To this end, we rely on a flexible model code that allows a specification of a wide range of model structures, from simple to complex. Second, using global optimization each model structure is calibrated to a record of rainfall-runoff data, yielding optimal parameter values for each model structure. Third, accuracy of each model structure is determined by estimating model prediction errors using independent validation and statistical theory. Fourth, parameter identifiability of each calibrated model structure is estimated by means of Monte Carlo Markov Chain simulation. Finally, an assessment is made about each model structure in terms of its accuracy of mimicking rainfall-runoff processes (step 3), and the uniqueness of its parameters (step 4). The procedure results in the identification of the most complex and accurate model supported by the data, without causing parameter equifinality. As such, it provides insight into the information content of the data for identifying nonlinear rainfall-runoff models. We illustrate the method using rainfall-runoff data records from several MOPEX basins in the US.
Accurate lattice parameter measurements of stoichiometric uranium dioxide
NASA Astrophysics Data System (ADS)
Leinders, Gregory; Cardinaels, Thomas; Binnemans, Koen; Verwerft, Marc
2015-04-01
The paper presents and discusses lattice parameter analyses of pure, stoichiometric UO2. Attention was paid to prepare stoichiometric samples and to maintain stoichiometry throughout the analyses. The lattice parameter of UO2.000±0.001 was evaluated as being 547.127 ± 0.008 pm at 20 °C, which is substantially higher than many published values for the UO2 lattice constant and has an improved precision by about one order of magnitude. The higher value of the lattice constant is mainly attributed to the avoidance of hyperstoichiometry in the present study and to a minor extent to the use of the currently accepted Cu Kα1 X-ray wavelength value. Many of the early studies used Cu Kα1 wavelength values that differ from the currently accepted value, which also contributed to an underestimation of the true lattice parameter.
Clinically accurate fetal ECG parameters acquired from maternal abdominal sensors
CLIFFORD, Gari; SAMENI, Reza; WARD, Mr. Jay; ROBINSON, Julian; WOLFBERG, Adam J.
2011-01-01
OBJECTIVE To evaluate the accuracy of a novel system for measuring fetal heart rate and ST-segment changes using non-invasive electrodes on the maternal abdomen. STUDY DESIGN Fetal ECGs were recorded using abdominal sensors from 32 term laboring women who had a fetal scalp electrode (FSE) placed for a clinical indication. RESULTS Good quality data for FHR estimation was available in 91.2% of the FSE segments, and 89.9% of the abdominal electrode segments. The root mean square (RMS) error between the FHR data calculated by both methods over all processed segments was 0.36 beats per minute. ST deviation from the isoelectric point ranged from 0 to 14.2% of R-wave amplitude. The RMS error between the ST change calculated by both methods averaged over all processed segments was 3.2%. CONCLUSION FHR and ST change acquired from the maternal abdomen is highly accurate and on average is clinically indistinguishable from FHR and ST change calculated using FSE data. PMID:21514560
Isomerism of Cyanomethanimine: Accurate Structural, Energetic, and Spectroscopic Characterization.
Puzzarini, Cristina
2015-11-25
The structures, relative stabilities, and rotational and vibrational parameters of the Z-C-, E-C-, and N-cyanomethanimine isomers have been evaluated using state-of-the-art quantum-chemical approaches. Equilibrium geometries have been calculated by means of a composite scheme based on coupled-cluster calculations that accounts for the extrapolation to the complete basis set limit and core-correlation effects. The latter approach is proved to provide molecular structures with an accuracy of 0.001-0.002 Å and 0.05-0.1° for bond lengths and angles, respectively. Systematically extrapolated ab initio energies, accounting for electron correlation through coupled-cluster theory, including up to single, double, triple, and quadruple excitations, and corrected for core-electron correlation and anharmonic zero-point vibrational energy, have been used to accurately determine relative energies and the Z-E isomerization barrier with an accuracy of about 1 kJ/mol. Vibrational and rotational spectroscopic parameters have been investigated by means of hybrid schemes that allow us to obtain rotational constants accurate to about a few megahertz and vibrational frequencies with a mean absolute error of ∼1%. Where available, for all properties considered, a very good agreement with experimental data has been observed. PMID:26529434
Distillation tray structural parameter study: Phase 1
NASA Technical Reports Server (NTRS)
Winter, J. Ronald
1991-01-01
The purpose here is to identify the structural parameters (plate thickness, liquid level, beam size, number of beams, tray diameter, etc.) that affect the structural integrity of distillation trays in distillation columns. Once the sensitivity of the trays' dynamic response to these parameters has been established, the designer will be able to use this information to prepare more accurate specifications for the construction of new trays. Information is given on both static and dynamic analysis, modal response, and tray failure details.
Predicting accurate line shape parameters for CO2 transitions
NASA Astrophysics Data System (ADS)
Gamache, Robert R.; Lamouroux, Julien
2013-11-01
The vibrational dependence of CO2 half-widths and line shifts are given by a modification of the model proposed by Gamache and Hartmann [Gamache R, Hartmann J-M. J Quant Spectrosc Radiat Transfer 2004;83:119]. This model allows the half-widths and line shifts for a ro-vibrational transition to be expressed in terms of the number of vibrational quanta exchanged in the transition raised to a power and a reference ro-vibrational transition. Calculations were made for 24 bands for lower rotational quantum numbers from 0 to 160 for N2-, O2-, air-, and self-collisions with CO2. These data were extrapolated to J″=200 to accommodate several databases. Comparison of the CRB calculations with measurement gives very high confidence in the data. In the model a Quantum Coordinate is defined by (c1 |Δν1|+c2 |Δν2|+c3|Δν3|)p. The power p is adjusted and a linear least-squares fit to the data by the model expression is made. The procedure is iterated on the correlation coefficient, R, until [|R|-1] is less than a threshold. The results demonstrate the appropriateness of the model. The model allows the determination of the slope and intercept as a function of rotational transition, broadening gas, and temperature. From the data of the fits, the half-width, line shift, and the temperature dependence of the half-width can be estimated for any ro-vibrational transition, allowing spectroscopic CO2 databases to have complete information for the line shape parameters.
Loewe, Axel; Wilhelms, Mathias; Schmid, Jochen; Krause, Mathias J.; Fischer, Fathima; Thomas, Dierk; Scholz, Eberhard P.; Dössel, Olaf; Seemann, Gunnar
2016-01-01
Computational models of cardiac electrophysiology provided insights into arrhythmogenesis and paved the way toward tailored therapies in the last years. To fully leverage in silico models in future research, these models need to be adapted to reflect pathologies, genetic alterations, or pharmacological effects, however. A common approach is to leave the structure of established models unaltered and estimate the values of a set of parameters. Today’s high-throughput patch clamp data acquisition methods require robust, unsupervised algorithms that estimate parameters both accurately and reliably. In this work, two classes of optimization approaches are evaluated: gradient-based trust-region-reflective and derivative-free particle swarm algorithms. Using synthetic input data and different ion current formulations from the Courtemanche et al. electrophysiological model of human atrial myocytes, we show that neither of the two schemes alone succeeds to meet all requirements. Sequential combination of the two algorithms did improve the performance to some extent but not satisfactorily. Thus, we propose a novel hybrid approach coupling the two algorithms in each iteration. This hybrid approach yielded very accurate estimates with minimal dependency on the initial guess using synthetic input data for which a ground truth parameter set exists. When applied to measured data, the hybrid approach yielded the best fit, again with minimal variation. Using the proposed algorithm, a single run is sufficient to estimate the parameters. The degree of superiority over the other investigated algorithms in terms of accuracy and robustness depended on the type of current. In contrast to the non-hybrid approaches, the proposed method proved to be optimal for data of arbitrary signal to noise ratio. The hybrid algorithm proposed in this work provides an important tool to integrate experimental data into computational models both accurately and robustly allowing to assess the often non
Reinbolt, Jeffrey A.; Haftka, Raphael T.; Chmielewski, Terese L.; Fregly, Benjamin J.
2013-01-01
Variations in joint parameter values (axis positions and orientations in body segments) and inertial parameter values (segment masses, mass centers, and moments of inertia) as well as kinematic noise alter the results of inverse dynamics analyses of gait. Three-dimensional linkage models with joint constraints have been proposed as one way to minimize the effects of noisy kinematic data. Such models can also be used to perform gait optimizations to predict post-treatment function given pre-treatment gait data. This study evaluates whether accurate patient-specific joint and inertial parameter values are needed in three-dimensional linkage models to produce accurate inverse dynamics results for gait. The study was performed in two stages. First, we used optimization analyses to evaluate whether patient-specific joint and inertial parameter values can be calibrated accurately from noisy kinematic data, and second, we used Monte Carlo analyses to evaluate how errors in joint and inertial parameter values affect inverse dynamics calculations. Both stages were performed using a dynamic, 27 degree-of-freedom, full-body linkage model and synthetic (i.e., computer generated) gait data corresponding to a nominal experimental gait motion. In general, joint but not inertial parameter values could be found accurately from noisy kinematic data. Root-mean-square (RMS) errors were 3° and 4 mm for joint parameter values and 1 kg, 22 mm, and 74,500 kg*mm2 for inertial parameter values. Furthermore, errors in joint but not inertial parameter values had a significant effect on calculated lower-extremity inverse dynamics joint torques. The worst RMS torque error averaged 4% bodyweight*height (BW*H) due to joint parameter variations but less than 0.25% BW*H due to inertial parameter variations. These results suggest that inverse dynamics analyses of gait utilizing linkage models with joint constraints should calibrate the model’s joint parameter values to obtain accurate joint
Structure parameters in molecular tunneling ionization theory
NASA Astrophysics Data System (ADS)
Wang, Jun-Ping; Li, Wei; Zhao, Song-Feng
2014-04-01
We extracted the accurate structure parameters in molecular tunneling ionization theory (so called MO-ADK theory) for 22 selected linear molecules including some inner orbitals. The molecular wave functions with the correct asymptotic behavior are obtained by solving the time-independent Schrödinger equation with B-spline functions and molecular potentials numerically constructed using the modified Leeuwen-Baerends (LBα) model.
A simple and accurate resist parameter extraction method for sub-80-nm DRAM patterns
NASA Astrophysics Data System (ADS)
Lee, Sook; Hwang, Chan; Park, Dong-Woon; Kim, In-Sung; Kim, Ho-Chul; Woo, Sang-Gyun; Cho, Han-Ku; Moon, Joo-Tae
2004-05-01
Due to the polarization effect of high NA lithography, the consideration of resist effect in lithography simulation becomes increasingly important. In spite of the importance of resist simulation, many process engineers are reluctant to consider resist effect in lithography simulation due to time-consuming procedure to extract required resist parameters and the uncertainty of measurement of some parameters. Weiss suggested simplified development model, and this model does not require the complex kinetic parameters. For the device fabrication engineers, there is a simple and accurate parameter extraction and optimizing method using Weiss model. This method needs refractive index, Dill"s parameters and development rate monitoring (DRM) data in parameter extraction. The parameters extracted using referred sequence is not accurate, so that we have to optimize the parameters to fit the critical dimension scanning electron microscopy (CD SEM) data of line and space patterns. Hence, the FiRM of Sigma-C is utilized as a resist parameter-optimizing program. According to our study, the illumination shape, the aberration and the pupil mesh point have a large effect on the accuracy of resist parameter in optimization. To obtain the optimum parameters, we need to find the saturated mesh points in terms of normalized intensity log slope (NILS) prior to an optimization. The simulation results using the optimized parameters by this method shows good agreement with experiments for iso-dense bias, Focus-Exposure Matrix data and sub 80nm device pattern simulation.
Measuring accurate body parameters of dressed humans with large-scale motion using a Kinect sensor.
Xu, Huanghao; Yu, Yao; Zhou, Yu; Li, Yang; Du, Sidan
2013-01-01
Non-contact human body measurement plays an important role in surveillance, physical healthcare, on-line business and virtual fitting. Current methods for measuring the human body without physical contact usually cannot handle humans wearing clothes, which limits their applicability in public environments. In this paper, we propose an effective solution that can measure accurate parameters of the human body with large-scale motion from a Kinect sensor, assuming that the people are wearing clothes. Because motion can drive clothes attached to the human body loosely or tightly, we adopt a space-time analysis to mine the information across the posture variations. Using this information, we recover the human body, regardless of the effect of clothes, and measure the human body parameters accurately. Experimental results show that our system can perform more accurate parameter estimation on the human body than state-of-the-art methods. PMID:24064597
Accurate equilibrium structures of fluoro- and chloroderivatives of methane
NASA Astrophysics Data System (ADS)
Vogt, Natalja; Demaison, Jean; Rudolph, Heinz Dieter
2014-11-01
This work is a systematic study of molecular structure of fluoro-, chloro-, and fluorochloromethanes. For the first time, the accurate ab initio structure is computed for 10 molecules (CF4, CClF3, CCl2F2, CCl3F, CHClF2, CHCl2F, CH2F2, CH2ClF, CH2Cl2, and CCl4) at the coupled cluster level of electronic structure theory including single and double excitations augmented by a perturbational estimate of the effects of connected triple excitations [CCSD(T)] with all electrons being correlated and Gaussian basis sets of at least quadruple-ζ quality. Furthermore, when possible, namely for the molecules CH2F2, CH2Cl2, CH2ClF, CHClF2, and CCl2F2, accurate semi-experimental equilibrium (rSEe) structure has also been determined. This is achieved through a least-squares structural refinement procedure based on the equilibrium rotational constants of all available isotopomers, determined by correcting the experimental ground-state rotational constants with computed ab initio vibration-rotation interaction constants and electronic g-factors. The computed and semi-experimental equilibrium structures are in excellent agreement with each other, but the rSEe structure is generally more accurate, in particular for the CF and CCl bond lengths. The carbon-halogen bond length is discussed within the framework of the ligand close-packing model as a function of the atomic charges. For this purpose, the accurate equilibrium structures of some other molecules with alternative ligands, such as CH3Li, CF3CCH, and CF3CN, are also computed.
Accurate and efficient reconstruction of deep phylogenies from structured RNAs
Stocsits, Roman R.; Letsch, Harald; Hertel, Jana; Misof, Bernhard; Stadler, Peter F.
2009-01-01
Ribosomal RNA (rRNA) genes are probably the most frequently used data source in phylogenetic reconstruction. Individual columns of rRNA alignments are not independent as a consequence of their highly conserved secondary structures. Unless explicitly taken into account, these correlation can distort the phylogenetic signal and/or lead to gross overestimates of tree stability. Maximum likelihood and Bayesian approaches are of course amenable to using RNA-specific substitution models that treat conserved base pairs appropriately, but require accurate secondary structure models as input. So far, however, no accurate and easy-to-use tool has been available for computing structure-aware alignments and consensus structures that can deal with the large rRNAs. The RNAsalsa approach is designed to fill this gap. Capitalizing on the improved accuracy of pairwise consensus structures and informed by a priori knowledge of group-specific structural constraints, the tool provides both alignments and consensus structures that are of sufficient accuracy for routine phylogenetic analysis based on RNA-specific substitution models. The power of the approach is demonstrated using two rRNA data sets: a mitochondrial rRNA set of 26 Mammalia, and a collection of 28S nuclear rRNAs representative of the five major echinoderm groups. PMID:19723687
Accurate and efficient reconstruction of deep phylogenies from structured RNAs.
Stocsits, Roman R; Letsch, Harald; Hertel, Jana; Misof, Bernhard; Stadler, Peter F
2009-10-01
Ribosomal RNA (rRNA) genes are probably the most frequently used data source in phylogenetic reconstruction. Individual columns of rRNA alignments are not independent as a consequence of their highly conserved secondary structures. Unless explicitly taken into account, these correlation can distort the phylogenetic signal and/or lead to gross overestimates of tree stability. Maximum likelihood and Bayesian approaches are of course amenable to using RNA-specific substitution models that treat conserved base pairs appropriately, but require accurate secondary structure models as input. So far, however, no accurate and easy-to-use tool has been available for computing structure-aware alignments and consensus structures that can deal with the large rRNAs. The RNAsalsa approach is designed to fill this gap. Capitalizing on the improved accuracy of pairwise consensus structures and informed by a priori knowledge of group-specific structural constraints, the tool provides both alignments and consensus structures that are of sufficient accuracy for routine phylogenetic analysis based on RNA-specific substitution models. The power of the approach is demonstrated using two rRNA data sets: a mitochondrial rRNA set of 26 Mammalia, and a collection of 28S nuclear rRNAs representative of the five major echinoderm groups. PMID:19723687
NASA Astrophysics Data System (ADS)
Peng, Liang-You; Gong, Qihuang
2010-12-01
The accurate computations of hydrogenic continuum wave functions are very important in many branches of physics such as electron-atom collisions, cold atom physics, and atomic ionization in strong laser fields, etc. Although there already exist various algorithms and codes, most of them are only reliable in a certain ranges of parameters. In some practical applications, accurate continuum wave functions need to be calculated at extremely low energies, large radial distances and/or large angular momentum number. Here we provide such a code, which can generate accurate hydrogenic continuum wave functions and corresponding Coulomb phase shifts at a wide range of parameters. Without any essential restrict to angular momentum number, the present code is able to give reliable results at the electron energy range [10,10] eV for radial distances of [10,10] a.u. We also find the present code is very efficient, which should find numerous applications in many fields such as strong field physics. Program summaryProgram title: HContinuumGautchi Catalogue identifier: AEHD_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEHD_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 1233 No. of bytes in distributed program, including test data, etc.: 7405 Distribution format: tar.gz Programming language: Fortran90 in fixed format Computer: AMD Processors Operating system: Linux RAM: 20 MBytes Classification: 2.7, 4.5 Nature of problem: The accurate computation of atomic continuum wave functions is very important in many research fields such as strong field physics and cold atom physics. Although there have already existed various algorithms and codes, most of them can only be applicable and reliable in a certain range of parameters. We present here an accurate FORTRAN program for
FAST TRACK COMMUNICATION Accurate estimate of α variation and isotope shift parameters in Na and Mg+
NASA Astrophysics Data System (ADS)
Sahoo, B. K.
2010-12-01
We present accurate calculations of fine-structure constant variation coefficients and isotope shifts in Na and Mg+ using the relativistic coupled-cluster method. In our approach, we are able to discover the roles of various correlation effects explicitly to all orders in these calculations. Most of the results, especially for the excited states, are reported for the first time. It is possible to ascertain suitable anchor and probe lines for the studies of possible variation in the fine-structure constant by using the above results in the considered systems.
Najafizadeh, Laleh; Gandjbakhche, Amir H.; Pourrezaei, Kambiz; Daryoush, Afshin
2013-01-01
Abstract. Modeling behavior of broadband (30 to 1000 MHz) frequency modulated near-infrared (NIR) photons through a phantom is the basis for accurate extraction of optical absorption and scattering parameters of biological turbid media. Photon dynamics in a phantom are predicted using both analytical and numerical simulation and are related to the measured insertion loss (IL) and insertion phase (IP) for a given geometry based on phantom optical parameters. Accuracy of the extracted optical parameters using finite element method (FEM) simulation is compared to baseline analytical calculations from the diffusion equation (DE) for homogenous brain phantoms. NIR spectroscopy is performed using custom-designed, broadband, free-space optical transmitter (Tx) and receiver (Rx) modules that are developed for photon migration at wavelengths of 680, 780, and 820 nm. Differential detection between two optical Rx locations separated by 0.3 cm is employed to eliminate systemic artifacts associated with interfaces of the optical Tx and Rx with the phantoms. Optical parameter extraction is achieved for four solid phantom samples using the least-square-error method in MATLAB (for DE) and COMSOL (for FEM) simulation by fitting data to measured results over broadband and narrowband frequency modulation. Confidence in numerical modeling of the photonic behavior using FEM has been established here by comparing the transmission mode’s experimental results with the predictions made by DE and FEM for known commercial solid brain phantoms. PMID:23322361
Accurate estimation of motion blur parameters in noisy remote sensing image
NASA Astrophysics Data System (ADS)
Shi, Xueyan; Wang, Lin; Shao, Xiaopeng; Wang, Huilin; Tao, Zhong
2015-05-01
The relative motion between remote sensing satellite sensor and objects is one of the most common reasons for remote sensing image degradation. It seriously weakens image data interpretation and information extraction. In practice, point spread function (PSF) should be estimated firstly for image restoration. Identifying motion blur direction and length accurately is very crucial for PSF and restoring image with precision. In general, the regular light-and-dark stripes in the spectrum can be employed to obtain the parameters by using Radon transform. However, serious noise existing in actual remote sensing images often causes the stripes unobvious. The parameters would be difficult to calculate and the error of the result relatively big. In this paper, an improved motion blur parameter identification method to noisy remote sensing image is proposed to solve this problem. The spectrum characteristic of noisy remote sensing image is analyzed firstly. An interactive image segmentation method based on graph theory called GrabCut is adopted to effectively extract the edge of the light center in the spectrum. Motion blur direction is estimated by applying Radon transform on the segmentation result. In order to reduce random error, a method based on whole column statistics is used during calculating blur length. Finally, Lucy-Richardson algorithm is applied to restore the remote sensing images of the moon after estimating blur parameters. The experimental results verify the effectiveness and robustness of our algorithm.
BIOFILM IMAGE RECONSTRUCTION FOR ASSESSING STRUCTURAL PARAMETERS
Renslow, Ryan; Lewandowski, Zbigniew; Beyenal, Haluk
2011-01-01
The structure of biofilms can be numerically quantified from microscopy images using structural parameters. These parameters are used in biofilm image analysis to compare biofilms, to monitor temporal variation in biofilm structure, to quantify the effects of antibiotics on biofilm structure and to determine the effects of environmental conditions on biofilm structure. It is often hypothesized that biofilms with similar structural parameter values will have similar structures; however, this hypothesis has never been tested. The main goal was to test the hypothesis that the commonly used structural parameters can characterize the differences or similarities between biofilm structures. To achieve this goal 1) biofilm image reconstruction was developed as a new tool for assessing structural parameters, 2) independent reconstructions using the same starting structural parameters were tested to see how they differed from each other, 3) the effect of the original image parameter values on reconstruction success was evaluated and 4) the effect of the number and type of the parameters on reconstruction success was evaluated. It was found that two biofilms characterized by identical commonly used structural parameter values may look different, that the number and size of clusters in the original biofilm image affect image reconstruction success and that, in general, a small set of arbitrarily selected parameters may not reveal relevant differences between biofilm structures. PMID:21280029
Damon, Bruce M.; Heemskerk, Anneriet M.; Ding, Zhaohua
2012-01-01
Fiber curvature is a functionally significant muscle structural property, but its estimation from diffusion-tensor MRI fiber tracking data may be confounded by noise. The purpose of this study was to investigate the use of polynomial fitting of fiber tracts for improving the accuracy and precision of fiber curvature (κ) measurements. Simulated image datasets were created in order to provide data with known values for κ and pennation angle (θ). Simulations were designed to test the effects of increasing inherent fiber curvature (3.8, 7.9, 11.8, and 15.3 m−1), signal-to-noise ratio (50, 75, 100, and 150), and voxel geometry (13.8 and 27.0 mm3 voxel volume with isotropic resolution; 13.5 mm3 volume with an aspect ratio of 4.0) on κ and θ measurements. In the originally reconstructed tracts, θ was estimated accurately under most curvature and all imaging conditions studied; however, the estimates of κ were imprecise and inaccurate. Fitting the tracts to 2nd order polynomial functions provided accurate and precise estimates of κ for all conditions except very high curvature (κ=15.3 m−1), while preserving the accuracy of the θ estimates. Similarly, polynomial fitting of in vivo fiber tracking data reduced the κ values of fitted tracts from those of unfitted tracts and did not change the θ values. Polynomial fitting of fiber tracts allows accurate estimation of physiologically reasonable values of κ, while preserving the accuracy of θ estimation. PMID:22503094
Wang, Mingyu; Han, Lijuan; Liu, Shasha; Zhao, Xuebing; Yang, Jinghua; Loh, Soh Kheang; Sun, Xiaomin; Zhang, Chenxi; Fang, Xu
2015-09-01
Renewable energy from lignocellulosic biomass has been deemed an alternative to depleting fossil fuels. In order to improve this technology, we aim to develop robust mathematical models for the enzymatic lignocellulose degradation process. By analyzing 96 groups of previously published and newly obtained lignocellulose saccharification results and fitting them to Weibull distribution, we discovered Weibull statistics can accurately predict lignocellulose saccharification data, regardless of the type of substrates, enzymes and saccharification conditions. A mathematical model for enzymatic lignocellulose degradation was subsequently constructed based on Weibull statistics. Further analysis of the mathematical structure of the model and experimental saccharification data showed the significance of the two parameters in this model. In particular, the λ value, defined the characteristic time, represents the overall performance of the saccharification system. This suggestion was further supported by statistical analysis of experimental saccharification data and analysis of the glucose production levels when λ and n values change. In conclusion, the constructed Weibull statistics-based model can accurately predict lignocellulose hydrolysis behavior and we can use the λ parameter to assess the overall performance of enzymatic lignocellulose degradation. Advantages and potential applications of the model and the λ value in saccharification performance assessment were discussed. PMID:26121186
NASA Astrophysics Data System (ADS)
Lachaume, Regis; Rabus, Markus; Jordan, Andres
2015-08-01
In stellar interferometry, the assumption that the observables can be seen as Gaussian, independent variables is the norm. In particular, neither the optical interferometry FITS (OIFITS) format nor the most popular fitting software in the field, LITpro, offer means to specify a covariance matrix or non-Gaussian uncertainties. Interferometric observables are correlated by construct, though. Also, the calibration by an instrumental transfer function ensures that the resulting observables are not Gaussian, even if uncalibrated ones happened to be so.While analytic frameworks have been published in the past, they are cumbersome and there is no generic implementation available. We propose here a relatively simple way of dealing with correlated errors without the need to extend the OIFITS specification or making some Gaussian assumptions. By repeatedly picking at random which interferograms, which calibrator stars, and which are the errors on their diameters, and performing the data processing on the bootstrapped data, we derive a sampling of p(O), the multivariate probability density function (PDF) of the observables O. The results can be stored in a normal OIFITS file. Then, given a model m with parameters P predicting observables O = m(P), we can estimate the PDF of the model parameters f(P) = p(m(P)) by using a density estimation of the observables' PDF p.With observations repeated over different baselines, on nights several days apart, and with a significant set of calibrators systematic errors are de facto taken into account. We apply the technique to a precise and accurate assessment of stellar diameters obtained at the Very Large Telescope Interferometer with PIONIER.
Jeong, Hyunjo; Zhang, Shuzeng; Li, Xiongbing; Barnard, Dan
2015-09-15
The accurate measurement of acoustic nonlinearity parameter β for fluids or solids generally requires making corrections for diffraction effects due to finite size geometry of transmitter and receiver. These effects are well known in linear acoustics, while those for second harmonic waves have not been well addressed and therefore not properly considered in previous studies. In this work, we explicitly define the attenuation and diffraction corrections using the multi-Gaussian beam (MGB) equations which were developed from the quasilinear solutions of the KZK equation. The effects of making these corrections are examined through the simulation of β determination in water. Diffraction corrections are found to have more significant effects than attenuation corrections, and the β values of water can be estimated experimentally with less than 5% errors when the exact second harmonic diffraction corrections are used together with the negligible attenuation correction effects on the basis of linear frequency dependence between attenuation coefficients, α{sub 2} ≃ 2α{sub 1}.
Accurate parameters for HD 209458 and its planet from HST spectrophotometry
NASA Astrophysics Data System (ADS)
del Burgo, C.; Allende Prieto, C.
2016-08-01
We present updated parameters for the star HD 209458 and its transiting giant planet. The stellar angular diameter θ=0.2254±0.0017 mas is obtained from the average ratio between the absolute flux observed with the Hubble Space Telescope and that of the best-fitting Kurucz model atmosphere. This angular diameter represents an improvement in precision of more than four times compared to available interferometric determinations. The stellar radius R⋆=1.20±0.05 R⊙ is ascertained by combining the angular diameter with the Hipparcos trigonometric parallax, which is the main contributor to its uncertainty, and therefore the radius accuracy should be significantly improved with Gaia's measurements. The radius of the exoplanet Rp=1.41±0.06 RJ is derived from the corresponding transit depth in the light curve and our stellar radius. From the model fitting, we accurately determine the effective temperature, Teff=6071±20 K, which is in perfect agreement with the value of 6070±24 K calculated from the angular diameter and the integrated spectral energy distribution. We also find precise values from recent Padova Isochrones, such as R⋆=1.20±0.06 R⊙ and Teff=6099±41 K. We arrive at a consistent picture from these methods and compare the results with those from the literature.
Parameter identification of civil engineering structures
NASA Technical Reports Server (NTRS)
Juang, J. N.; Sun, C. T.
1980-01-01
This paper concerns the development of an identification method required in determining structural parameter variations for systems subjected to an extended exposure to the environment. The concept of structural identifiability of a large scale structural system in the absence of damping is presented. Three criteria are established indicating that a large number of system parameters (the coefficient parameters of the differential equations) can be identified by a few actuators and sensors. An eight-bay-fifteen-story frame structure is used as example. A simple model is employed for analyzing the dynamic response of the frame structure.
Parameter estimation for distributed parameter models of complex, flexible structures
NASA Technical Reports Server (NTRS)
Taylor, Lawrence W., Jr.
1991-01-01
Distributed parameter modeling of structural dynamics has been limited to simple spacecraft configurations because of the difficulty of handling several distributed parameter systems linked at their boundaries. Although there is other computer software able to generate such models or complex, flexible spacecraft, unfortunately, neither is suitable for parameter estimation. Because of this limitation the computer software PDEMOD is being developed for the express purposes of modeling, control system analysis, parameter estimation and structure optimization. PDEMOD is capable of modeling complex, flexible spacecraft which consist of a three-dimensional network of flexible beams and rigid bodies. Each beam has bending (Bernoulli-Euler or Timoshenko) in two directions, torsion, and elongation degrees of freedom. The rigid bodies can be attached to the beam ends at any angle or body location. PDEMOD is also capable of performing parameter estimation based on matching experimental modal frequencies and static deflection test data. The underlying formulation and the results of using this approach for test data of the Mini-MAST truss will be discussed. The resulting accuracy of the parameter estimates when using such limited data can impact significantly the instrumentation requirements for on-orbit tests.
Ralph, Duncan K.; Matsen, Frederick A.
2016-01-01
VDJ rearrangement and somatic hypermutation work together to produce antibody-coding B cell receptor (BCR) sequences for a remarkable diversity of antigens. It is now possible to sequence these BCRs in high throughput; analysis of these sequences is bringing new insight into how antibodies develop, in particular for broadly-neutralizing antibodies against HIV and influenza. A fundamental step in such sequence analysis is to annotate each base as coming from a specific one of the V, D, or J genes, or from an N-addition (a.k.a. non-templated insertion). Previous work has used simple parametric distributions to model transitions from state to state in a hidden Markov model (HMM) of VDJ recombination, and assumed that mutations occur via the same process across sites. However, codon frame and other effects have been observed to violate these parametric assumptions for such coding sequences, suggesting that a non-parametric approach to modeling the recombination process could be useful. In our paper, we find that indeed large modern data sets suggest a model using parameter-rich per-allele categorical distributions for HMM transition probabilities and per-allele-per-position mutation probabilities, and that using such a model for inference leads to significantly improved results. We present an accurate and efficient BCR sequence annotation software package using a novel HMM “factorization” strategy. This package, called partis (https://github.com/psathyrella/partis/), is built on a new general-purpose HMM compiler that can perform efficient inference given a simple text description of an HMM. PMID:26751373
Ralph, Duncan K; Matsen, Frederick A
2016-01-01
VDJ rearrangement and somatic hypermutation work together to produce antibody-coding B cell receptor (BCR) sequences for a remarkable diversity of antigens. It is now possible to sequence these BCRs in high throughput; analysis of these sequences is bringing new insight into how antibodies develop, in particular for broadly-neutralizing antibodies against HIV and influenza. A fundamental step in such sequence analysis is to annotate each base as coming from a specific one of the V, D, or J genes, or from an N-addition (a.k.a. non-templated insertion). Previous work has used simple parametric distributions to model transitions from state to state in a hidden Markov model (HMM) of VDJ recombination, and assumed that mutations occur via the same process across sites. However, codon frame and other effects have been observed to violate these parametric assumptions for such coding sequences, suggesting that a non-parametric approach to modeling the recombination process could be useful. In our paper, we find that indeed large modern data sets suggest a model using parameter-rich per-allele categorical distributions for HMM transition probabilities and per-allele-per-position mutation probabilities, and that using such a model for inference leads to significantly improved results. We present an accurate and efficient BCR sequence annotation software package using a novel HMM "factorization" strategy. This package, called partis (https://github.com/psathyrella/partis/), is built on a new general-purpose HMM compiler that can perform efficient inference given a simple text description of an HMM. PMID:26751373
Effects of Structural Errors on Parameter Estimates
NASA Technical Reports Server (NTRS)
Hadaegh, F. Y.; Bekey, G. A.
1987-01-01
Paper introduces concept of near equivalence in probability between different parameters or mathematical models of physical system. One in series of papers, each establishes different part of rigorous theory of mathematical modeling based on concepts of structural error, identifiability, and equivalence. This installment focuses upon effects of additive structural errors on degree of bias in estimates parameters.
Accurate protein structure modeling using sparse NMR data and homologous structure information
Thompson, James M.; Sgourakis, Nikolaos G.; Liu, Gaohua; Rossi, Paolo; Tang, Yuefeng; Mills, Jeffrey L.; Szyperski, Thomas; Montelione, Gaetano T.; Baker, David
2012-01-01
While information from homologous structures plays a central role in X-ray structure determination by molecular replacement, such information is rarely used in NMR structure determination because it can be incorrect, both locally and globally, when evolutionary relationships are inferred incorrectly or there has been considerable evolutionary structural divergence. Here we describe a method that allows robust modeling of protein structures of up to 225 residues by combining , 13C, and 15N backbone and 13Cβ chemical shift data, distance restraints derived from homologous structures, and a physically realistic all-atom energy function. Accurate models are distinguished from inaccurate models generated using incorrect sequence alignments by requiring that (i) the all-atom energies of models generated using the restraints are lower than models generated in unrestrained calculations and (ii) the low-energy structures converge to within 2.0 Å backbone rmsd over 75% of the protein. Benchmark calculations on known structures and blind targets show that the method can accurately model protein structures, even with very remote homology information, to a backbone rmsd of 1.2–1.9 Å relative to the conventional determined NMR ensembles and of 0.9–1.6 Å relative to X-ray structures for well-defined regions of the protein structures. This approach facilitates the accurate modeling of protein structures using backbone chemical shift data without need for side-chain resonance assignments and extensive analysis of NOESY cross-peak assignments. PMID:22665781
Accurate analytical method for the extraction of solar cell model parameters
NASA Astrophysics Data System (ADS)
Phang, J. C. H.; Chan, D. S. H.; Phillips, J. R.
1984-05-01
Single diode solar cell model parameters are rapidly extracted from experimental data by means of the presently derived analytical expressions. The parameter values obtained have a less than 5 percent error for most solar cells, in light of the extraction of model parameters for two cells of differing quality which were compared with parameters extracted by means of the iterative method.
NASA Astrophysics Data System (ADS)
Ghezzi, Luan; Dutra-Ferreira, Letícia; Lorenzo-Oliveira, Diego; Porto de Mello, Gustavo F.; Santiago, Basílio X.; De Lee, Nathan; Lee, Brian L.; da Costa, Luiz N.; Maia, Marcio A. G.; Ogando, Ricardo L. C.; Wisniewski, John P.; González Hernández, Jonay I.; Stassun, Keivan G.; Fleming, Scott W.; Schneider, Donald P.; Mahadevan, Suvrath; Cargile, Phillip; Ge, Jian; Pepper, Joshua; Wang, Ji; Paegert, Martin
2014-12-01
Studies of Galactic chemical, and dynamical evolution in the solar neighborhood depend on the availability of precise atmospheric parameters (effective temperature T eff, metallicity [Fe/H], and surface gravity log g) for solar-type stars. Many large-scale spectroscopic surveys operate at low to moderate spectral resolution for efficiency in observing large samples, which makes the stellar characterization difficult due to the high degree of blending of spectral features. Therefore, most surveys employ spectral synthesis, which is a powerful technique, but relies heavily on the completeness and accuracy of atomic line databases and can yield possibly correlated atmospheric parameters. In this work, we use an alternative method based on spectral indices to determine the atmospheric parameters of a sample of nearby FGK dwarfs and subgiants observed by the MARVELS survey at moderate resolving power (R ~ 12,000). To avoid a time-consuming manual analysis, we have developed three codes to automatically normalize the observed spectra, measure the equivalent widths of the indices, and, through a comparison of those with values calculated with predetermined calibrations, estimate the atmospheric parameters of the stars. The calibrations were derived using a sample of 309 stars with precise stellar parameters obtained from the analysis of high-resolution FEROS spectra, permitting the low-resolution equivalent widths to be directly related to the stellar parameters. A validation test of the method was conducted with a sample of 30 MARVELS targets that also have reliable atmospheric parameters derived from the high-resolution spectra and spectroscopic analysis based on the excitation and ionization equilibria method. Our approach was able to recover the parameters within 80 K for T eff, 0.05 dex for [Fe/H], and 0.15 dex for log g, values that are lower than or equal to the typical external uncertainties found between different high-resolution analyses. An additional test was
Ghezzi, Luan; Da Costa, Luiz N.; Maia, Marcio A. G.; Ogando, Ricardo L. C.; Dutra-Ferreira, Letícia; Lorenzo-Oliveira, Diego; Porto de Mello, Gustavo F.; Santiago, Basílio X.; De Lee, Nathan; Lee, Brian L.; Ge, Jian; Wisniewski, John P.; González Hernández, Jonay I.; Stassun, Keivan G.; Cargile, Phillip; Pepper, Joshua; Fleming, Scott W.; Schneider, Donald P.; Mahadevan, Suvrath; Wang, Ji; and others
2014-12-01
Studies of Galactic chemical, and dynamical evolution in the solar neighborhood depend on the availability of precise atmospheric parameters (effective temperature T {sub eff}, metallicity [Fe/H], and surface gravity log g) for solar-type stars. Many large-scale spectroscopic surveys operate at low to moderate spectral resolution for efficiency in observing large samples, which makes the stellar characterization difficult due to the high degree of blending of spectral features. Therefore, most surveys employ spectral synthesis, which is a powerful technique, but relies heavily on the completeness and accuracy of atomic line databases and can yield possibly correlated atmospheric parameters. In this work, we use an alternative method based on spectral indices to determine the atmospheric parameters of a sample of nearby FGK dwarfs and subgiants observed by the MARVELS survey at moderate resolving power (R ∼ 12,000). To avoid a time-consuming manual analysis, we have developed three codes to automatically normalize the observed spectra, measure the equivalent widths of the indices, and, through a comparison of those with values calculated with predetermined calibrations, estimate the atmospheric parameters of the stars. The calibrations were derived using a sample of 309 stars with precise stellar parameters obtained from the analysis of high-resolution FEROS spectra, permitting the low-resolution equivalent widths to be directly related to the stellar parameters. A validation test of the method was conducted with a sample of 30 MARVELS targets that also have reliable atmospheric parameters derived from the high-resolution spectra and spectroscopic analysis based on the excitation and ionization equilibria method. Our approach was able to recover the parameters within 80 K for T {sub eff}, 0.05 dex for [Fe/H], and 0.15 dex for log g, values that are lower than or equal to the typical external uncertainties found between different high-resolution analyses. An
Robert, Stéphane; Battie, Yann; Jamon, Damien; Royer, Francois
2007-04-10
Optimal performances of integrated optical devices are obtained by the use of an accurate and reliable characterization method. The parameters of interest, i.e., optical indices and thickness of the waveguide structure, are calculated from effective indices by means of an inversion procedure. We demonstrate how an artificial neural network can achieve such a process. The artificial neural network used is a multilayer perceptron. The first result concerns a simulated anisotropic waveguide. The accuracy in the determination of optical indices and waveguide thickness is 5 x 10(-5) and 4 nm, respectively. Then an experimental application on a silica-titania thin film is performed. In addition, effective indices are measured by m-lines spectroscopy. Finally, a comparison with a classical optimization algorithm demonstrates the robustness of the neural method. PMID:17384718
Accurate nuclear masses from a three parameter Kohn-Sham DFT approach (BCPM)
Baldo, M.; Robledo, L. M.; Schuck, P.; Vinas, X.
2012-10-20
Given the promising features of the recently proposed Barcelona-Catania-Paris (BCP) functional [1], it is the purpose of this work to still improve on it. It is, for instance, shown that the number of open parameters can be reduced from 4-5 to 2-3, i.e. by practically a factor of two without deteriorating the results.
Accurate parameters of the oldest known rocky-exoplanet hosting system: Kepler-10 revisited
Fogtmann-Schulz, Alexandra; Hinrup, Brian; Van Eylen, Vincent; Christensen-Dalsgaard, Jørgen; Kjeldsen, Hans; Silva Aguirre, Víctor; Tingley, Brandon
2014-02-01
Since the discovery of Kepler-10, the system has received considerable interest because it contains a small, rocky planet which orbits the star in less than a day. The system's parameters, announced by the Kepler team and subsequently used in further research, were based on only five months of data. We have reanalyzed this system using the full span of 29 months of Kepler photometric data, and obtained improved information about its star and the planets. A detailed asteroseismic analysis of the extended time series provides a significant improvement on the stellar parameters: not only can we state that Kepler-10 is the oldest known rocky-planet-harboring system at 10.41 ± 1.36 Gyr, but these parameters combined with improved planetary parameters from new transit fits gives us the radius of Kepler-10b to within just 125 km. A new analysis of the full planetary phase curve leads to new estimates on the planetary temperature and albedo, which remain degenerate in the Kepler band. Our modeling suggests that the flux level during the occultation is slightly lower than at the transit wings, which would imply that the nightside of this planet has a non-negligible temperature.
NASA Astrophysics Data System (ADS)
Hochlaf, M.; Puzzarini, C.; Senent, M. L.
2015-07-01
We present multi-component computations for rotational constants, vibrational and torsional levels of medium-sized molecules. Through the treatment of two organic sulphur molecules, ethyl mercaptan and dimethyl sulphide, which are relevant for atmospheric and astrophysical media, we point out the outstanding capabilities of explicitly correlated coupled clusters (CCSD(T)-F12) method in conjunction with the cc-pVTZ-F12 basis set for the accurate predictions of such quantities. Indeed, we show that the CCSD(T)-F12/cc-pVTZ-F12 equilibrium rotational constants are in good agreement with those obtained by means of a composite scheme based on CCSD(T) calculations that accounts for the extrapolation to the complete basis set (CBS) limit and core-correlation effects [CCSD(T)/CBS+CV], thus leading to values of ground-state rotational constants rather close to the corresponding experimental data. For vibrational and torsional levels, our analysis reveals that the anharmonic frequencies derived from CCSD(T)-F12/cc-pVTZ-F12 harmonic frequencies and anharmonic corrections (Δν = ω - ν) at the CCSD/cc-pVTZ level closely agree with experimental results. The pattern of the torsional transitions and the shape of the potential energy surfaces along the torsional modes are also well reproduced using the CCSD(T)-F12/cc-pVTZ-F12 energies. Interestingly, this good accuracy is accompanied with a strong reduction of the computational costs. This makes the procedures proposed here as schemes of choice for effective and accurate prediction of spectroscopic properties of organic compounds. Finally, popular density functional approaches are compared with the coupled cluster (CC) methodologies in torsional studies. The long-range CAM-B3LYP functional of Handy and co-workers is recommended for large systems.
Lower bound on reliability for Weibull distribution when shape parameter is not estimated accurately
NASA Technical Reports Server (NTRS)
Huang, Zhaofeng; Porter, Albert A.
1991-01-01
The mathematical relationships between the shape parameter Beta and estimates of reliability and a life limit lower bound for the two parameter Weibull distribution are investigated. It is shown that under rather general conditions, both the reliability lower bound and the allowable life limit lower bound (often called a tolerance limit) have unique global minimums over a range of Beta. Hence lower bound solutions can be obtained without assuming or estimating Beta. The existence and uniqueness of these lower bounds are proven. Some real data examples are given to show how these lower bounds can be easily established and to demonstrate their practicality. The method developed here has proven to be extremely useful when using the Weibull distribution in analysis of no-failure or few-failures data. The results are applicable not only in the aerospace industry but anywhere that system reliabilities are high.
Lower bound on reliability for Weibull distribution when shape parameter is not estimated accurately
NASA Technical Reports Server (NTRS)
Huang, Zhaofeng; Porter, Albert A.
1990-01-01
The mathematical relationships between the shape parameter Beta and estimates of reliability and a life limit lower bound for the two parameter Weibull distribution are investigated. It is shown that under rather general conditions, both the reliability lower bound and the allowable life limit lower bound (often called a tolerance limit) have unique global minimums over a range of Beta. Hence lower bound solutions can be obtained without assuming or estimating Beta. The existence and uniqueness of these lower bounds are proven. Some real data examples are given to show how these lower bounds can be easily established and to demonstrate their practicality. The method developed here has proven to be extremely useful when using the Weibull distribution in analysis of no-failure or few-failures data. The results are applicable not only in the aerospace industry but anywhere that system reliabilities are high.
Proline puckering parameters for collagen structure simulations
Wu, Di
2015-03-15
Collagen is made of triple helices rich in proline residues, and hence is influenced by the conformational motions of prolines. Because the backbone motions of prolines are restricted by the helical structures, the only side chain motion—proline puckering—becomes an influential factor that may affect the stability of collagen structures. In molecular simulations, a proper proline puckering population is desired so to yield valid results of the collagen properties. Here we design the proline puckering parameters in order to yield suitable proline puckering populations as demonstrated in the experimental results. We test these parameters in collagen and the proline dipeptide simulations. Compared with the results of the PDB and the quantum calculations, we propose the proline puckering parameters for the selected collagen model simulations.
Beyond Ellipse(s): Accurately Modelling the Isophotal Structure of Galaxies with ISOFIT and CMODEL
NASA Astrophysics Data System (ADS)
Ciambur, B. C.
2015-09-01
This work introduces a new fitting formalism for isophotes that enables more accurate modeling of galaxies with non-elliptical shapes, such as disk galaxies viewed edge-on or galaxies with X-shaped/peanut bulges. Within this scheme, the angular parameter that defines quasi-elliptical isophotes is transformed from the commonly used, but inappropriate, polar coordinate to the “eccentric anomaly.” This provides a superior description of deviations from ellipticity, better capturing the true isophotal shape. Furthermore, this makes it possible to accurately recover both the surface brightness profile, using the correct azimuthally averaged isophote, and the two-dimensional model of any galaxy: the hitherto ubiquitous, but artificial, cross-like features in residual images are completely removed. The formalism has been implemented into the Image Reduction and Analysis Facility tasks Ellipse and Bmodel to create the new tasks “Isofit,” and “Cmodel.” The new tools are demonstrated here with application to five galaxies, chosen to be representative case-studies for several areas where this technique makes it possible to gain new scientific insight. Specifically: properly quantifying boxy/disky isophotes via the fourth harmonic order in edge-on galaxies, quantifying X-shaped/peanut bulges, higher-order Fourier moments for modeling bars in disks, and complex isophote shapes. Higher order (n > 4) harmonics now become meaningful and may correlate with structural properties, as boxyness/diskyness is known to do. This work also illustrates how the accurate construction, and subtraction, of a model from a galaxy image facilitates the identification and recovery of over-lapping sources such as globular clusters and the optical counterparts of X-ray sources.
Andrianaki, Maria; Azariadis, Kalliopi; Kampouri, Errika; Theodoropoulou, Katerina; Lavrentaki, Katerina; Kastrinakis, Stelios; Kampa, Marilena; Agouridakis, Panagiotis; Pirintsos, Stergios; Castanas, Elias
2015-01-01
Severe allergic reactions of unknown etiology,necessitating a hospital visit, have an important impact in the life of affected individuals and impose a major economic burden to societies. The prediction of clinically severe allergic reactions would be of great importance, but current attempts have been limited by the lack of a well-founded applicable methodology and the wide spatiotemporal distribution of allergic reactions. The valid prediction of severe allergies (and especially those needing hospital treatment) in a region, could alert health authorities and implicated individuals to take appropriate preemptive measures. In the present report we have collecterd visits for serious allergic reactions of unknown etiology from two major hospitals in the island of Crete, for two distinct time periods (validation and test sets). We have used the Normalized Difference Vegetation Index (NDVI), a satellite-based, freely available measurement, which is an indicator of live green vegetation at a given geographic area, and a set of meteorological data to develop a model capable of describing and predicting severe allergic reaction frequency. Our analysis has retained NDVI and temperature as accurate identifiers and predictors of increased hospital severe allergic reactions visits. Our approach may contribute towards the development of satellite-based modules, for the prediction of severe allergic reactions in specific, well-defined geographical areas. It could also probably be used for the prediction of other environment related diseases and conditions. PMID:25794106
An Accurate and Generic Testing Approach to Vehicle Stability Parameters Based on GPS and INS
Miao, Zhibin; Zhang, Hongtian; Zhang, Jinzhu
2015-01-01
With the development of the vehicle industry, controlling stability has become more and more important. Techniques of evaluating vehicle stability are in high demand. As a common method, usually GPS sensors and INS sensors are applied to measure vehicle stability parameters by fusing data from the two system sensors. Although prior model parameters should be recognized in a Kalman filter, it is usually used to fuse data from multi-sensors. In this paper, a robust, intelligent and precise method to the measurement of vehicle stability is proposed. First, a fuzzy interpolation method is proposed, along with a four-wheel vehicle dynamic model. Second, a two-stage Kalman filter, which fuses the data from GPS and INS, is established. Next, this approach is applied to a case study vehicle to measure yaw rate and sideslip angle. The results show the advantages of the approach. Finally, a simulation and real experiment is made to verify the advantages of this approach. The experimental results showed the merits of this method for measuring vehicle stability, and the approach can meet the design requirements of a vehicle stability controller. PMID:26690154
An Accurate and Generic Testing Approach to Vehicle Stability Parameters Based on GPS and INS.
Miao, Zhibin; Zhang, Hongtian; Zhang, Jinzhu
2015-01-01
With the development of the vehicle industry, controlling stability has become more and more important. Techniques of evaluating vehicle stability are in high demand. As a common method, usually GPS sensors and INS sensors are applied to measure vehicle stability parameters by fusing data from the two system sensors. Although prior model parameters should be recognized in a Kalman filter, it is usually used to fuse data from multi-sensors. In this paper, a robust, intelligent and precise method to the measurement of vehicle stability is proposed. First, a fuzzy interpolation method is proposed, along with a four-wheel vehicle dynamic model. Second, a two-stage Kalman filter, which fuses the data from GPS and INS, is established. Next, this approach is applied to a case study vehicle to measure yaw rate and sideslip angle. The results show the advantages of the approach. Finally, a simulation and real experiment is made to verify the advantages of this approach. The experimental results showed the merits of this method for measuring vehicle stability, and the approach can meet the design requirements of a vehicle stability controller. PMID:26690154
Accurate motion parameter estimation for colonoscopy tracking using a regression method
NASA Astrophysics Data System (ADS)
Liu, Jianfei; Subramanian, Kalpathi R.; Yoo, Terry S.
2010-03-01
Co-located optical and virtual colonoscopy images have the potential to provide important clinical information during routine colonoscopy procedures. In our earlier work, we presented an optical flow based algorithm to compute egomotion from live colonoscopy video, permitting navigation and visualization of the corresponding patient anatomy. In the original algorithm, motion parameters were estimated using the traditional Least Sum of squares(LS) procedure which can be unstable in the context of optical flow vectors with large errors. In the improved algorithm, we use the Least Median of Squares (LMS) method, a robust regression method for motion parameter estimation. Using the LMS method, we iteratively analyze and converge toward the main distribution of the flow vectors, while disregarding outliers. We show through three experiments the improvement in tracking results obtained using the LMS method, in comparison to the LS estimator. The first experiment demonstrates better spatial accuracy in positioning the virtual camera in the sigmoid colon. The second and third experiments demonstrate the robustness of this estimator, resulting in longer tracked sequences: from 300 to 1310 in the ascending colon, and 410 to 1316 in the transverse colon.
Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations
NASA Technical Reports Server (NTRS)
Anderson, W. Kyle; Batina, John T.
1988-01-01
Parameter studies are conducted using the Euler and potential flow equation models for steady and unsteady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux vector splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flow results are made, as well as with experimental data where available.
Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations
NASA Technical Reports Server (NTRS)
Anderson, W. Kyle; Batina, John T.
1988-01-01
Parameter studies are conducted using the Euler and potential flow equation models for unsteady and steady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux-vector-splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flows results are made as well as with experimental data where available.
Cartwright, Michael S; Dupuis, Janae E; Bargoil, Jessica M; Foster, Dana C
2015-09-01
Mild traumatic brain injury, often referred to as concussion, is a common, potentially debilitating, and costly condition. One of the main challenges in diagnosing and managing concussion is that there is not currently an objective test to determine the presence of a concussion and to guide return-to-play decisions for athletes. Traditional neuroimaging tests, such as brain magnetic resonance imaging, are normal in concussion, and therefore diagnosis and management are guided by reported symptoms. Some athletes will under-report symptoms to accelerate their return-to-play and others will over-report symptoms out of fear of further injury or misinterpretation of underlying conditions, such as migraine headache. Therefore, an objective measure is needed to assist in several facets of concussion management. Limited data in animal and human testing indicates that intracranial pressure increases slightly and cerebrovascular reactivity (the ability of the cerebral arteries to auto-regulate in response to changes in carbon dioxide) decreases slightly following mild traumatic brain injury. We hypothesize that a combination of ultrasonographic measurements (optic nerve sheath diameter and transcranial Doppler assessment of cerebrovascular reactivity) into a single index will allow for an accurate and non-invasive measurement of intracranial pressure and cerebrovascular reactivity, and this index will be clinically relevant and useful for guiding concussion diagnosis and management. Ultrasound is an ideal modality for the evaluation of concussion because it is portable (allowing for evaluation in many settings, such as on the playing field or in a combat zone), radiation-free (making repeat scans safe), and relatively inexpensive (resulting in nearly universal availability). This paper reviews the literature supporting our hypothesis that an ultrasonographic index can assist in the diagnosis and management of concussion, and it also presents limited data regarding the
Natural frequencies of structures with interval parameters
NASA Astrophysics Data System (ADS)
Sofi, A.; Muscolino, G.; Elishakoff, I.
2015-07-01
This paper deals with the evaluation of the lower and upper bounds of the natural frequencies of structures with uncertain-but-bounded parameters. The solution of the generalized interval eigenvalue problem is pursued by taking into account the actual variability and dependencies of uncertain structural parameters affecting the mass and stiffness matrices. To this aim, interval uncertainties are handled by applying the improved interval analysis via extra unitary interval (EUI), recently introduced by the first two authors. By associating an EUI to each uncertain-but-bounded parameter, the cases of mass and stiffness matrices affected by fully disjoint, completely or partially coincident uncertainties are considered. Then, based on sensitivity analysis, it is shown that the bounds of the interval eigenvalues can be evaluated as solution of two appropriate deterministic eigenvalue problems without requiring any combinatorial procedure. If the eigenvalues are monotonic functions of the uncertain parameters, then the exact bounds are obtained. The accuracy of the proposed method is demonstrated by numerical results concerning truss and beam structures with material and/or geometrical uncertainties.
Sela, Itamar; Ashkenazy, Haim; Katoh, Kazutaka; Pupko, Tal
2015-01-01
Inference of multiple sequence alignments (MSAs) is a critical part of phylogenetic and comparative genomics studies. However, from the same set of sequences different MSAs are often inferred, depending on the methodologies used and the assumed parameters. Much effort has recently been devoted to improving the ability to identify unreliable alignment regions. Detecting such unreliable regions was previously shown to be important for downstream analyses relying on MSAs, such as the detection of positive selection. Here we developed GUIDANCE2, a new integrative methodology that accounts for: (i) uncertainty in the process of indel formation, (ii) uncertainty in the assumed guide tree and (iii) co-optimal solutions in the pairwise alignments, used as building blocks in progressive alignment algorithms. We compared GUIDANCE2 with seven methodologies to detect unreliable MSA regions using extensive simulations and empirical benchmarks. We show that GUIDANCE2 outperforms all previously developed methodologies. Furthermore, GUIDANCE2 also provides a set of alternative MSAs which can be useful for downstream analyses. The novel algorithm is implemented as a web-server, available at: http://guidance.tau.ac.il. PMID:25883146
NASA Astrophysics Data System (ADS)
Martínez, M. J.; Marco, F. J.; López, J. A.
2009-02-01
The Hipparcos catalog provides a reference frame at optical wavelengths for the new International Celestial Reference System (ICRS). This new reference system was adopted following the resolution agreed at the 23rd IAU General Assembly held in Kyoto in 1997. Differences in the Hipparcos system of proper motions and the previous materialization of the reference frame, the FK5, are expected to be caused only by the combined effects of the motion of the equinox of the FK5 and the precession of the equator and the ecliptic. Several authors have pointed out an inconsistency between the differences in proper motion of the Hipparcos-FK5 and the correction of the precessional values derived from VLBI and lunar laser ranging (LLR) observations. Most of them have claimed that these discrepancies are due to slightly biased proper motions in the FK5 catalog. The different mathematical models that have been employed to explain these errors have not fully accounted for the discrepancies in the correction of the precessional parameters. Our goal here is to offer an explanation for this fact. We propose the use of independent parametric and nonparametric models. The introduction of a nonparametric model, combined with the inner product in the square integrable functions over the unitary sphere, would give us values which do not depend on the possible interdependencies existing in the data set. The evidence shows that zonal studies are needed. This would lead us to introduce a local nonparametric model. All these models will provide independent corrections to the precessional values, which could then be compared in order to study the reliability in each case. Finally, we obtain values for the precession corrections that are very consistent with those that are currently adopted.
Discrete state model and accurate estimation of loop entropy of RNA secondary structures.
Zhang, Jian; Lin, Ming; Chen, Rong; Wang, Wei; Liang, Jie
2008-03-28
Conformational entropy makes important contribution to the stability and folding of RNA molecule, but it is challenging to either measure or compute conformational entropy associated with long loops. We develop optimized discrete k-state models of RNA backbone based on known RNA structures for computing entropy of loops, which are modeled as self-avoiding walks. To estimate entropy of hairpin, bulge, internal loop, and multibranch loop of long length (up to 50), we develop an efficient sampling method based on the sequential Monte Carlo principle. Our method considers excluded volume effect. It is general and can be applied to calculating entropy of loops with longer length and arbitrary complexity. For loops of short length, our results are in good agreement with a recent theoretical model and experimental measurement. For long loops, our estimated entropy of hairpin loops is in excellent agreement with the Jacobson-Stockmayer extrapolation model. However, for bulge loops and more complex secondary structures such as internal and multibranch loops, we find that the Jacobson-Stockmayer extrapolation model has large errors. Based on estimated entropy, we have developed empirical formulae for accurate calculation of entropy of long loops in different secondary structures. Our study on the effect of asymmetric size of loops suggest that loop entropy of internal loops is largely determined by the total loop length, and is only marginally affected by the asymmetric size of the two loops. Our finding suggests that the significant asymmetric effects of loop length in internal loops measured by experiments are likely to be partially enthalpic. Our method can be applied to develop improved energy parameters important for studying RNA stability and folding, and for predicting RNA secondary and tertiary structures. The discrete model and the program used to calculate loop entropy can be downloaded at http://gila.bioengr.uic.edu/resources/RNA.html. PMID:18376982
Accurate response surface approximations for weight equations based on structural optimization
NASA Astrophysics Data System (ADS)
Papila, Melih
Accurate weight prediction methods are vitally important for aircraft design optimization. Therefore, designers seek weight prediction techniques with low computational cost and high accuracy, and usually require a compromise between the two. The compromise can be achieved by combining stress analysis and response surface (RS) methodology. While stress analysis provides accurate weight information, RS techniques help to transmit effectively this information to the optimization procedure. The focus of this dissertation is structural weight equations in the form of RS approximations and their accuracy when fitted to results of structural optimizations that are based on finite element analyses. Use of RS methodology filters out the numerical noise in structural optimization results and provides a smooth weight function that can easily be used in gradient-based configuration optimization. In engineering applications RS approximations of low order polynomials are widely used, but the weight may not be modeled well by low-order polynomials, leading to bias errors. In addition, some structural optimization results may have high-amplitude errors (outliers) that may severely affect the accuracy of the weight equation. Statistical techniques associated with RS methodology are sought in order to deal with these two difficulties: (1) high-amplitude numerical noise (outliers) and (2) approximation model inadequacy. The investigation starts with reducing approximation error by identifying and repairing outliers. A potential reason for outliers in optimization results is premature convergence, and outliers of such nature may be corrected by employing different convergence settings. It is demonstrated that outlier repair can lead to accuracy improvements over the more standard approach of removing outliers. The adequacy of approximation is then studied by a modified lack-of-fit approach, and RS errors due to the approximation model are reduced by using higher order polynomials. In
Estimating Building Simulation Parameters via Bayesian Structure Learning
Edwards, Richard E; New, Joshua Ryan; Parker, Lynne Edwards
2013-01-01
Many key building design policies are made using sophisticated computer simulations such as EnergyPlus (E+), the DOE flagship whole-building energy simulation engine. E+ and other sophisticated computer simulations have several major problems. The two main issues are 1) gaps between the simulation model and the actual structure, and 2) limitations of the modeling engine's capabilities. Currently, these problems are addressed by having an engineer manually calibrate simulation parameters to real world data or using algorithmic optimization methods to adjust the building parameters. However, some simulations engines, like E+, are computationally expensive, which makes repeatedly evaluating the simulation engine costly. This work explores addressing this issue by automatically discovering the simulation's internal input and output dependencies from 20 Gigabytes of E+ simulation data, future extensions will use 200 Terabytes of E+ simulation data. The model is validated by inferring building parameters for E+ simulations with ground truth building parameters. Our results indicate that the model accurately represents parameter means with some deviation from the means, but does not support inferring parameter values that exist on the distribution's tail.
Identification of bolted lap joints parameters in assembled structures
NASA Astrophysics Data System (ADS)
Ahmadian, Hamid; Jalali, Hassan
2007-02-01
Bolted lap joints have significant influence on the dynamical behaviour of the assembled structures due to creation of strong local flexibility and damping. In modelling the dynamical behaviour of assembled structures the joint interface model must be represented accurately. A nonlinear model for bolted lap joints and interfaces is proposed capable of representing the dominant physics involved in the joint such as micro/macro-slip. The joint interface is modelled using a combination of linear and nonlinear springs and a damper to simulate the damping effects of the joint. An estimate of the response of the structure with a nonlinear model for the bolted joint under external excitations is obtained using the method of multiple scales. The parameters of the model, i.e. the spring constants and the damper coefficient, are functions of normal and tangential stresses at the joint interface and are identified by minimizing the difference between the model predictions and the experimentally measured data.
Subramanian, Swetha; Mast, T Douglas
2015-10-01
Computational finite element models are commonly used for the simulation of radiofrequency ablation (RFA) treatments. However, the accuracy of these simulations is limited by the lack of precise knowledge of tissue parameters. In this technical note, an inverse solver based on the unscented Kalman filter (UKF) is proposed to optimize values for specific heat, thermal conductivity, and electrical conductivity resulting in accurately simulated temperature elevations. A total of 15 RFA treatments were performed on ex vivo bovine liver tissue. For each RFA treatment, 15 finite-element simulations were performed using a set of deterministically chosen tissue parameters to estimate the mean and variance of the resulting tissue ablation. The UKF was implemented as an inverse solver to recover the specific heat, thermal conductivity, and electrical conductivity corresponding to the measured area of the ablated tissue region, as determined from gross tissue histology. These tissue parameters were then employed in the finite element model to simulate the position- and time-dependent tissue temperature. Results show good agreement between simulated and measured temperature. PMID:26352462
NASA Astrophysics Data System (ADS)
Subramanian, Swetha; Mast, T. Douglas
2015-09-01
Computational finite element models are commonly used for the simulation of radiofrequency ablation (RFA) treatments. However, the accuracy of these simulations is limited by the lack of precise knowledge of tissue parameters. In this technical note, an inverse solver based on the unscented Kalman filter (UKF) is proposed to optimize values for specific heat, thermal conductivity, and electrical conductivity resulting in accurately simulated temperature elevations. A total of 15 RFA treatments were performed on ex vivo bovine liver tissue. For each RFA treatment, 15 finite-element simulations were performed using a set of deterministically chosen tissue parameters to estimate the mean and variance of the resulting tissue ablation. The UKF was implemented as an inverse solver to recover the specific heat, thermal conductivity, and electrical conductivity corresponding to the measured area of the ablated tissue region, as determined from gross tissue histology. These tissue parameters were then employed in the finite element model to simulate the position- and time-dependent tissue temperature. Results show good agreement between simulated and measured temperature.
Bellili, A; Linguerri, R; Hochlaf, M; Puzzarini, C
2015-11-14
In an effort to provide an accurate structural and spectroscopic characterization of acetyl cyanide, its two enolic isomers and the corresponding cationic species, state-of-the-art computational methods, and approaches have been employed. The coupled-cluster theory including single and double excitations together with a perturbative treatment of triples has been used as starting point in composite schemes accounting for extrapolation to the complete basis-set limit as well as core-valence correlation effects to determine highly accurate molecular structures, fundamental vibrational frequencies, and rotational parameters. The available experimental data for acetyl cyanide allowed us to assess the reliability of our computations: structural, energetic, and spectroscopic properties have been obtained with an overall accuracy of about, or better than, 0.001 Å, 2 kcal/mol, 1-10 MHz, and 11 cm(-1) for bond distances, adiabatic ionization potentials, rotational constants, and fundamental vibrational frequencies, respectively. We are therefore confident that the highly accurate spectroscopic data provided herein can be useful for guiding future experimental investigations and/or astronomical observations. PMID:26567669
Gaggiotti, Oscar E
2010-11-01
Ever since the introduction of allozymes in the 1960s, evolutionary biologists and ecologists have continued to search for more powerful molecular markers to estimate important parameters such as effective population size and migration rates and to make inferences about the demographic history of populations, the relationships between individuals and the genetic architecture of phenotypic variation (Bensch & Akesson 2005; Bonin et al. 2007). Choosing a marker requires a thorough consideration of the trade-offs associated with the different techniques and the type of data obtained from them. Some markers can be very informative but require substantial amounts of start-up time (e.g. microsatellites), while others require very little time but are much less polymorphic. Amplified fragment length polymorphism (AFLP) is a firmly established molecular marker technique that falls in this latter category. AFLPs are widely distributed throughout the genome and can be used on organisms for which there is no a priori sequence information (Meudt & Clarke 2007). These properties together with their moderate cost and short start-up time have made them the method of choice for many molecular ecology studies of wild species (Bensch & Akesson 2005). However, they have a major disadvantage, they are dominant. This represents a very important limitation because many statistical genetics methods appropriate for molecular ecology studies require the use of codominant markers. In this issue, Foll et al. (2010) present an innovative hierarchical Bayesian method that overcomes this limitation. The proposed approach represents a comprehensive statistical treatment of the fluorescence of AFLP bands and leads to accurate inferences about the genetic structure of natural populations. Besides allowing a quasi-codominant treatment of AFLPs, this new method also solves the difficult problems posed by subjectivity in the scoring of AFLP bands. PMID:20958811
NASA Astrophysics Data System (ADS)
Iorio, L.
2016-01-01
By using the most recently published Doppler tomography measurements and accurate theoretical modelling of the oblateness-driven orbital precessions, we tightly constrain some of the physical and orbital parameters of the planetary system hosted by the fast rotating star WASP-33. In particular, the measurements of the orbital inclination ip to the plane of the sky and of the sky-projected spin-orbit misalignment λ at two epochs about six years apart allowed for the determination of the longitude of the ascending node Ω and of the orbital inclination I to the apparent equatorial plane at the same epochs. As a consequence, average rates of change dot{Ω }_exp, dot{I}_exp of this two orbital elements, accurate to a ≈10-2 deg yr-1 level, were calculated as well. By comparing them to general theoretical expressions dot{Ω }_{J_2}, dot{I}_{J_2} for their precessions induced by an oblate star whose symmetry axis is arbitrarily oriented, we were able to determine the angle i⋆ between the line of sight the star's spin {S}^{star } and its first even zonal harmonic J_2^{star } obtaining i^{star } = {142}^{+10}_{-11} deg, J_2^{star } = 2.1^{+0.8}_{-0.5}times; 10^{-4}. As a by-product, the angle between {S}^{star } and the orbital angular momentum L is as large as about ψ ≈ 100 ° psi; ^{2008} = 99^{+5}_{-4} deg, ψ ^{{2014}} = 103^{+5}_{-4} deg and changes at a rate dot{ψ }= 0.{7}^{+1.5}_{-1.6} deg {yr}^{-1}. The predicted general relativistic Lense-Thirring precessions, of the order of ≈10-3deg yr-1, are, at present, about one order of magnitude below the measurability threshold.
Accurate structural study of langasite-family Ca3TaGa3Si2O14 crystal
NASA Astrophysics Data System (ADS)
Dudka, A. P.
2016-03-01
An accurate X-ray diffraction study of Ca3TaGa3Si2O14 single crystal has been performed using two datasets obtained on a diffractometer equipped with a CCD area detector ( a = 8.1056(2) Å, c = 4.9800(1) Å, sp. gr. P321, Z = 1, R/ wR = 0.486/0.488%). A model structure is determined which is characterized by a high degree of reproducibility of structural parameters. Each site in Ca3TaGa3Si2O14 is occupied by atoms of only one type. Nevertheless, its structural feature is asymmetric disordering of sites of Ca, Ta, Ga, and two out of three oxygen atoms occupying special and general sites. A transition of some part of Ca atoms (~3%) from 3 e sites on the twofold symmetry axis to general 6 g sites is revealed.
Accurate Electron Affinity of Iron and Fine Structures of Negative Iron ions
Chen, Xiaolin; Luo, Zhihong; Li, Jiaming; Ning, Chuangang
2016-01-01
Ionization potential (IP) is defined as the amount of energy required to remove the most loosely bound electron of an atom, while electron affinity (EA) is defined as the amount of energy released when an electron is attached to a neutral atom. Both IP and EA are critical for understanding chemical properties of an element. In contrast to accurate IPs and structures of neutral atoms, EAs and structures of negative ions are relatively unexplored, especially for the transition metal anions. Here, we report the accurate EA value of Fe and fine structures of Fe− using the slow electron velocity imaging method. These measurements yield a very accurate EA value of Fe, 1235.93(28) cm−1 or 153.236(34) meV. The fine structures of Fe− were also successfully resolved. The present work provides a reliable benchmark for theoretical calculations, and also paves the way for improving the EA measurements of other transition metal atoms to the sub cm−1 accuracy. PMID:27138292
Accurate Electron Affinity of Iron and Fine Structures of Negative Iron ions.
Chen, Xiaolin; Luo, Zhihong; Li, Jiaming; Ning, Chuangang
2016-01-01
Ionization potential (IP) is defined as the amount of energy required to remove the most loosely bound electron of an atom, while electron affinity (EA) is defined as the amount of energy released when an electron is attached to a neutral atom. Both IP and EA are critical for understanding chemical properties of an element. In contrast to accurate IPs and structures of neutral atoms, EAs and structures of negative ions are relatively unexplored, especially for the transition metal anions. Here, we report the accurate EA value of Fe and fine structures of Fe(-) using the slow electron velocity imaging method. These measurements yield a very accurate EA value of Fe, 1235.93(28) cm(-1) or 153.236(34) meV. The fine structures of Fe(-) were also successfully resolved. The present work provides a reliable benchmark for theoretical calculations, and also paves the way for improving the EA measurements of other transition metal atoms to the sub cm(-1) accuracy. PMID:27138292
Accurate Electron Affinity of Iron and Fine Structures of Negative Iron ions
NASA Astrophysics Data System (ADS)
Chen, Xiaolin; Luo, Zhihong; Li, Jiaming; Ning, Chuangang
2016-05-01
Ionization potential (IP) is defined as the amount of energy required to remove the most loosely bound electron of an atom, while electron affinity (EA) is defined as the amount of energy released when an electron is attached to a neutral atom. Both IP and EA are critical for understanding chemical properties of an element. In contrast to accurate IPs and structures of neutral atoms, EAs and structures of negative ions are relatively unexplored, especially for the transition metal anions. Here, we report the accurate EA value of Fe and fine structures of Fe‑ using the slow electron velocity imaging method. These measurements yield a very accurate EA value of Fe, 1235.93(28) cm‑1 or 153.236(34) meV. The fine structures of Fe‑ were also successfully resolved. The present work provides a reliable benchmark for theoretical calculations, and also paves the way for improving the EA measurements of other transition metal atoms to the sub cm‑1 accuracy.
Determination of structure parameters in molecular tunnelling ionisation model
NASA Astrophysics Data System (ADS)
Wang, Jun-Ping; Zhao, Song-Feng; Zhang, Cai-Rong; Li, Wei; Zhou, Xiao-Xin
2014-04-01
We extracted the accurate structure parameters in a molecular tunnelling ionisation model (the so-called MO-ADK model) for 23 selected linear molecules including some inner orbitals. The molecular wave functions with the correct asymptotic behaviour are obtained by solving the time-independent Schrödinger equation with B-spline functions and molecular potentials numerically constructed using the modified Leeuwen-Baerends (LBα) model. We show that the orientation-dependent ionisation rate reflects the shape of the ionising orbitals in general. The influences of the Stark shifts of the energy levels on the orientation-dependent ionisation rates of the polar molecules are studied. We also examine the angle-dependent ionisation rates (or probabilities) based on the MO-ADK model by comparing with the molecular strong-field approximation calculations and with recent experimental measurements.
NASA Astrophysics Data System (ADS)
Skone, Jonathan; Govoni, Marco; Galli, Giulia
Dielectric-dependent hybrid [DDH] functionals have recently been shown to yield highly accurate energy gaps and dielectric constants for a wide variety of solids, at a computational cost considerably less than standard GW calculations. The fraction of exact exchange included in the definition of DDH functionals depends (self-consistently) on the dielectric constant of the material. In the present talk we introduce a range-separated (RS) version of DDH functionals where short and long-range components are matched using material dependent, non-empirical parameters. Comparing with state of the art GW calculations and experiment, we show that such RS hybrids yield accurate electronic properties of both molecules and solids, including energy gaps, photoelectron spectra and absolute ionization potentials. This work was supported by NSF-CCI Grant Number NSF-CHE-0802907 and DOE-BES.
Parameter identification of material constants in a composite shell structure
Martinez, D.R.; Carne, T.G.
1988-01-01
One of the basic requirements in engineering analysis is the development of a mathematical model describing the system. Frequently, comparisons with test data are used as a measurement of the adequacy of the model. An attempt is typically made to update or improve the model to provide a test-verified analysis tool. System identification provides a systematic procedure for accomplishing this task. The terms system identification, parameter estimation, and model correlation all refer to techniques that use test information to update or verify mathematical models. The goal of system identification is to improve the correlation of model predictions with measured test data, and produce accurate, predictive models. For nonmetallic structures the modeling task is often difficult due to uncertainties in the elastic constants. In this work a parameter identification procedure was used to determine the elastic constants of a cylindrical, graphite epoxy composite shell. A finite element model of the shell was created, which included uncertain orthotropic elastic constants. A modal survey test was then performed on the shell. The resulting modal data, along with the finite element model of the shell, were used in a Bayes estimation algorithm. This permitted the use of covariance matrices to weight the confidence in the initial parameter values as well as confidence in the measured test data. The estimation procedure also employed the concept of successive linearization to obtain an approximate solution to the original nonlinear estimation problem. 17 refs., 7 figs.
CoMOGrad and PHOG: From Computer Vision to Fast and Accurate Protein Tertiary Structure Retrieval
Karim, Rezaul; Aziz, Mohd. Momin Al; Shatabda, Swakkhar; Rahman, M. Sohel; Mia, Md. Abul Kashem; Zaman, Farhana; Rakin, Salman
2015-01-01
The number of entries in a structural database of proteins is increasing day by day. Methods for retrieving protein tertiary structures from such a large database have turn out to be the key to comparative analysis of structures that plays an important role to understand proteins and their functions. In this paper, we present fast and accurate methods for the retrieval of proteins having tertiary structures similar to a query protein from a large database. Our proposed methods borrow ideas from the field of computer vision. The speed and accuracy of our methods come from the two newly introduced features- the co-occurrence matrix of the oriented gradient and pyramid histogram of oriented gradient- and the use of Euclidean distance as the distance measure. Experimental results clearly indicate the superiority of our approach in both running time and accuracy. Our method is readily available for use from this website: http://research.buet.ac.bd:8080/Comograd/. PMID:26293226
CoMOGrad and PHOG: From Computer Vision to Fast and Accurate Protein Tertiary Structure Retrieval.
Karim, Rezaul; Aziz, Mohd Momin Al; Shatabda, Swakkhar; Rahman, M Sohel; Mia, Md Abul Kashem; Zaman, Farhana; Rakin, Salman
2015-01-01
The number of entries in a structural database of proteins is increasing day by day. Methods for retrieving protein tertiary structures from such a large database have turn out to be the key to comparative analysis of structures that plays an important role to understand proteins and their functions. In this paper, we present fast and accurate methods for the retrieval of proteins having tertiary structures similar to a query protein from a large database. Our proposed methods borrow ideas from the field of computer vision. The speed and accuracy of our methods come from the two newly introduced features- the co-occurrence matrix of the oriented gradient and pyramid histogram of oriented gradient- and the use of Euclidean distance as the distance measure. Experimental results clearly indicate the superiority of our approach in both running time and accuracy. Our method is readily available for use from this website: http://research.buet.ac.bd:8080/Comograd/. PMID:26293226
EXAFS Energy Shift and Structural Parameters
NASA Astrophysics Data System (ADS)
Kelly, Shelly D.; Ravel, Bruce
2007-02-01
In EXAFS analysis, the energy shift parameter is used to align the theoretical calculated spectrum to the energy grid of the measured spectrum. Unrealistically large energy shift values, sometimes in excess of 20 eV, are at times published in research articles. We therefore see the need for a concise discussion of the EXAFS energy shift parameter. This paper is intended as a learning tool for the proper alignment of theory to measured EXAFS spectra and proper interpretation of the energy shift parameter.
Reconstruction of biofilm images: combining local and global structural parameters.
Resat, Haluk; Renslow, Ryan S; Beyenal, Haluk
2014-10-01
Digitized images can be used for quantitative comparison of biofilms grown under different conditions. Using biofilm image reconstruction, it was previously found that biofilms with a completely different look can have nearly identical structural parameters and that the most commonly utilized global structural parameters were not sufficient to uniquely define these biofilms. Here, additional local and global parameters are introduced to show that these parameters considerably increase the reliability of the image reconstruction process. Assessment using human evaluators indicated that the correct identification rate of the reconstructed images increased from 50% to 72% with the introduction of the new parameters into the reconstruction procedure. An expanded set of parameters especially improved the identification of biofilm structures with internal orientational features and of structures in which colony sizes and spatial locations varied. Hence, the newly introduced structural parameter sets helped to better classify the biofilms by incorporating finer local structural details into the reconstruction process. PMID:25377487
Reconstruction of biofilm images: combining local and global structural parameters
Resat, Haluk; Renslow, Ryan S.; Beyenal, Haluk
2014-11-07
Digitized images can be used for quantitative comparison of biofilms grown under different conditions. Using biofilm image reconstruction, it was previously found that biofilms with a completely different look can have nearly identical structural parameters and that the most commonly utilized global structural parameters were not sufficient to uniquely define these biofilms. Here, additional local and global parameters are introduced to show that these parameters considerably increase the reliability of the image reconstruction process. Assessment using human evaluators indicated that the correct identification rate of the reconstructed images increased from 50% to 72% with the introduction of the new parameters into the reconstruction procedure. An expanded set of parameters especially improved the identification of biofilm structures with internal orientational features and of structures in which colony sizes and spatial locations varied. Hence, the newly introduced structural parameter sets helped to better classify the biofilms by incorporating finer local structural details into the reconstruction process.
Geng, Hao; Jiang, Fan; Wu, Yun-Dong
2016-05-19
Cyclic peptides (CPs) are promising candidates for drugs, chemical biology tools, and self-assembling nanomaterials. However, the development of reliable and accurate computational methods for their structure prediction has been challenging. Here, 20 all-trans CPs of 5-12 residues selected from Cambridge Structure Database have been simulated using replica-exchange molecular dynamics with four different force fields. Our recently developed residue-specific force fields RSFF1 and RSFF2 can correctly identify the crystal-like conformations of more than half CPs as the most populated conformation. The RSFF2 performs the best, which consistently predicts the crystal structures of 17 out of 20 CPs with rmsd < 1.1 Å. We also compared the backbone (ϕ, ψ) sampling of residues in CPs with those in short linear peptides and in globular proteins. In general, unlike linear peptides, CPs have local conformational free energies and entropies quite similar to globular proteins. PMID:27128113
Brandenburg, Jan Gerit; Caldeweyher, Eike; Grimme, Stefan
2016-06-21
We extend the recently introduced PBEh-3c global hybrid density functional [S. Grimme et al., J. Chem. Phys., 2015, 143, 054107] by a screened Fock exchange variant based on the Henderson-Janesko-Scuseria exchange hole model. While the excellent performance of the global hybrid is maintained for small covalently bound molecules, its performance for computed condensed phase mass densities is further improved. Most importantly, a speed up of 30 to 50% can be achieved and especially for small orbital energy gap cases, the method is numerically much more robust. The latter point is important for many applications, e.g., for metal-organic frameworks, organic semiconductors, or protein structures. This enables an accurate density functional based electronic structure calculation of a full DNA helix structure on a single core desktop computer which is presented as an example in addition to comprehensive benchmark results. PMID:27240749
Parameter identification of material constants in a composite shell structure
NASA Technical Reports Server (NTRS)
Martinez, David R.; Carne, Thomas G.
1988-01-01
One of the basic requirements in engineering analysis is the development of a mathematical model describing the system. Frequently comparisons with test data are used as a measurement of the adequacy of the model. An attempt is typically made to update or improve the model to provide a test verified analysis tool. System identification provides a systematic procedure for accomplishing this task. The terms system identification, parameter estimation, and model correlation all refer to techniques that use test information to update or verify mathematical models. The goal of system identification is to improve the correlation of model predictions with measured test data, and produce accurate, predictive models. For nonmetallic structures the modeling task is often difficult due to uncertainties in the elastic constants. A finite element model of the shell was created, which included uncertain orthotropic elastic constants. A modal survey test was then performed on the shell. The resulting modal data, along with the finite element model of the shell, were used in a Bayes estimation algorithm. This permitted the use of covariance matrices to weight the confidence in the initial parameter values as well as confidence in the measured test data. The estimation procedure also employed the concept of successive linearization to obtain an approximate solution to the original nonlinear estimation problem.
NASA Astrophysics Data System (ADS)
Eaton, M.; Pearson, M.; Lee, W.; Pullin, R.
2015-07-01
The ability to accurately locate damage in any given structure is a highly desirable attribute for an effective structural health monitoring system and could help to reduce operating costs and improve safety. This becomes a far greater challenge in complex geometries and materials, such as modern composite airframes. The poor translation of promising laboratory based SHM demonstrators to industrial environments forms a barrier to commercial up take of technology. The acoustic emission (AE) technique is a passive NDT method that detects elastic stress waves released by the growth of damage. It offers very sensitive damage detection, using a sparse array of sensors to detect and globally locate damage within a structure. However its application to complex structures commonly yields poor accuracy due to anisotropic wave propagation and the interruption of wave propagation by structural features such as holes and thickness changes. This work adopts an empirical mapping technique for AE location, known as Delta T Mapping, which uses experimental training data to account for such structural complexities. The technique is applied to a complex geometry composite aerospace structure undergoing certification testing. The component consists of a carbon fibre composite tube with varying wall thickness and multiple holes, that was loaded under bending. The damage location was validated using X-ray CT scanning and the Delta T Mapping technique was shown to improve location accuracy when compared with commercial algorithms. The onset and progression of damage were monitored throughout the test and used to inform future design iterations.
Accurate macromolecular structures using minimal measurements from X-ray free-electron lasers
Hattne, Johan; Echols, Nathaniel; Tran, Rosalie; Kern, Jan; Gildea, Richard J.; Brewster, Aaron S.; Alonso-Mori, Roberto; Glöckner, Carina; Hellmich, Julia; Laksmono, Hartawan; Sierra, Raymond G.; Lassalle-Kaiser, Benedikt; Lampe, Alyssa; Han, Guangye; Gul, Sheraz; DiFiore, Dörte; Milathianaki, Despina; Fry, Alan R.; Miahnahri, Alan; White, William E.; Schafer, Donald W.; Seibert, M. Marvin; Koglin, Jason E.; Sokaras, Dimosthenis; Weng, Tsu-Chien; Sellberg, Jonas; Latimer, Matthew J.; Glatzel, Pieter; Zwart, Petrus H.; Grosse-Kunstleve, Ralf W.; Bogan, Michael J.; Messerschmidt, Marc; Williams, Garth J.; Boutet, Sébastien; Messinger, Johannes; Zouni, Athina; Yano, Junko; Bergmann, Uwe; Yachandra, Vittal K.; Adams, Paul D.; Sauter, Nicholas K.
2014-01-01
X-ray free-electron laser (XFEL) sources enable the use of crystallography to solve three-dimensional macromolecular structures under native conditions and free from radiation damage. Results to date, however, have been limited by the challenge of deriving accurate Bragg intensities from a heterogeneous population of microcrystals, while at the same time modeling the X-ray spectrum and detector geometry. Here we present a computational approach designed to extract statistically significant high-resolution signals from fewer diffraction measurements. PMID:24633409
Reeves, Patrick A.; Richards, Christopher M.
2009-01-01
Background Accurate inference of genetic discontinuities between populations is an essential component of intraspecific biodiversity and evolution studies, as well as associative genetics. The most widely-used methods to infer population structure are model-based, Bayesian MCMC procedures that minimize Hardy-Weinberg and linkage disequilibrium within subpopulations. These methods are useful, but suffer from large computational requirements and a dependence on modeling assumptions that may not be met in real data sets. Here we describe the development of a new approach, PCO-MC, which couples principal coordinate analysis to a clustering procedure for the inference of population structure from multilocus genotype data. Methodology/Principal Findings PCO-MC uses data from all principal coordinate axes simultaneously to calculate a multidimensional “density landscape”, from which the number of subpopulations, and the membership within subpopulations, is determined using a valley-seeking algorithm. Using extensive simulations, we show that this approach outperforms a Bayesian MCMC procedure when many loci (e.g. 100) are sampled, but that the Bayesian procedure is marginally superior with few loci (e.g. 10). When presented with sufficient data, PCO-MC accurately delineated subpopulations with population Fst values as low as 0.03 (G'st>0.2), whereas the limit of resolution of the Bayesian approach was Fst = 0.05 (G'st>0.35). Conclusions/Significance We draw a distinction between population structure inference for describing biodiversity as opposed to Type I error control in associative genetics. We suggest that discrete assignments, like those produced by PCO-MC, are appropriate for circumscribing units of biodiversity whereas expression of population structure as a continuous variable is more useful for case-control correction in structured association studies. PMID:19172174
NASA Astrophysics Data System (ADS)
Kimura, H.; Asano, Y.; Matsumoto, T.
2012-12-01
The rapid determination of hypocentral parameters and their transmission to the public are valuable components of disaster mitigation. We have operated an automatic system for this purpose—termed the Accurate and QUick Analysis system for source parameters (AQUA)—since 2005 (Matsumura et al., 2006). In this system, the initial hypocenter, the moment tensor (MT), and the centroid moment tensor (CMT) solutions are automatically determined and posted on the NIED Hi-net Web site (www.hinet.bosai.go.jp). This paper describes improvements made to the AQUA to overcome limitations that became apparent after the 2011 Tohoku Earthquake (05:46:17, March 11, 2011 in UTC). The improvements included the processing of NIED F-net velocity-type strong motion records, because NIED F-net broadband seismographs are saturated for great earthquakes such as the 2011 Tohoku Earthquake. These velocity-type strong motion seismographs provide unsaturated records not only for the 2011 Tohoku Earthquake, but also for recording stations located close to the epicenters of M>7 earthquakes. We used 0.005-0.020 Hz records for M>7.5 earthquakes, in contrast to the 0.01-0.05 Hz records employed in the original system. The initial hypocenters determined based on arrival times picked by using seismograms recorded by NIED Hi-net stations can have large errors in terms of magnitude and hypocenter location, especially for great earthquakes or earthquakes located far from the onland Hi-net network. The size of the 2011 Tohoku Earthquake was initially underestimated in the AQUA to be around M5 at the initial stage of rupture. Numerous aftershocks occurred at the outer rise east of the Japan trench, where a great earthquake is anticipated to occur. Hence, we modified the system to repeat the MT analyses assuming a larger size, for all earthquakes for which the magnitude was initially underestimated. We also broadened the search range of centroid depth for earthquakes located far from the onland Hi
"Cosmological Parameters from Large Scale Structure"
NASA Technical Reports Server (NTRS)
Hamilton, A. J. S.
2005-01-01
This grant has provided primary support for graduate student Mark Neyrinck, and some support for the PI and for colleague Nick Gnedin, who helped co-supervise Neyrinck. This award had two major goals. First, to continue to develop and apply methods for measuring galaxy power spectra on large, linear scales, with a view to constraining cosmological parameters. And second, to begin try to understand galaxy clustering at smaller. nonlinear scales well enough to constrain cosmology from those scales also. Under this grant, the PI and collaborators, notably Max Tegmark. continued to improve their technology for measuring power spectra from galaxy surveys at large, linear scales. and to apply the technology to surveys as the data become available. We believe that our methods are best in the world. These measurements become the foundation from which we and other groups measure cosmological parameters.
NASA Astrophysics Data System (ADS)
Skone, Jonathan; Govoni, Marco; Galli, Giulia
2015-03-01
Building upon a recently proposed self-consistent hybrid (sc-hybrid) functional, where the optimal dielectric screening is included self-consistently, we propose an improved form by incorporating range-separation of the exchange part. We discuss the choice of the non-empirical parameters defining range separation, and we present results for condensed media including semiconductors, amorphous insulators, and molecular crystals. We find that the range-separated sc-hybrid functional further improves upon the electronic gaps obtained with full-range sc-hybrids, thus providing an accurate functional for high throughput band gap engineering. This work was supported by NSF-CCI Grant Number NSF-CHE-0802907 and ARL Grant Number W911NF-12-2-0023.
DR-TAMAS: Diffeomorphic Registration for Tensor Accurate Alignment of Anatomical Structures.
Irfanoglu, M Okan; Nayak, Amritha; Jenkins, Jeffrey; Hutchinson, Elizabeth B; Sadeghi, Neda; Thomas, Cibu P; Pierpaoli, Carlo
2016-05-15
In this work, we propose DR-TAMAS (Diffeomorphic Registration for Tensor Accurate alignMent of Anatomical Structures), a novel framework for intersubject registration of Diffusion Tensor Imaging (DTI) data sets. This framework is optimized for brain data and its main goal is to achieve an accurate alignment of all brain structures, including white matter (WM), gray matter (GM), and spaces containing cerebrospinal fluid (CSF). Currently most DTI-based spatial normalization algorithms emphasize alignment of anisotropic structures. While some diffusion-derived metrics, such as diffusion anisotropy and tensor eigenvector orientation, are highly informative for proper alignment of WM, other tensor metrics such as the trace or mean diffusivity (MD) are fundamental for a proper alignment of GM and CSF boundaries. Moreover, it is desirable to include information from structural MRI data, e.g., T1-weighted or T2-weighted images, which are usually available together with the diffusion data. The fundamental property of DR-TAMAS is to achieve global anatomical accuracy by incorporating in its cost function the most informative metrics locally. Another important feature of DR-TAMAS is a symmetric time-varying velocity-based transformation model, which enables it to account for potentially large anatomical variability in healthy subjects and patients. The performance of DR-TAMAS is evaluated with several data sets and compared with other widely-used diffeomorphic image registration techniques employing both full tensor information and/or DTI-derived scalar maps. Our results show that the proposed method has excellent overall performance in the entire brain, while being equivalent to the best existing methods in WM. PMID:26931817
Cross-sectional structural parameters from densitometry
NASA Technical Reports Server (NTRS)
Cleek, Tammy M.; Whalen, Robert T.
2002-01-01
Bone densitometry has previously been used to obtain cross-sectional properties of bone from a single X-ray projection across the bone width. Using three unique projections, we have extended the method to obtain the principal area moments of inertia and orientations of the principal axes at each scan cross-section along the length of the scan. Various aluminum phantoms were used to examine scanner characteristics to develop the highest accuracy possible for in vitro non-invasive analysis of cross-sectional properties. Factors considered included X-ray photon energy, initial scan orientation, the angle spanned by the three scans (included angle), and I(min)/I(max) ratios. Principal moments of inertia were accurate to within +/-3.1% and principal angles were within +/-1 degrees of the expected value for phantoms scanned with included angles of 60 degrees and 90 degrees at the higher X-ray photon energy (140 kVp). Low standard deviations in the error (0.68-1.84%) also indicate high precision of calculated measurements with these included angles. Accuracy and precision decreased slightly when the included angle was reduced to 30 degrees. The method was then successfully applied to a pair of excised cadaveric tibiae. The accuracy and insensitivity of the algorithms to cross-sectional shape and changing isotropy (I(min)/I(max)) values when various included angles are used make this technique viable for future in vivo studies.
Holton, James M; Classen, Scott; Frankel, Kenneth A; Tainer, John A
2014-01-01
In macromolecular crystallography, the agreement between observed and predicted structure factors (Rcryst and Rfree) is seldom better than 20%. This is much larger than the estimate of experimental error (Rmerge). The difference between Rcryst and Rmerge is the R-factor gap. There is no such gap in small-molecule crystallography, for which calculated structure factors are generally considered more accurate than the experimental measurements. Perhaps the true noise level of macromolecular data is higher than expected? Or is the gap caused by inaccurate phases that trap refined models in local minima? By generating simulated diffraction patterns using the program MLFSOM, and including every conceivable source of experimental error, we show that neither is the case. Processing our simulated data yielded values that were indistinguishable from those of real data for all crystallographic statistics except the final Rcryst and Rfree. These values decreased to 3.8% and 5.5% for simulated data, suggesting that the reason for high R-factors in macromolecular crystallography is neither experimental error nor phase bias, but rather an underlying inadequacy in the models used to explain our observations. The present inability to accurately represent the entire macromolecule with both its flexibility and its protein-solvent interface may be improved by synergies between small-angle X-ray scattering, computational chemistry and crystallography. The exciting implication of our finding is that macromolecular data contain substantial hidden and untapped potential to resolve ambiguities in the true nature of the nanoscale, a task that the second century of crystallography promises to fulfill. Database Coordinates and structure factors for the real data have been submitted to the Protein Data Bank under accession 4tws. PMID:25040949
NASA Technical Reports Server (NTRS)
Dominquez, Jesus A.; Tate, Lanetra C.; Wright, M. Clara; Caraccio, Anne
2013-01-01
Accomplishing the best-performing composite matrix (resin) requires that not only the processing method but also the cure cycle generate low-void-content structures. If voids are present, the performance of the composite matrix will be significantly reduced. This is usually noticed by significant reductions in matrix-dominated properties, such as compression and shear strength. Voids in composite materials are areas that are absent of the composite components: matrix and fibers. The characteristics of the voids and their accurate estimation are critical to determine for high performance composite structures. One widely used method of performing void analysis on a composite structure sample is acquiring optical micrographs or Scanning Electron Microscope (SEM) images of lateral sides of the sample and retrieving the void areas within the micrographs/images using an image analysis technique. Segmentation for the retrieval and subsequent computation of void areas within the micrographs/images is challenging as the gray-scaled values of the void areas are close to the gray-scaled values of the matrix leading to the need of manually performing the segmentation based on the histogram of the micrographs/images to retrieve the void areas. The use of an algorithm developed by NASA and based on Fuzzy Reasoning (FR) proved to overcome the difficulty of suitably differentiate void and matrix image areas with similar gray-scaled values leading not only to a more accurate estimation of void areas on composite matrix micrographs but also to a faster void analysis process as the algorithm is fully autonomous.
Accurate airway segmentation based on intensity structure analysis and graph-cut
NASA Astrophysics Data System (ADS)
Meng, Qier; Kitsaka, Takayuki; Nimura, Yukitaka; Oda, Masahiro; Mori, Kensaku
2016-03-01
This paper presents a novel airway segmentation method based on intensity structure analysis and graph-cut. Airway segmentation is an important step in analyzing chest CT volumes for computerized lung cancer detection, emphysema diagnosis, asthma diagnosis, and pre- and intra-operative bronchoscope navigation. However, obtaining a complete 3-D airway tree structure from a CT volume is quite challenging. Several researchers have proposed automated algorithms basically based on region growing and machine learning techniques. However these methods failed to detect the peripheral bronchi branches. They caused a large amount of leakage. This paper presents a novel approach that permits more accurate extraction of complex bronchial airway region. Our method are composed of three steps. First, the Hessian analysis is utilized for enhancing the line-like structure in CT volumes, then a multiscale cavity-enhancement filter is employed to detect the cavity-like structure from the previous enhanced result. In the second step, we utilize the support vector machine (SVM) to construct a classifier for removing the FP regions generated. Finally, the graph-cut algorithm is utilized to connect all of the candidate voxels to form an integrated airway tree. We applied this method to sixteen cases of 3D chest CT volumes. The results showed that the branch detection rate of this method can reach about 77.7% without leaking into the lung parenchyma areas.
conSSert: Consensus SVM Model for Accurate Prediction of Ordered Secondary Structure.
Kieslich, Chris A; Smadbeck, James; Khoury, George A; Floudas, Christodoulos A
2016-03-28
Accurate prediction of protein secondary structure remains a crucial step in most approaches to the protein-folding problem, yet the prediction of ordered secondary structure, specifically beta-strands, remains a challenge. We developed a consensus secondary structure prediction method, conSSert, which is based on support vector machines (SVM) and provides exceptional accuracy for the prediction of beta-strands with QE accuracy of over 0.82 and a Q2-EH of 0.86. conSSert uses as input probabilities for the three types of secondary structure (helix, strand, and coil) that are predicted by four top performing methods: PSSpred, PSIPRED, SPINE-X, and RAPTOR. conSSert was trained/tested using 4261 protein chains from PDBSelect25, and 8632 chains from PISCES. Further validation was performed using targets from CASP9, CASP10, and CASP11. Our data suggest that poor performance in strand prediction is likely a result of training bias and not solely due to the nonlocal nature of beta-sheet contacts. conSSert is freely available for noncommercial use as a webservice: http://ares.tamu.edu/conSSert/ . PMID:26928531
Critical parameters of superconducting materials and structures
Fluss, M.J.; Howell, R.H.; Sterne, P.A.; Dykes, J.W.; Mosley, W.D.; Chaiken, A.; Ralls, K.; Radousky, H.
1995-02-01
We report here the completion of a one year project to investigate the synthesis, electronic structure, defect structure, and physical transport properties of high temperature superconducting oxide materials. During the course of this project we produced some of the finest samples of single crystal detwinned YBa{sub 2}Cu{sub 3}O{sub 7}, and stoichiometrically perfect (Ba,K)BiO{sub 3}. We deduced the Fermi surface of YBa{sub 2}Cu{sub 3}O{sub 7}, (La,Sr){sub 2}CuO{sub 4}, and (Ba,K)BiO{sub 3} through the recording of the electron momentum density in these materials as measured by positron annihilation spectroscopy and angle resolved photoemission. We also performed extensive studies on Pr substituted (Y,Pr)Ba{sub 2}Cu{sub 3}O{sub 7} so as to further understand the origin of the electron pairing leading to superconductivity.
Kieslich, Chris A.; Tamamis, Phanourios; Guzman, Yannis A.; Onel, Melis; Floudas, Christodoulos A.
2016-01-01
HIV-1 entry into host cells is mediated by interactions between the V3-loop of viral glycoprotein gp120 and chemokine receptor CCR5 or CXCR4, collectively known as HIV-1 coreceptors. Accurate genotypic prediction of coreceptor usage is of significant clinical interest and determination of the factors driving tropism has been the focus of extensive study. We have developed a method based on nonlinear support vector machines to elucidate the interacting residue pairs driving coreceptor usage and provide highly accurate coreceptor usage predictions. Our models utilize centroid-centroid interaction energies from computationally derived structures of the V3-loop:coreceptor complexes as primary features, while additional features based on established rules regarding V3-loop sequences are also investigated. We tested our method on 2455 V3-loop sequences of various lengths and subtypes, and produce a median area under the receiver operator curve of 0.977 based on 500 runs of 10-fold cross validation. Our study is the first to elucidate a small set of specific interacting residue pairs between the V3-loop and coreceptors capable of predicting coreceptor usage with high accuracy across major HIV-1 subtypes. The developed method has been implemented as a web tool named CRUSH, CoReceptor USage prediction for HIV-1, which is available at http://ares.tamu.edu/CRUSH/. PMID:26859389
Advances in parameter estimation techniques applied to flexible structures
NASA Technical Reports Server (NTRS)
Maben, Egbert; Zimmerman, David C.
1994-01-01
In this work, various parameter estimation techniques are investigated in the context of structural system identification utilizing distributed parameter models and 'measured' time-domain data. Distributed parameter models are formulated using the PDEMOD software developed by Taylor. Enhancements made to PDEMOD for this work include the following: (1) a Wittrick-Williams based root solving algorithm; (2) a time simulation capability; and (3) various parameter estimation algorithms. The parameter estimations schemes will be contrasted using the NASA Mini-Mast as the focus structure.
Distributed parameter modeling of repeated truss structures
NASA Technical Reports Server (NTRS)
Wang, Han-Ching
1994-01-01
A new approach to find homogeneous models for beam-like repeated flexible structures is proposed which conceptually involves two steps. The first step involves the approximation of 3-D non-homogeneous model by a 1-D periodic beam model. The structure is modeled as a 3-D non-homogeneous continuum. The displacement field is approximated by Taylor series expansion. Then, the cross sectional mass and stiffness matrices are obtained by energy equivalence using their additive properties. Due to the repeated nature of the flexible bodies, the mass, and stiffness matrices are also periodic. This procedure is systematic and requires less dynamics detail. The first step involves the homogenization from a 1-D periodic beam model to a 1-D homogeneous beam model. The periodic beam model is homogenized into an equivalent homogeneous beam model using the additive property of compliance along the generic axis. The major departure from previous approaches in literature is using compliance instead of stiffness in homogenization. An obvious justification is that the stiffness is additive at each cross section but not along the generic axis. The homogenized model preserves many properties of the original periodic model.
NASA Astrophysics Data System (ADS)
Bloßfeld, Mathis; Panzetta, Francesca; Müller, Horst; Gerstl, Michael
2016-04-01
The GGOS vision is to integrate geometric and gravimetric observation techniques to estimate consistent geodetic-geophysical parameters. In order to reach this goal, the common estimation of station coordinates, Stokes coefficients and Earth Orientation Parameters (EOP) is necessary. Satellite Laser Ranging (SLR) provides the ability to study correlations between the different parameter groups since the observed satellite orbit dynamics are sensitive to the above mentioned geodetic parameters. To decrease the correlations, SLR observations to multiple satellites have to be combined. In this paper, we compare the estimated EOP of (i) single satellite SLR solutions and (ii) multi-satellite SLR solutions. Therefore, we jointly estimate station coordinates, EOP, Stokes coefficients and orbit parameters using different satellite constellations. A special focus in this investigation is put on the de-correlation of different geodetic parameter groups due to the combination of SLR observations. Besides SLR observations to spherical satellites (commonly used), we discuss the impact of SLR observations to non-spherical satellites such as, e.g., the JASON-2 satellite. The goal of this study is to discuss the existing parameter interactions and to present a strategy how to obtain reliable estimates of station coordinates, EOP, orbit parameter and Stokes coefficients in one common adjustment. Thereby, the benefits of a multi-satellite SLR solution are evaluated.
Fitmunk: improving protein structures by accurate, automatic modeling of side-chain conformations.
Porebski, Przemyslaw Jerzy; Cymborowski, Marcin; Pasenkiewicz-Gierula, Marta; Minor, Wladek
2016-02-01
Improvements in crystallographic hardware and software have allowed automated structure-solution pipelines to approach a near-`one-click' experience for the initial determination of macromolecular structures. However, in many cases the resulting initial model requires a laborious, iterative process of refinement and validation. A new method has been developed for the automatic modeling of side-chain conformations that takes advantage of rotamer-prediction methods in a crystallographic context. The algorithm, which is based on deterministic dead-end elimination (DEE) theory, uses new dense conformer libraries and a hybrid energy function derived from experimental data and prior information about rotamer frequencies to find the optimal conformation of each side chain. In contrast to existing methods, which incorporate the electron-density term into protein-modeling frameworks, the proposed algorithm is designed to take advantage of the highly discriminatory nature of electron-density maps. This method has been implemented in the program Fitmunk, which uses extensive conformational sampling. This improves the accuracy of the modeling and makes it a versatile tool for crystallographic model building, refinement and validation. Fitmunk was extensively tested on over 115 new structures, as well as a subset of 1100 structures from the PDB. It is demonstrated that the ability of Fitmunk to model more than 95% of side chains accurately is beneficial for improving the quality of crystallographic protein models, especially at medium and low resolutions. Fitmunk can be used for model validation of existing structures and as a tool to assess whether side chains are modeled optimally or could be better fitted into electron density. Fitmunk is available as a web service at http://kniahini.med.virginia.edu/fitmunk/server/ or at http://fitmunk.bitbucket.org/. PMID:26894674
Fitmunk: improving protein structures by accurate, automatic modeling of side-chain conformations
Porebski, Przemyslaw Jerzy; Cymborowski, Marcin; Pasenkiewicz-Gierula, Marta; Minor, Wladek
2016-01-01
Improvements in crystallographic hardware and software have allowed automated structure-solution pipelines to approach a near-‘one-click’ experience for the initial determination of macromolecular structures. However, in many cases the resulting initial model requires a laborious, iterative process of refinement and validation. A new method has been developed for the automatic modeling of side-chain conformations that takes advantage of rotamer-prediction methods in a crystallographic context. The algorithm, which is based on deterministic dead-end elimination (DEE) theory, uses new dense conformer libraries and a hybrid energy function derived from experimental data and prior information about rotamer frequencies to find the optimal conformation of each side chain. In contrast to existing methods, which incorporate the electron-density term into protein-modeling frameworks, the proposed algorithm is designed to take advantage of the highly discriminatory nature of electron-density maps. This method has been implemented in the program Fitmunk, which uses extensive conformational sampling. This improves the accuracy of the modeling and makes it a versatile tool for crystallographic model building, refinement and validation. Fitmunk was extensively tested on over 115 new structures, as well as a subset of 1100 structures from the PDB. It is demonstrated that the ability of Fitmunk to model more than 95% of side chains accurately is beneficial for improving the quality of crystallographic protein models, especially at medium and low resolutions. Fitmunk can be used for model validation of existing structures and as a tool to assess whether side chains are modeled optimally or could be better fitted into electron density. Fitmunk is available as a web service at http://kniahini.med.virginia.edu/fitmunk/server/ or at http://fitmunk.bitbucket.org/. PMID:26894674
NASA Astrophysics Data System (ADS)
Moore, Christopher; Hopkins, Matthew; Moore, Stan; Boerner, Jeremiah; Cartwright, Keith
2015-09-01
Simulation of breakdown is important for understanding and designing a variety of applications such as mitigating undesirable discharge events. Such simulations need to be accurate through early time arc initiation to late time stable arc behavior. Here we examine constraints on the timestep and mesh size required for arc simulations using the particle-in-cell (PIC) method with direct simulation Monte Carlo (DMSC) collisions. Accurate simulation of electron avalanche across a fixed voltage drop and constant neutral density (reduced field of 1000 Td) was found to require a timestep ~ 1/100 of the mean time between collisions and a mesh size ~ 1/25 the mean free path. These constraints are much smaller than the typical PIC-DSMC requirements for timestep and mesh size. Both constraints are related to the fact that charged particles are accelerated by the external field. Thus gradients in the electron energy distribution function can exist at scales smaller than the mean free path and these must be resolved by the mesh size for accurate collision rates. Additionally, the timestep must be small enough that the particle energy change due to the fields be small in order to capture gradients in the cross sections versus energy. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. DOE's National Nuclear Security Administration under Contract DE-AC04-94AL85000.
Bakker, Chris J G; de Leeuw, Hendrik; van de Maat, Gerrit H; van Gorp, Jetse S; Bouwman, Job G; Seevinck, Peter R
2013-01-01
Lack of spatial accuracy is a recognized problem in magnetic resonance imaging (MRI) which severely detracts from its value as a stand-alone modality for applications that put high demands on geometric fidelity, such as radiotherapy treatment planning and stereotactic neurosurgery. In this paper, we illustrate the potential and discuss the limitations of spectroscopic imaging as a tool for generating purely phase-encoded MR images and parameter maps that preserve the geometry of an object and allow localization of object features in world coordinates. Experiments were done on a clinical system with standard facilities for imaging and spectroscopy. Images were acquired with a regular spin echo sequence and a corresponding spectroscopic imaging sequence. In the latter, successive samples of the acquired echo were used for the reconstruction of a series of evenly spaced images in the time and frequency domain. Experiments were done with a spatial linearity phantom and a series of test objects representing a wide range of susceptibility- and chemical-shift-induced off-resonance conditions. In contrast to regular spin echo imaging, spectroscopic imaging was shown to be immune to off-resonance effects, such as those caused by field inhomogeneity, susceptibility, chemical shift, f(0) offset and field drift, and to yield geometrically accurate images and parameter maps that allowed object structures to be localized in world coordinates. From these illustrative examples and a discussion of the limitations of purely phase-encoded imaging techniques, it is concluded that spectroscopic imaging offers a fundamental solution to the geometric deficiencies of MRI which may evolve toward a practical solution when full advantage will be taken of current developments with regard to scan time reduction. This perspective is backed up by a demonstration of the significant scan time reduction that may be achieved by the use of compressed sensing for a simple phantom. PMID:22898694
Discretely disordered photonic bandgap structures: a more accurate invariant measure calculation
NASA Astrophysics Data System (ADS)
Kissel, Glen J.
2009-02-01
In the one-dimensional optical analog to Anderson localization, a periodically layered medium has one or more parameters randomly disordered. Such a randomized system can be modeled by an infinite product of 2x2 random transfer matrices with the upper Lyapunov exponent of the matrix product identified as the localization factor (inverse localization length) for the model. The theorem of Furstenberg allows us, at least theoretically, to calculate this upper Lyapunov exponent. In Furstenberg's formula we not only integrate with respect to the probability measure of the random matrices, but also with respect to the invariant probability measure of the direction of the vector propagated by the random matrices. This invariant measure is difficult to find analytically, and, as a result, the most successful approach is to determine the invariant measure numerically. A Monte Carlo simulation which uses accumulated bin counts to track the direction of the propagated vector through a long chain of random matrices does a good job of estimating the invariant probability measure, but with a level of uncertainty. A potentially more accurate numerical technique by Froyland and Aihara obtains the invariant measure as a left eigenvector of a large sparse matrix containing probability values determined by the action of the random matrices on input vectors. We first apply these two techniques to a random Fibonacci sequence whose Lyapunov exponent was determined by Viswanath. We then demonstrate these techniques on a quarter-wave stack model with binary discrete disorder in layer thickness, and compare results to the continuously disordered counterpart.
Can structured data fields accurately measure quality of care? The example of falls.
Ganz, David A; Almeida, Shone; Roth, Carol P; Reuben, David B; Wenger, Neil S
2012-01-01
By automating collection of data elements, electronic health records may simplify the process of measuring the quality of medical care. Using data from a quality improvement initiative in primary care medical groups, we sought to determine whether the quality of care for falls and fear of falling in outpatients aged 75 and older could be accurately measured solely from codable (non-free-text) data in a structured visit note. A traditional medical record review by trained abstractors served as the criterion standard. Among 215 patient records reviewed, we found a structured visit note in 54% of charts within 3 mo of the date patients had been identified as having falls or fear of falling. The reliability of an algorithm based on codable data was at least good (kappa of at least 0.61) compared with full medical record review for three care processes recommended for patients with two falls or one fall with injury in the past year: orthostatic vital signs, vision test/eye examination, and home safety evaluation. However, the automated algorithm routinely underestimated quality of care. Performance standards based on automated measurement of quality of care from electronic health records need to account for documentation occurring in nonstructured form. PMID:23408222
Toward Hamiltonian Adaptive QM/MM: Accurate Solvent Structures Using Many-Body Potentials.
Boereboom, Jelle M; Potestio, Raffaello; Donadio, Davide; Bulo, Rosa E
2016-08-01
Adaptive quantum mechanical (QM)/molecular mechanical (MM) methods enable efficient molecular simulations of chemistry in solution. Reactive subregions are modeled with an accurate QM potential energy expression while the rest of the system is described in a more approximate manner (MM). As solvent molecules diffuse in and out of the reactive region, they are gradually included into (and excluded from) the QM expression. It would be desirable to model such a system with a single adaptive Hamiltonian, but thus far this has resulted in distorted structures at the boundary between the two regions. Solving this long outstanding problem will allow microcanonical adaptive QM/MM simulations that can be used to obtain vibrational spectra and dynamical properties. The difficulty lies in the complex QM potential energy expression, with a many-body expansion that contains higher order terms. Here, we outline a Hamiltonian adaptive multiscale scheme within the framework of many-body potentials. The adaptive expressions are entirely general, and complementary to all standard (nonadaptive) QM/MM embedding schemes available. We demonstrate the merit of our approach on a molecular system defined by two different MM potentials (MM/MM'). For the long-range interactions a numerical scheme is used (particle mesh Ewald), which yields energy expressions that are many-body in nature. Our Hamiltonian approach is the first to provide both energy conservation and the correct solvent structure everywhere in this system. PMID:27332140
NASA Astrophysics Data System (ADS)
Brawand, Nicholas; Vörös, Márton; Govoni, Marco; Galli, Giulia
The accurate prediction of optoelectronic properties of molecules and solids is a persisting challenge for current density functional theory (DFT) based methods. We propose a hybrid functional where the mixing fraction of exact and local exchange is determined by a non-empirical, system dependent function. This functional yields ionization potentials, fundamental and optical gaps of many, diverse systems in excellent agreement with experiments, including organic and inorganic molecules and nanocrystals. We further demonstrate that the newly defined hybrid functional gives the correct alignment between the energy level of the exemplary TTF-TCNQ donor-acceptor system. DOE-BES: DE-FG02-06ER46262.
Jiang, Bin; Guo, Hua
2016-08-01
In search for an accurate description of the dissociative chemisorption of water on the Ni(111) surface, we report a new nine-dimensional potential energy surface (PES) based on a large number of density functional theory points using the RPBE functional. Seven-dimensional quantum dynamical calculations have been carried out on the RPBE PES, followed by site averaging and lattice effect corrections, yielding sticking probabilities that are compared with both the previous theoretical results based on a PW91 PES and experiment. It is shown that the RPBE functional increases the reaction barrier, but has otherwise a minor impact on the PES topography. Better agreement with experimental results is obtained with the new PES, but the agreement is still not quantitative. Possible sources of the remaining discrepancies are discussed. PMID:27436348
NASA Astrophysics Data System (ADS)
Deb, S.; Maitra, K.; Roychoudhuri, A.
1985-06-01
In the wake of the energy crisis, attempts are being made to develop a variety of energy conversion devices, such as solar cells. The single most important operational characteristic for a conversion element generating electricity is the V against I curve. Three points on this characteristic curve are of paramount importance, including the short-circuit, the open-circuit, and the maximum power point. The present paper has the objective to propose a new simple and accurate method of determining the maximum power point (Vm, Im) of the V against I characteristics, based on a geometrical interpretation. The method is general enough to be applicable to any energy conversion device having a nonlinear V against I characteristic. The paper provides also a method for determining the fill factor (FF), the series resistance (Rs), and the diode ideality factor (A) from a single set of connected observations.
NASA Astrophysics Data System (ADS)
Behmanesh, Iman; Moaveni, Babak
2016-07-01
This paper presents a Hierarchical Bayesian model updating framework to account for the effects of ambient temperature and excitation amplitude. The proposed approach is applied for model calibration, response prediction and damage identification of a footbridge under changing environmental/ambient conditions. The concrete Young's modulus of the footbridge deck is the considered updating structural parameter with its mean and variance modeled as functions of temperature and excitation amplitude. The identified modal parameters over 27 months of continuous monitoring of the footbridge are used to calibrate the updating parameters. One of the objectives of this study is to show that by increasing the levels of information in the updating process, the posterior variation of the updating structural parameter (concrete Young's modulus) is reduced. To this end, the calibration is performed at three information levels using (1) the identified modal parameters, (2) modal parameters and ambient temperatures, and (3) modal parameters, ambient temperatures, and excitation amplitudes. The calibrated model is then validated by comparing the model-predicted natural frequencies and those identified from measured data after deliberate change to the structural mass. It is shown that accounting for modeling error uncertainties is crucial for reliable response prediction, and accounting only the estimated variability of the updating structural parameter is not sufficient for accurate response predictions. Finally, the calibrated model is used for damage identification of the footbridge.
Yamagata, Koichi; Yamanishi, Ayako; Kokubu, Chikara; Takeda, Junji; Sese, Jun
2016-01-01
An important challenge in cancer genomics is precise detection of structural variations (SVs) by high-throughput short-read sequencing, which is hampered by the high false discovery rates of existing analysis tools. Here, we propose an accurate SV detection method named COSMOS, which compares the statistics of the mapped read pairs in tumor samples with isogenic normal control samples in a distinct asymmetric manner. COSMOS also prioritizes the candidate SVs using strand-specific read-depth information. Performance tests on modeled tumor genomes revealed that COSMOS outperformed existing methods in terms of F-measure. We also applied COSMOS to an experimental mouse cell-based model, in which SVs were induced by genome engineering and gamma-ray irradiation, followed by polymerase chain reaction-based confirmation. The precision of COSMOS was 84.5%, while the next best existing method was 70.4%. Moreover, the sensitivity of COSMOS was the highest, indicating that COSMOS has great potential for cancer genome analysis. PMID:26833260
Structural stability augmentation system design using BODEDIRECT: A quick and accurate approach
NASA Technical Reports Server (NTRS)
Goslin, T. J.; Ho, J. K.
1989-01-01
A methodology is presented for a modal suppression control law design using flight test data instead of mathematical models to obtain the required gain and phase information about the flexible airplane. This approach is referred to as BODEDIRECT. The purpose of the BODEDIRECT program is to provide a method of analyzing the modal phase relationships measured directly from the airplane. These measurements can be achieved with a frequency sweep at the control surface input while measuring the outputs of interest. The measured Bode-models can be used directly for analysis in the frequency domain, and for control law design. Besides providing a more accurate representation for the system inputs and outputs of interest, this method is quick and relatively inexpensive. To date, the BODEDIRECT program has been tested and verified for computational integrity. Its capabilities include calculation of series, parallel and loop closure connections between Bode-model representations. System PSD, together with gain and phase margins of stability may be calculated for successive loop closures of multi-input/multi-output systems. Current plans include extensive flight testing to obtain a Bode-model representation of a commercial aircraft for design of a structural stability augmentation system.
Fang, Tao; Li, Wei; Gu, Fangwei; Li, Shuhua
2015-01-13
We extend the generalized energy-based fragmentation (GEBF) approach to molecular crystals under periodic boundary conditions (PBC), and we demonstrate the performance of the method for a variety of molecular crystals. With this approach, the lattice energy of a molecular crystal can be obtained from the energies of a series of embedded subsystems, which can be computed with existing advanced molecular quantum chemistry methods. The use of the field compensation method allows the method to take long-range electrostatic interaction of the infinite crystal environment into account and make the method almost translationally invariant. The computational cost of the present method scales linearly with the number of molecules in the unit cell. Illustrative applications demonstrate that the PBC-GEBF method with explicitly correlated quantum chemistry methods is capable of providing accurate descriptions on the lattice energies and structures for various types of molecular crystals. In addition, this approach can be employed to quantify the contributions of various intermolecular interactions to the theoretical lattice energy. Such qualitative understanding is very useful for rational design of molecular crystals. PMID:26574207
Sharma, Virag; Elghafari, Anas; Hiller, Michael
2016-01-01
Identifying coding genes is an essential step in genome annotation. Here, we utilize existing whole genome alignments to detect conserved coding exons and then map gene annotations from one genome to many aligned genomes. We show that genome alignments contain thousands of spurious frameshifts and splice site mutations in exons that are truly conserved. To overcome these limitations, we have developed CESAR (Coding Exon-Structure Aware Realigner) that realigns coding exons, while considering reading frame and splice sites of each exon. CESAR effectively avoids spurious frameshifts in conserved genes and detects 91% of shifted splice sites. This results in the identification of thousands of additional conserved exons and 99% of the exons that lack inactivating mutations match real exons. Finally, to demonstrate the potential of using CESAR for comparative gene annotation, we applied it to 188 788 exons of 19 865 human genes to annotate human genes in 99 other vertebrates. These comparative gene annotations are available as a resource (http://bds.mpi-cbg.de/hillerlab/CESAR/). CESAR (https://github.com/hillerlab/CESAR/) can readily be applied to other alignments to accurately annotate coding genes in many other vertebrate and invertebrate genomes. PMID:27016733
Straightforward and accurate technique for post-coupler stabilization in drift tube linac structures
NASA Astrophysics Data System (ADS)
Khalvati, Mohammad Reza; Ramberger, Suitbert
2016-04-01
The axial electric field of Alvarez drift tube linacs (DTLs) is known to be susceptible to variations due to static and dynamic effects like manufacturing tolerances and beam loading. Post-couplers are used to stabilize the accelerating fields of DTLs against tuning errors. Tilt sensitivity and its slope have been introduced as measures for the stability right from the invention of post-couplers but since then the actual stabilization has mostly been done by tedious iteration. In the present article, the local tilt-sensitivity slope TSn' is established as the principal measure for stabilization instead of tilt sensitivity or some visual slope, and its significance is developed on the basis of an equivalent-circuit diagram of the DTL. Experimental and 3D simulation results are used to analyze its behavior and to define a technique for stabilization that allows finding the best post-coupler settings with just four tilt-sensitivity measurements. CERN's Linac4 DTL Tank 2 and Tank 3 have been stabilized successfully using this technique. The final tilt-sensitivity error has been reduced from ±100 %/MHz down to ±3 %/MHz for Tank 2 and down to ±1 %/MHz for Tank 3. Finally, an accurate procedure for tuning the structure using slug tuners is discussed.
Sharma, Virag; Elghafari, Anas; Hiller, Michael
2016-06-20
Identifying coding genes is an essential step in genome annotation. Here, we utilize existing whole genome alignments to detect conserved coding exons and then map gene annotations from one genome to many aligned genomes. We show that genome alignments contain thousands of spurious frameshifts and splice site mutations in exons that are truly conserved. To overcome these limitations, we have developed CESAR (Coding Exon-Structure Aware Realigner) that realigns coding exons, while considering reading frame and splice sites of each exon. CESAR effectively avoids spurious frameshifts in conserved genes and detects 91% of shifted splice sites. This results in the identification of thousands of additional conserved exons and 99% of the exons that lack inactivating mutations match real exons. Finally, to demonstrate the potential of using CESAR for comparative gene annotation, we applied it to 188 788 exons of 19 865 human genes to annotate human genes in 99 other vertebrates. These comparative gene annotations are available as a resource (http://bds.mpi-cbg.de/hillerlab/CESAR/). CESAR (https://github.com/hillerlab/CESAR/) can readily be applied to other alignments to accurately annotate coding genes in many other vertebrate and invertebrate genomes. PMID:27016733
Yamagata, Koichi; Yamanishi, Ayako; Kokubu, Chikara; Takeda, Junji; Sese, Jun
2016-05-01
An important challenge in cancer genomics is precise detection of structural variations (SVs) by high-throughput short-read sequencing, which is hampered by the high false discovery rates of existing analysis tools. Here, we propose an accurate SV detection method named COSMOS, which compares the statistics of the mapped read pairs in tumor samples with isogenic normal control samples in a distinct asymmetric manner. COSMOS also prioritizes the candidate SVs using strand-specific read-depth information. Performance tests on modeled tumor genomes revealed that COSMOS outperformed existing methods in terms of F-measure. We also applied COSMOS to an experimental mouse cell-based model, in which SVs were induced by genome engineering and gamma-ray irradiation, followed by polymerase chain reaction-based confirmation. The precision of COSMOS was 84.5%, while the next best existing method was 70.4%. Moreover, the sensitivity of COSMOS was the highest, indicating that COSMOS has great potential for cancer genome analysis. PMID:26833260
NASA Astrophysics Data System (ADS)
Marin, Andrew T.; Musselman, Kevin P.; MacManus-Driscoll, Judith L.
2013-04-01
This work shows that when a Schottky barrier is present in a photovoltaic device, such as in a device with an ITO/ZnO contact, equivalent circuit analysis must be performed with admittance spectroscopy to accurately determine the pn junction interface recombination parameters (i.e., capture cross section and density of trap states). Without equivalent circuit analysis, a Schottky barrier can produce an error of ˜4-orders of magnitude in the capture cross section and ˜50% error in the measured density of trap states. Using a solution processed ZnO/Cu2O photovoltaic test system, we apply our analysis to clearly separate the contributions of interface states at the pn junction from the Schottky barrier at the ITO/ZnO contact so that the interface state recombination parameters can be accurately characterized. This work is widely applicable to the multitude of photovoltaic devices, which use ZnO adjacent to ITO.
NASA Astrophysics Data System (ADS)
Zhang, Lili
This work aims to improve the capability of accurate information extraction from high-dimensional data, with a specific neural learning paradigm, the Self-Organizing Map (SOM). The SOM is an unsupervised learning algorithm that can faithfully sense the manifold structure and support supervised learning of relevant information from the data. Yet open problems regarding SOM learning exist. We focus on the following two issues. (1) Evaluation of topology preservation. Topology preservation is essential for SOMs in faithful representation of manifold structure. However, in reality, topology violations are not unusual, especially when the data have complicated structure. Measures capable of accurately quantifying and informatively expressing topology violations are lacking. One contribution of this work is a new measure, the Weighted Differential Topographic Function (WDTF), which differentiates an existing measure, the Topographic Function (TF), and incorporates detailed data distribution as an importance weighting of violations to distinguish severe violations from insignificant ones. Another contribution is an interactive visual tool, TopoView, which facilitates the visual inspection of violations on the SOM lattice. We show the effectiveness of the combined use of the WDTF and TopoView through a simple two-dimensional data set and two hyperspectral images. (2) Learning multiple latent variables from high-dimensional data. We use an existing two-layer SOM-hybrid supervised architecture, which captures the manifold structure in its SOM hidden layer, and then, uses its output layer to perform the supervised learning of latent variables. In the customary way, the output layer only uses the strongest output of the SOM neurons. This severely limits the learning capability. We allow multiple, k, strongest responses of the SOM neurons for the supervised learning. Moreover, the fact that different latent variables can be best learned with different values of k motivates a
Kostylev, Maxim; Wilson, David
2014-01-01
Lignocellulosic biomass is a potential source of renewable, low-carbon-footprint liquid fuels. Biomass recalcitrance and enzyme cost are key challenges associated with the large-scale production of cellulosic fuel. Kinetic modeling of enzymatic cellulose digestion has been complicated by the heterogeneous nature of the substrate and by the fact that a true steady state cannot be attained. We present a two-parameter kinetic model based on the Michaelis-Menten scheme (Michaelis L and Menten ML. (1913) Biochem Z 49:333–369), but with a time-dependent activity coefficient analogous to fractal-like kinetics formulated by Kopelman (Kopelman R. (1988) Science 241:1620–1626). We provide a mathematical derivation and experimental support to show that one of the parameters is a total activity coefficient and the other is an intrinsic constant that reflects the ability of the cellulases to overcome substrate recalcitrance. The model is applicable to individual cellulases and their mixtures at low-to-medium enzyme loads. Using biomass degrading enzymes from a cellulolytic bacterium Thermobifida fusca we show that the model can be used for mechanistic studies of enzymatic cellulose digestion. We also demonstrate that it applies to the crude supernatant of the widely studied cellulolytic fungus Trichoderma reesei and can thus be used to compare cellulases from different organisms. The two parameters may serve a similar role to Vmax, KM, and kcat in classical kinetics. A similar approach may be applicable to other enzymes with heterogeneous substrates and where a steady state is not achievable. PMID:23837567
Nicolotti, Orazio; Gillet, Valerie J; Fleming, Peter J; Green, Darren V S
2002-11-01
Deriving quantitative structure-activity relationship (QSAR) models that are accurate, reliable, and easily interpretable is a difficult task. In this study, two new methods have been developed that aim to find useful QSAR models that represent an appropriate balance between model accuracy and complexity. Both methods are based on genetic programming (GP). The first method, referred to as genetic QSAR (or GPQSAR), uses a penalty function to control model complexity. GPQSAR is designed to derive a single linear model that represents an appropriate balance between the variance and the number of descriptors selected for the model. The second method, referred to as multiobjective genetic QSAR (MoQSAR), is based on multiobjective GP and represents a new way of thinking of QSAR. Specifically, QSAR is considered as a multiobjective optimization problem that comprises a number of competitive objectives. Typical objectives include model fitting, the total number of terms, and the occurrence of nonlinear terms. MoQSAR results in a family of equivalent QSAR models where each QSAR represents a different tradeoff in the objectives. A practical consideration often overlooked in QSAR studies is the need for the model to promote an understanding of the biochemical response under investigation. To accomplish this, chemically intuitive descriptors are needed but do not always give rise to statistically robust models. This problem is addressed by the addition of a further objective, called chemical desirability, that aims to reward models that consist of descriptors that are easily interpretable by chemists. GPQSAR and MoQSAR have been tested on various data sets including the Selwood data set and two different solubility data sets. The study demonstrates that the MoQSAR method is able to find models that are at least as good as models derived using standard statistical approaches and also yields models that allow a medicinal chemist to trade statistical robustness for chemical
Lee, Yu Ran; Kang, Do Won; Kim, Hong Lae E-mail: hlkim@kangwon.ac.kr; Kwon, Chan Ho E-mail: hlkim@kangwon.ac.kr
2014-11-07
Ionization energies and cationic structures of pyridine were intensively investigated utilizing one-photon mass-analyzed threshold ionization (MATI) spectroscopy with vacuum ultraviolet radiation generated by four-wave difference frequency mixing in Kr. The present one-photon high-resolution MATI spectrum of pyridine demonstrated a much finer and richer vibrational structure than that of the previously reported two-photon MATI spectrum. From the MATI spectrum and photoionization efficiency curve, the accurate ionization energy of the ionic ground state of pyridine was confidently determined to be 73 570 ± 6 cm{sup −1} (9.1215 ± 0.0007 eV). The observed spectrum was almost completely assigned by utilizing Franck-Condon factors and vibrational frequencies calculated through adjustments of the geometrical parameters of cationic pyridine at the B3LYP/cc-pVTZ level. A unique feature unveiled through rigorous analysis was the prominent progression of the 10 vibrational mode, which corresponds to in-plane ring bending, and the combination of other totally symmetric fundamentals with the ring bending overtones, which contribute to the geometrical change upon ionization. Notably, the remaining peaks originate from the upper electronic state ({sup 2}A{sub 2}), as predicted by high-resolution photoelectron spectroscopy studies and symmetry-adapted cluster configuration interaction calculations. Based on the quantitatively good agreement between the experimental and calculated results, it was concluded that upon ionization the pyridine cation in the ground electronic state should have a planar structure of C{sub 2v} symmetry through the C-N axis.
NASA Astrophysics Data System (ADS)
Chiu, Yung-Chia
2014-12-01
Parameter structure identification is formulated in terms of solving an inverse problem, which allows for a determination of an appropriate level of parameter structure complexity, and the identification of its pattern and the associated parameter values. With the increasing complexity of parameter structure identification in groundwater modeling, demand for robust, fast, and accurate optimizers is on the rise among researchers from groundwater hydrology fields. A novel global optimizer, differential evolution (DE), has been proposed to solve the parameter-structure-identification problem. The Voronoi tessellation is adopted for the automatic parameterization. The stepwise regression method and the error covariance matrix are used to determine the optimal structure complexity. Numerical experiments with a continuous hydraulic conductivity distribution are conducted to demonstrate the proposed methodology. The results indicate that the DE can identify the global optimum effectively and efficiently. A sensitivity analysis of the control parameters and mutation schemes implemented in the DE is employed to examine their influence on the objective function. The comparison between DE and genetic algorithm shows the advantage of DE in terms of robustness and efficiency. The proposed methodology is also applied to a real groundwater system, Pingtung Plain in Taiwan, and the properties of aquifers are successfully identified.
NASA Astrophysics Data System (ADS)
Roy Choudhury, Raja; Roy Choudhury, Arundhati; Kanti Ghose, Mrinal
2013-01-01
A semi-analytical model with three optimizing parameters and a novel non-Gaussian function as the fundamental modal field solution has been proposed to arrive at an accurate solution to predict various propagation parameters of graded-index fibers with less computational burden than numerical methods. In our semi analytical formulation the optimization of core parameter U which is usually uncertain, noisy or even discontinuous, is being calculated by Nelder-Mead method of nonlinear unconstrained minimizations as it is an efficient and compact direct search method and does not need any derivative information. Three optimizing parameters are included in the formulation of fundamental modal field of an optical fiber to make it more flexible and accurate than other available approximations. Employing variational technique, Petermann I and II spot sizes have been evaluated for triangular and trapezoidal-index fibers with the proposed fundamental modal field. It has been demonstrated that, the results of the proposed solution identically match with the numerical results over a wide range of normalized frequencies. This approximation can also be used in the study of doped and nonlinear fiber amplifier.
Effective Parameters on Seismic Design of Rectangular Underground Structures
Amiri, G. Ghodrati; Maddah, N.; Mohebi, B.
2008-07-08
Underground structures are a significant part of the transportation in the modern society and in the seismic zones should withstand against both seismic and static loadings. Embedded structures should conform to ground deformations during the earthquake but almost exact evaluation of structure to ground distortion is critical. Several two-dimensional finite difference models are used to find effective parameters on racking ratio (structure to ground distortion) including flexibility ratio, various cross sections, embedment depth, and Poisson's ratio of soil. Results show that influence of different cross sections, by themselves is negligible but embedment depth in addition to flexibility ratio and Poisson's ratio is known as a consequential parameter. A comparison with pseudo-static method (simplified frame analysis) is also performed. The results show that for a stiffer structure than soil, racking ratio decreases as the depth of burial decreases; on the other hand, shallow and flexible structures can suffer greater distortion than deeper ones up to 30 percents.
NASA Astrophysics Data System (ADS)
Ryu, Jaiyoung; Hu, Xiao; Shadden, Shawn C.
2014-11-01
The cerebral circulation is unique in its ability to maintain blood flow to the brain under widely varying physiologic conditions. Incorporating this autoregulatory response is critical to cerebral blood flow modeling, as well as investigations into pathological conditions. We discuss a one-dimensional nonlinear model of blood flow in the cerebral arteries that includes coupling of autoregulatory lumped parameter networks. The model is tested to reproduce a common clinical test to assess autoregulatory function - the carotid artery compression test. The change in the flow velocity at the middle cerebral artery (MCA) during carotid compression and release demonstrated strong agreement with published measurements. The model is then used to investigate vasospasm of the MCA, a common clinical concern following subarachnoid hemorrhage. Vasospasm was modeled by prescribing vessel area reduction in the middle portion of the MCA. Our model showed similar increases in velocity for moderate vasospasms, however, for serious vasospasm (~ 90% area reduction), the blood flow velocity demonstrated decrease due to blood flow rerouting. This demonstrates a potentially important phenomenon, which otherwise would lead to false-negative decisions on clinical vasospasm if not properly anticipated.
Sandala, Gregory M.; Hopmann, Kathrin H.; Ghosh, Abhik
2011-01-01
structure. Significant improvements to the isomer shift calibrations are obtained for B3LYP and B3LYP* when geometries obtained with the OLYP functional are used. In addition, greatly improved performance of these functionals is found if the complete test set is grouped separately into Fe–NO and Fe–S complexes. Calibration fits including only Fe–NO complexes are found to be excellent, while those containing the non-nitrosyl Fe–S complexes alone are found to demonstrate less accurate correlations. Similar trends are also found with OLYP, OPBE, PW91, and BP86. Correlations between experimental and calculated QSs were also investigated. Generally, universal and separate Fe–NO and Fe–S fit parameters obtained to determine QSs are found to be of good to excellent quality for every density functional examined, especially if [Fe4(NO)4(μ3-S)4]− is removed from the test set. PMID:22039359
Uncertainty in dual permeability model parameters for structured soils
Arora, B.; Mohanty, B. P.; McGuire, J. T.
2013-01-01
Successful application of dual permeability models (DPM) to predict contaminant transport is contingent upon measured or inversely estimated soil hydraulic and solute transport parameters. The difficulty in unique identification of parameters for the additional macropore- and matrix-macropore interface regions, and knowledge about requisite experimental data for DPM has not been resolved to date. Therefore, this study quantifies uncertainty in dual permeability model parameters of experimental soil columns with different macropore distributions (single macropore, and low- and high-density multiple macropores). Uncertainty evaluation is conducted using adaptive Markov chain Monte Carlo (AMCMC) and conventional Metropolis-Hastings (MH) algorithms while assuming 10 out of 17 parameters to be uncertain or random. Results indicate that AMCMC resolves parameter correlations and exhibits fast convergence for all DPM parameters while MH displays large posterior correlations for various parameters. This study demonstrates that the choice of parameter sampling algorithms is paramount in obtaining unique DPM parameters when information on covariance structure is lacking, or else additional information on parameter correlations must be supplied to resolve the problem of equifinality of DPM parameters. This study also highlights the placement and significance of matrix-macropore interface in flow experiments of soil columns with different macropore densities. Histograms for certain soil hydraulic parameters display tri-modal characteristics implying that macropores are drained first followed by the interface region and then by pores of the matrix domain in drainage experiments. Results indicate that hydraulic properties and behavior of the matrix-macropore interface is not only a function of saturated hydraulic conductivity of the macroporematrix interface (Ksa) and macropore tortuosity (lf) but also of other parameters of the matrix and macropore domains. PMID:24478531
Crystalline structure and symmetry dependence of acoustic nonlinearity parameters
NASA Technical Reports Server (NTRS)
Cantrell, John H.
1994-01-01
A quantitative measure of elastic wave nonlinearity in crystals is provided by the acoustic nonlinearity parameters. The nonlinearity parameters are defined for arbitrary propagation modes for solids of arbitrary crystalline symmetry and are determined along the pure mode propagation directions for 33 crystals of cubic symmetry from data reported in the literature. The magnitudes of the nonlinearity parameters are found to exhibit a strong dependence on the crystalline structure and symmetries associated with the modal direction in the solid. Calculations based on the Born-Mayer potential for crystals having a dominant repulsive contribution to the elastic constants from the interatomic pair potential suggest that the origin of the structure dependence is associated with the shape rather than the strength of the potential. Considerations based on variations in crystal symmetry during loading along pure mode propagation directions of face-centered-cubic solids provide a qualitative explanation for the dependence of the acoustic nonlinearity parameters on modal direction.
Baker, Christopher M.; Lopes, Pedro E. M.; Zhu, Xiao; Roux, Benoît; MacKerell, Alexander D.
2010-01-01
Lennard-Jones (LJ) parameters for a variety of model compounds have previously been optimized within the CHARMM Drude polarizable force field to reproduce accurately pure liquid phase thermodynamic properties as well as additional target data. While the polarizable force field resulting from this optimization procedure has been shown to satisfactorily reproduce a wide range of experimental reference data across numerous series of small molecules, a slight but systematic overestimate of the hydration free energies has also been noted. Here, the reproduction of experimental hydration free energies is greatly improved by the introduction of pair-specific LJ parameters between solute heavy atoms and water oxygen atoms that override the standard LJ parameters obtained from combining rules. The changes are small and a systematic protocol is developed for the optimization of pair-specific LJ parameters and applied to the development of pair-specific LJ parameters for alkanes, alcohols and ethers. The resulting parameters not only yield hydration free energies in good agreement with experimental values, but also provide a framework upon which other pair-specific LJ parameters can be added as new compounds are parametrized within the CHARMM Drude polarizable force field. Detailed analysis of the contributions to the hydration free energies reveals that the dispersion interaction is the main source of the systematic errors in the hydration free energies. This information suggests that the systematic error may result from problems with the LJ combining rules and is combined with analysis of the pair-specific LJ parameters obtained in this work to identify a preliminary improved combining rule. PMID:20401166
Identification of uncertain structural parameters in flexible boosters
NASA Astrophysics Data System (ADS)
Sawai, Shujiro; Kawaguchi, Jun'ichiro; Matsuo, Hiroki
A novel algorithm for identification of unknown structural parameters in flexible boosters is presented here. In this method, first of all coefficient matrices in recursive form are estimated based on sensor outputs, which is followed by identifying structural parameters utilizing matrix relations between those and system description form. This algorithm is relatively robust and promising, since neither modal coordinates nor state reconstruction is required, in which tremendous amount of efforts have been to be concentrated to. Numerical discussion here was carried out and it demonstrated its practical effectiveness.
NASA Astrophysics Data System (ADS)
Fang, J.; Sheng, L.; Li, D.; Zhao, J.; Li, Sh.; Qin, W.; Fan, Y.; Zheng, Q. L.; Zhang, W.
A novel High Temperature Superconductor Linear Induction Motor (HTS LIM) is researched in this paper. Since the critical current and the electromagnetic force of the motor are determined mainly by the primary slot leakage flux, the main magnetic flux and eddy current respectively, in order to research the influence of structural parameters and operating parameters on electromagnetic properties of HTS LIM, the motor was analyzed by 2D transient Finite Element Method (FEM). The properties of the motor, such as the maximum slot leakage flux density, motor thrust, motor vertical force and critical current are analyzed with different structural parameters and operating parameters. In addition, an experimental investigation was carried out on prototype HTS motor. Electrical parameters were deduced from these tests and also compared with the analysis results from FEM. AC losses of one HTS coil in the motor were measured and AC losses of all HTS coils in HTS LIM were estimated. The results in this paper could provide reference for the design and research on the HTS LIM.
Precision and Accuracy Parameters in Structured Light 3-D Scanning
NASA Astrophysics Data System (ADS)
Eiríksson, E. R.; Wilm, J.; Pedersen, D. B.; Aanæs, H.
2016-04-01
Structured light systems are popular in part because they can be constructed from off-the-shelf low cost components. In this paper we quantitatively show how common design parameters affect precision and accuracy in such systems, supplying a much needed guide for practitioners. Our quantitative measure is the established VDI/VDE 2634 (Part 2) guideline using precision made calibration artifacts. Experiments are performed on our own structured light setup, consisting of two cameras and a projector. We place our focus on the influence of calibration design parameters, the calibration procedure and encoding strategy and present our findings. Finally, we compare our setup to a state of the art metrology grade commercial scanner. Our results show that comparable, and in some cases better, results can be obtained using the parameter settings determined in this study.
Using a quadratic parameter sinusoid model to characterize the structure of EEG sleep spindles
Palliyali, Abdul J.; Ahmed, Mohammad N.; Ahmed, Beena
2015-01-01
Sleep spindles are essentially non-stationary signals that display time and frequency-varying characteristics within their envelope, which makes it difficult to accurately identify its instantaneous frequency and amplitude. To allow a better parameterization of the structure of spindle, we propose modeling spindles using a Quadratic Parameter Sinusoid (QPS). The QPS is well suited to model spindle activity as it utilizes a quadratic representation to capture the inherent duration and frequency variations within spindles. The effectiveness of our proposed model and estimation technique was quantitatively evaluated in parameter determination experiments using simulated spindle-like signals and real spindles in the presence of background EEG. We used the QPS parameters to predict the energy and frequency of spindles with a mean accuracy of 92.34 and 97.73% respectively. We also show that the QPS parameters provide a quantification of the amplitude and frequency variations occurring within sleep spindles that can be observed visually and related to their characteristic “waxing and waning” shape. We analyze the variations in the parameters values to present how they can be used to understand the inter- and intra-participant variations in spindle structure. Finally, we present a comparison of the QPS parameters of spindles and non-spindles, which shows a substantial difference in parameter values between the two classes. PMID:25999833
Multiplicity Control in Structural Equation Modeling: Incorporating Parameter Dependencies
ERIC Educational Resources Information Center
Smith, Carrie E.; Cribbie, Robert A.
2013-01-01
When structural equation modeling (SEM) analyses are conducted, significance tests for all important model relationships (parameters including factor loadings, covariances, etc.) are typically conducted at a specified nominal Type I error rate ([alpha]). Despite the fact that many significance tests are often conducted in SEM, rarely is…
Bialasiewicz, J.T.
1995-06-01
The goal of this research is to develop advanced system identification techniques that can be used to accurately measure the frequency response functions of a wind-turbine structure immersed in wind noise. To allow for accurate identification, the authors have developed a special test signal called the Pseudo-Random Binary Sequence (PRBS). The Matlab program that generates this signal allows the user to interactively tailor its parameters for the frequency range of interest based on the response of the wind turbine under test. By controlling NREL`s Mobile Hydraulic Shaker System, which is attached to the wind turbine structure, the PRBS signal produces the wide-band excitation necessary to perform system identification in the presence of wind noise. The techniques presented here will enable researchers to obtain modal parameters from an operating wind turbine, including frequencies, damping coefficients, and mode shapes. More importantly, the algorithms they have developed and tested (so far using input-output data from a simulated structure) permit state-space representation of the system under test, particularly the modal state space representation. This is the only system description that reveals the internal behavior the system, such as the interaction between the physical parameters, and which, in contrast to transfer functions, is valid for non-zero initial conditions.
PROMALS3D web server for accurate multiple protein sequence and structure alignments.
Pei, Jimin; Tang, Ming; Grishin, Nick V
2008-07-01
Multiple sequence alignments are essential in computational sequence and structural analysis, with applications in homology detection, structure modeling, function prediction and phylogenetic analysis. We report PROMALS3D web server for constructing alignments for multiple protein sequences and/or structures using information from available 3D structures, database homologs and predicted secondary structures. PROMALS3D shows higher alignment accuracy than a number of other advanced methods. Input of PROMALS3D web server can be FASTA format protein sequences, PDB format protein structures and/or user-defined alignment constraints. The output page provides alignments with several formats, including a colored alignment augmented with useful information about sequence grouping, predicted secondary structures and consensus sequences. Intermediate results of sequence and structural database searches are also available. The PROMALS3D web server is available at: http://prodata.swmed.edu/promals3d/. PMID:18503087
Asmadi, Aldi; Neumann, Marcus A; Kendrick, John; Girard, Pascale; Perrin, Marc-Antoine; Leusen, Frank J J
2009-12-24
In the 2007 blind test of crystal structure prediction hosted by the Cambridge Crystallographic Data Centre (CCDC), a hybrid DFT/MM method correctly ranked each of the four experimental structures as having the lowest lattice energy of all the crystal structures predicted for each molecule. The work presented here further validates this hybrid method by optimizing the crystal structures (experimental and submitted) of the first three CCDC blind tests held in 1999, 2001, and 2004. Except for the crystal structures of compound IX, all structures were reminimized and ranked according to their lattice energies. The hybrid method computes the lattice energy of a crystal structure as the sum of the DFT total energy and a van der Waals (dispersion) energy correction. Considering all four blind tests, the crystal structure with the lowest lattice energy corresponds to the experimentally observed structure for 12 out of 14 molecules. Moreover, good geometrical agreement is observed between the structures determined by the hybrid method and those measured experimentally. In comparison with the correct submissions made by the blind test participants, all hybrid optimized crystal structures (apart from compound II) have the smallest calculated root mean squared deviations from the experimentally observed structures. It is predicted that a new polymorph of compound V exists under pressure. PMID:19950907
Seeking: Accurate Measurement Techniques for Deep-Bone Density and Structure
NASA Technical Reports Server (NTRS)
Sibonga, Jean
2009-01-01
We are seeking a clinically-useful technology with enough sensitivity to assess the microstructure of "spongy" bone that is found in the marrow cavities of whole bones. However, this technology must be for skeletal sites surrounded by layers of soft tissues, such as the spine and the hip. Soft tissue interferes with conventional imaging and using a more accessible area -- for example, the wrist or the ankle of limbs-- as a proxy for the less accessible skeletal regions, will not be accurate. A non-radioactive technology is strongly preferred.
Wheat streak mosaic virus-Structural parameters for a Potyvirus
Parker, Lauren; Kendall, Amy; Berger, P.H.; Shiel, P.J.; Stubbs, Gerald . E-mail: gerald.stubbs@vanderbilt.edu
2005-09-15
Wheat streak mosaic virus is a Tritimovirus, a member of the Potyviridae family, which includes the very large Potyvirus genus. We have examined wheat streak mosaic virus by electron microscopy and fiber diffraction from partially oriented sols, and analyzed the results to estimate the symmetry and structural parameters of the viral helix. The virions have an apparent radius of 63 {+-} 5 A. The viral helix has a pitch of 33.4 A {+-} 0.6 A. There appear to be 6.9 subunits per turn of the helix, although we cannot completely eliminate values of 5.9 or 7.9 for this parameter.
Order parameter in complex dipolar structures: Microscopic modeling
NASA Astrophysics Data System (ADS)
Prosandeev, S.; Bellaiche, L.
2008-02-01
Microscopic models have been used to reveal the existence of an order parameter that is associated with many complex dipolar structures in magnets and ferroelectrics. This order parameter involves a double cross product of the local dipoles with their positions. It provides a measure of subtle microscopic features, such as the helicity of the two domains inherent to onion states, curvature of the dipolar pattern in flower states, or characteristics of sets of vortices with opposite chirality (e.g., distance between the vortex centers and/or the magnitude of their local dipoles).
The role of structural parameters in DNA cyclization
Alexandrov, Ludmil B.; Bishop, Alan R.; Rasmussen, Kim O.; Alexandrov, Boian S.
2016-02-04
The intrinsic bendability of DNA plays an important role with relevance for myriad of essential cellular mechanisms. The flexibility of a DNA fragment can be experimentally and computationally examined by its propensity for cyclization, quantified by the Jacobson-Stockmayer J factor. In this paper, we use a well-established coarse-grained three-dimensional model of DNA and seven distinct sets of experimentally and computationally derived conformational parameters of the double helix to evaluate the role of structural parameters in calculating DNA cyclization.
PDBparam: Online Resource for Computing Structural Parameters of Proteins
Nagarajan, R.; Archana, A.; Thangakani, A. Mary; Jemimah, S.; Velmurugan, D.; Gromiha, M. Michael
2016-01-01
Understanding the structure–function relationship in proteins is a longstanding goal in molecular and computational biology. The development of structure-based parameters has helped to relate the structure with the function of a protein. Although several structural features have been reported in the literature, no single server can calculate a wide-ranging set of structure-based features from protein three-dimensional structures. In this work, we have developed a web-based tool, PDBparam, for computing more than 50 structure-based features for any given protein structure. These features are classified into four major categories: (i) interresidue interactions, which include short-, medium-, and long-range interactions, contact order, long-range order, total contact distance, contact number, and multiple contact index, (ii) secondary structure propensities such as α-helical propensity, β-sheet propensity, and propensity of amino acids to exist at various positions of α-helix and amino acid compositions in high B-value regions, (iii) physicochemical properties containing ionic interactions, hydrogen bond interactions, hydrophobic interactions, disulfide interactions, aromatic interactions, surrounding hydrophobicity, and buriedness, and (iv) identification of binding site residues in protein–protein, protein–nucleic acid, and protein–ligand complexes. The server can be freely accessed at http://www.iitm.ac.in/bioinfo/pdbparam/. We suggest the use of PDBparam as an effective tool for analyzing protein structures. PMID:27330281
Structural and parameter uncertainty in Bayesian cost-effectiveness models
Jackson, Christopher H; Sharples, Linda D; Thompson, Simon G
2010-01-01
Health economic decision models are subject to various forms of uncertainty, including uncertainty about the parameters of the model and about the model structure. These uncertainties can be handled within a Bayesian framework, which also allows evidence from previous studies to be combined with the data. As an example, we consider a Markov model for assessing the cost-effectiveness of implantable cardioverter defibrillators. Using Markov chain Monte Carlo posterior simulation, uncertainty about the parameters of the model is formally incorporated in the estimates of expected cost and effectiveness. We extend these methods to include uncertainty about the choice between plausible model structures. This is accounted for by averaging the posterior distributions from the competing models using weights that are derived from the pseudo-marginal-likelihood and the deviance information criterion, which are measures of expected predictive utility. We also show how these cost-effectiveness calculations can be performed efficiently in the widely used software WinBUGS. PMID:20383261
Controllability of brushite structural parameters using an applied magnetic field.
Kuznetsov, V N; Yanovska, A A; Stanislavov, A S; Danilchenko, S N; Kalinkevich, A N; Sukhodub, L F
2016-03-01
The paper studies the influence of low intensity static magnetic field on brushite structural and microstructural parameters using the X-ray diffraction and the transmission electron microscopy. This effect was shown to have various influences on DCPD (Dicalcium Phosphate Dihydrate) structure depending on a magnetic field configuration or time of synthesis, which allows achieving controllability of the main properties of an obtained material. The influence of the magnetic field leads mostly to the decrease of crystallite sizes with no impact on the crystal lattice parameters. In (0 2 0) and (1 5 0) planes the growth of crystallite sizes is observed after 2 and 3 days of crystallization, respectively. The analysis of different contributions to peak broadening in [0 b 0] direction showed a similar trend for the crystallite sizes with the lower lattice microstrains after 2 days of synthesis. The effect similar to the preferred orientation was observed and classified with the Harris method. PMID:26706562
Zimmermann, Olav; Hansmann, Ulrich H E
2008-09-01
Constraint generation for 3d structure prediction and structure-based database searches benefit from fine-grained prediction of local structure. In this work, we present LOCUSTRA, a novel scheme for the multiclass prediction of local structure that uses two layers of support vector machines (SVM). Using a 16-letter structural alphabet from de Brevern et al. (Proteins: Struct., Funct., Bioinf. 2000, 41, 271-287), we assess its prediction ability for an independent test set of 222 proteins and compare our method to three-class secondary structure prediction and direct prediction of dihedral angles. The prediction accuracy is Q16=61.0% for the 16 classes of the structural alphabet and Q3=79.2% for a simple mapping to the three secondary classes helix, sheet, and coil. We achieve a mean phi(psi) error of 24.74 degrees (38.35 degrees) and a median RMSDA (root-mean-square deviation of the (dihedral) angles) per protein chain of 52.1 degrees. These results compare favorably with related approaches. The LOCUSTRA web server is freely available to researchers at http://www.fz-juelich.de/nic/cbb/service/service.php. PMID:18763837
Roses, A D
2016-02-01
Structural variants (SVs) include all insertions, deletions, and rearrangements in the genome, with several common types of nucleotide repeats including single sequence repeats, short tandem repeats, and insertion-deletion length variants. Polyallelic SVs provide highly informative markers for association studies with well-phenotyped cohorts. SVs can influence gene regulation by affecting epigenetics, transcription, splicing, and/or translation. Accurate assays of polyallelic SV loci are required to define the range and allele frequency of variable length alleles. PMID:26517180
Ergül, Özgür; Gürel, Levent
2013-03-01
Accurate electromagnetic modeling of complicated optical structures poses several challenges. Optical metamaterial and plasmonic structures are composed of multiple coexisting dielectric and/or conducting parts. Such composite structures may possess diverse values of conductivities and dielectric constants, including negative permittivity and permeability. Further challenges are the large sizes of the structures with respect to wavelength and the complexities of the geometries. In order to overcome these challenges and to achieve rigorous and efficient electromagnetic modeling of three-dimensional optical composite structures, we have developed a parallel implementation of the multilevel fast multipole algorithm (MLFMA). Precise formulation of composite structures is achieved with the so-called "electric and magnetic current combined-field integral equation." Surface integral equations are carefully discretized with piecewise linear basis functions, and the ensuing dense matrix equations are solved iteratively with parallel MLFMA. The hierarchical strategy is used for the efficient parallelization of MLFMA on distributed-memory architectures. In this paper, fast and accurate solutions of large-scale canonical and complicated real-life problems, such as optical metamaterials, discretized with tens of millions of unknowns are presented in order to demonstrate the capabilities of the proposed electromagnetic solver. PMID:23456127
Arumugam, P.; Ferreira, L. S.; Maglione, E.
2008-10-15
With a proper formalism for proton emission from triaxially deformed nuclei, we perform exact calculations of decay widths for the decays to ground and first excited 2{sup +} states in the daughter nucleus. Our results for rotational spectrum, decay width and fine structure in the case of the nucleus {sup 145}Tm lead for the first time to an accurate identification of triaxial deformation using proton emission. This work also puts in evidence the advantage of proton emission over the conventional probes to study nuclear structure at the proton drip-line.
NASA Astrophysics Data System (ADS)
Maeda, Chiaki; Tasaki, Satoko; Kirihara, Soshu
2011-05-01
Computer graphic models of bioscaffolds with four-coordinate lattice structures of solid rods in artificial bones were designed by using a computer aided design. The scaffold models composed of acryl resin with hydroxyapatite particles at 45vol. % were fabricated by using stereolithography of a computer aided manufacturing. After dewaxing and sintering heat treatment processes, the ceramics scaffold models with four-coordinate lattices and fine hydroxyapatite microstructures were obtained successfully. By using a computer aided analysis, it was found that bio-fluids could flow extensively inside the sintered scaffolds. This result shows that the lattice structures will realize appropriate bio-fluid circulations and promote regenerations of new bones.
Braun, Tatjana; Koehler Leman, Julia; Lange, Oliver F.
2015-01-01
Recent work has shown that the accuracy of ab initio structure prediction can be significantly improved by integrating evolutionary information in form of intra-protein residue-residue contacts. Following this seminal result, much effort is put into the improvement of contact predictions. However, there is also a substantial need to develop structure prediction protocols tailored to the type of restraints gained by contact predictions. Here, we present a structure prediction protocol that combines evolutionary information with the resolution-adapted structural recombination approach of Rosetta, called RASREC. Compared to the classic Rosetta ab initio protocol, RASREC achieves improved sampling, better convergence and higher robustness against incorrect distance restraints, making it the ideal sampling strategy for the stated problem. To demonstrate the accuracy of our protocol, we tested the approach on a diverse set of 28 globular proteins. Our method is able to converge for 26 out of the 28 targets and improves the average TM-score of the entire benchmark set from 0.55 to 0.72 when compared to the top ranked models obtained by the EVFold web server using identical contact predictions. Using a smaller benchmark, we furthermore show that the prediction accuracy of our method is only slightly reduced when the contact prediction accuracy is comparatively low. This observation is of special interest for protein sequences that only have a limited number of homologs. PMID:26713437
Abdelnour, Farras; Voss, Henning U.; Raj, Ashish
2014-01-01
The relationship between anatomic connectivity of large-scale brain networks and their functional connectivity is of immense importance and an area of active research. Previous attempts have required complex simulations which model the dynamics of each cortical region, and explore the coupling between regions as derived by anatomic connections. While much insight is gained from these non-linear simulations, they can be computationally taxing tools for predicting functional from anatomic connectivities. Little attention has been paid to linear models. Here we show that a properly designed linear model appears to be superior to previous non-linear approaches in capturing the brain’s long-range second order correlation structure that governs the relationship between anatomic and functional connectivities. We derive a linear network of brain dynamics based on graph diffusion, whereby the diffusing quantity undergoes a random walk on a graph. We test our model using subjects who underwent diffusion MRI and resting state fMRI. The network diffusion model applied to the structural networks largely predicts the correlation structures derived from their fMRI data, to a greater extent than other approaches. The utility of the proposed approach is that it can routinely be used to infer functional correlation from anatomic connectivity. And since it is linear, anatomic connectivity can also be inferred from functional data. The success of our model confirms the linearity of ensemble average signals in the brain, and implies that their long-range correlation structure may percolate within the brain via purely mechanistic processes enacted on its structural connectivity pathways. PMID:24384152
Interplay between Peptide Bond Geometrical Parameters in Nonglobular Structural Contexts
Esposito, Luciana; De Simone, Alfonso; Vitagliano, Luigi
2013-01-01
Several investigations performed in the last two decades have unveiled that geometrical parameters of protein backbone show a remarkable variability. Although these studies have provided interesting insights into one of the basic aspects of protein structure, they have been conducted on globular and water-soluble proteins. We report here a detailed analysis of backbone geometrical parameters in nonglobular proteins/peptides. We considered membrane proteins and two distinct fibrous systems (amyloid-forming and collagen-like peptides). Present data show that in these systems the local conformation plays a major role in dictating the amplitude of the bond angle N-Cα-C and the propensity of the peptide bond to adopt planar/nonplanar states. Since the trends detected here are in line with the concept of the mutual influence of local geometry and conformation previously established for globular and water-soluble proteins, our analysis demonstrates that the interplay of backbone geometrical parameters is an intrinsic and general property of protein/peptide structures that is preserved also in nonglobular contexts. For amyloid-forming peptides significant distortions of the N-Cα-C bond angle, indicative of sterical hidden strain, may occur in correspondence with side chain interdigitation. The correlation between the dihedral angles Δω/ψ in collagen-like models may have interesting implications for triple helix stability. PMID:24455689
A new order parameter in complex dipolar structures
NASA Astrophysics Data System (ADS)
Prosandeev, Sergey; Bellaiche, Laurent
2008-03-01
Microscopic models have been used to reveal the existence of a new order parameter that is associated with many complex dipolar structures in magnets and ferroelectrics. This overlooked order parameter involves a double cross product of the local dipoles with their positions. It provides a measure of subtle microscopic features, such as the helicity of the two domains inherent to onion states, curvature of the dipolar pattern in flower states or characteristics of set of vortices with opposite chirality (e.g., distance between vortices' centers and/or magnitude of their local dipoles). This work is mostly supported by DOE grant DE-FG02-05ER46188. We also acknowledge support from ONR grant N00014-04-1-0413 and NSF grants DMR-0701558, DMR-0404335 and DMR-0080054 (C-SPIN). Some computations were made possible thanks to the MRI Grants 0421099 and 0722625 from NSF.
Pollastri, Gianluca; Martin, Alberto JM; Mooney, Catherine; Vullo, Alessandro
2007-01-01
Background Structural properties of proteins such as secondary structure and solvent accessibility contribute to three-dimensional structure prediction, not only in the ab initio case but also when homology information to known structures is available. Structural properties are also routinely used in protein analysis even when homology is available, largely because homology modelling is lower throughput than, say, secondary structure prediction. Nonetheless, predictors of secondary structure and solvent accessibility are virtually always ab initio. Results Here we develop high-throughput machine learning systems for the prediction of protein secondary structure and solvent accessibility that exploit homology to proteins of known structure, where available, in the form of simple structural frequency profiles extracted from sets of PDB templates. We compare these systems to their state-of-the-art ab initio counterparts, and with a number of baselines in which secondary structures and solvent accessibilities are extracted directly from the templates. We show that structural information from templates greatly improves secondary structure and solvent accessibility prediction quality, and that, on average, the systems significantly enrich the information contained in the templates. For sequence similarity exceeding 30%, secondary structure prediction quality is approximately 90%, close to its theoretical maximum, and 2-class solvent accessibility roughly 85%. Gains are robust with respect to template selection noise, and significant for marginal sequence similarity and for short alignments, supporting the claim that these improved predictions may prove beneficial beyond the case in which clear homology is available. Conclusion The predictive system are publicly available at the address . PMID:17570843
Constrained maximum likelihood modal parameter identification applied to structural dynamics
NASA Astrophysics Data System (ADS)
El-Kafafy, Mahmoud; Peeters, Bart; Guillaume, Patrick; De Troyer, Tim
2016-05-01
A new modal parameter estimation method to directly establish modal models of structural dynamic systems satisfying two physically motivated constraints will be presented. The constraints imposed in the identified modal model are the reciprocity of the frequency response functions (FRFs) and the estimation of normal (real) modes. The motivation behind the first constraint (i.e. reciprocity) comes from the fact that modal analysis theory shows that the FRF matrix and therefore the residue matrices are symmetric for non-gyroscopic, non-circulatory, and passive mechanical systems. In other words, such types of systems are expected to obey Maxwell-Betti's reciprocity principle. The second constraint (i.e. real mode shapes) is motivated by the fact that analytical models of structures are assumed to either be undamped or proportional damped. Therefore, normal (real) modes are needed for comparison with these analytical models. The work done in this paper is a further development of a recently introduced modal parameter identification method called ML-MM that enables us to establish modal model that satisfies such motivated constraints. The proposed constrained ML-MM method is applied to two real experimental datasets measured on fully trimmed cars. This type of data is still considered as a significant challenge in modal analysis. The results clearly demonstrate the applicability of the method to real structures with significant non-proportional damping and high modal densities.
A new class of atomic basis functions for accurate electronic structure calculations of molecules
NASA Astrophysics Data System (ADS)
Laikov, Dimitri N.
2005-11-01
A new general approach is developed for obtaining systematic sequences of atomic single-particle basis sets for use in correlated electronic structure calculations of molecules. All the constituent functions are defined as the solutions of variational problems and are of three types: a minimal Hartree-Fock set, additional functions to represent low-lying excited configurations, and general functions for describing electron correlation. The latter are determined to minimize a functional derived from the closed-shell second-order correlation energy expression. Generally-contracted Gaussian expansions are developed to approximate these general functions in the non-relativistic case and within a scalar-relativistic approximation.
Identification of structural parameters from helicopter dynamic test data
NASA Technical Reports Server (NTRS)
Giansante, N.; Flannelly, W. G.
1974-01-01
A method is presented for obtaining the mass, stiffness, and damping parameters of a linear mathematical model, having fewer degrees of freedom than the structure it represents, directly from dynamic response measurements on the actual helicopter without a priori knowledge of the physical characteristics of the fuselage. The only input information required in the formulation is the approximate natural frequency of each mode and mobility data measured proximate to these frequencies with sinusoidal force excitation applied at only one point on the vehicle. The practicality and numerical soundness of the theoretical development was demonstrated through a computer simulation of an experimental program.
Wills, John M; Mattsson, Ann E
2012-06-06
Brooks, Johansson, and Skriver, using the LMTO-ASA method and considerable insight, were able to explain many of the ground state properties of the actinides. In the many years since this work was done, electronic structure calculations of increasing sophistication have been applied to actinide elements and compounds, attempting to quantify the applicability of DFT to actinides and actinide compounds and to try to incorporate other methodologies (i.e. DMFT) into DFT calculations. Through these calculations, the limits of both available density functionals and ad hoc methodologies are starting to become clear. However, it has also become clear that approximations used to incorporate relativity are not adequate to provide rigorous tests of the underlying equations of DFT, not to mention ad hoc additions. In this talk, we describe the result of full-potential LMTO calculations for the elemental actinides, comparing results obtained with a full Dirac basis with those obtained from scalar-relativistic bases, with and without variational spin-orbit. This comparison shows that the scalar relativistic treatment of actinides does not have sufficient accuracy to provide a rigorous test of theory and that variational spin-orbit introduces uncontrolled errors in the results of electronic structure calculations on actinide elements.
Accurate structure prediction of peptide–MHC complexes for identifying highly immunogenic antigens
Park, Min-Sun; Park, Sung Yong; Miller, Keith R.; Collins, Edward J.; Lee, Ha Youn
2013-11-01
Designing an optimal HIV-1 vaccine faces the challenge of identifying antigens that induce a broad immune capacity. One factor to control the breadth of T cell responses is the surface morphology of a peptide–MHC complex. Here, we present an in silico protocol for predicting peptide–MHC structure. A robust signature of a conformational transition was identified during all-atom molecular dynamics, which results in a model with high accuracy. A large test set was used in constructing our protocol and we went another step further using a blind test with a wild-type peptide and two highly immunogenic mutants, which predicted substantial conformational changes in both mutants. The center residues at position five of the analogs were configured to be accessible to solvent, forming a prominent surface, while the residue of the wild-type peptide was to point laterally toward the side of the binding cleft. We then experimentally determined the structures of the blind test set, using high resolution of X-ray crystallography, which verified predicted conformational changes. Our observation strongly supports a positive association of the surface morphology of a peptide–MHC complex to its immunogenicity. Our study offers the prospect of enhancing immunogenicity of vaccines by identifying MHC binding immunogens.
Effects of the projectile electronic structure on Bethe-Bloch stopping parameters for Ag
NASA Astrophysics Data System (ADS)
Moussa, D.; Damache, S.; Ouichaoui, S.
2010-06-01
Energy losses of protons and alpha particles in silver have been accurately measured under the same experimental conditions over the velocity range E=(0.192-2.595) MeV/amu using the transmission method. Deduced S(E) stopping powers are compared to most accurate ones from the literature, to values generated by the SRIM-2008 computer code and to ICRU-49 compilation. They were analyzed in the framework of modified Bethe-Bloch theory for extracting Ag target mean excitation and ionization potential, I, and Barkas effect parameter, b. Values of ( 466±5) eV and 1.20±0.01 for these two parameters were inferred from the proton S(E) data while the alpha particle data yielded values of (438±4) eV and 1.38±0.01, respectively. The ( I, b) stopping parameters thus exhibit opposite variations as the projectile charge increases, similarly as we have found previously for nickel [6]. This can be ascribed only to an effect of the projectile electronic structure at low velocities. The obtained results are discussed in comparison to previous ones reported in the literature.
A FIB-nanotomography method for accurate 3D reconstruction of open nanoporous structures.
Mangipudi, K R; Radisch, V; Holzer, L; Volkert, C A
2016-04-01
We present an automated focused ion beam nanotomography method for nanoporous microstructures with open porosity, and apply it to reconstruct nanoporous gold (np-Au) structures with ligament sizes on the order of a few tens of nanometers. This method uses serial sectioning of a well-defined wedge-shaped geometry to determine the thickness of individual slices from the changes in the sample width in successive cross-sectional images. The pore space of a selected region of the np-Au is infiltrated with ion-beam-deposited Pt composite before serial sectioning. The cross-sectional images are binarized and stacked according to the individual slice thicknesses, and then processed using standard reconstruction methods. For the image conditions and sample geometry used here, we are able to determine the thickness of individual slices with an accuracy much smaller than a pixel. The accuracy of the new method based on actual slice thickness is assessed by comparing it with (i) a reconstruction using the same cross-sectional images but assuming a constant slice thickness, and (ii) a reconstruction using traditional FIB-tomography method employing constant slice thickness. The morphology and topology of the structures are characterized using ligament and pore size distributions, interface shape distribution functions, interface normal distributions, and genus. The results suggest that the morphology and topology of the final reconstructions are significantly influenced when a constant slice thickness is assumed. The study reveals grain-to-grain variations in the morphology and topology of np-Au. PMID:26906523
Structure, vibrations, and hydrogen bond parameters of dibenzotetraaza[14]annulene
NASA Astrophysics Data System (ADS)
Gawinkowski, S.; Eilmes, J.; Waluk, J.
2010-07-01
Geometry and vibrational structure of dibenzo[ b, i][1,4,8,11]tetraaza[14]annulene (TAA) have been studied using infrared and Raman spectroscopy combined with quantum-chemical calculations. The assignments were proposed for 106 out of the total of 108 TAA vibrations, based on comparison of the theoretical predictions with the experimental data obtained for the parent molecule and its isotopomer in which the NH protons were replaced by deuterons. Reassignments were suggesteded for the NH stretching and out-of-plane vibrations. The values of the parameters of the intramolecular NH⋯N hydrogen bonds were analysed in comparison with the corresponding data for porphyrin and porphycene, molecules with the same structural motif, a cavity composed of four nitrogen atoms and two inner protons. Both experiment and calculations suggest that the molecule of TAA is not planar and is present in a trans tautomeric form, with the protons located on the opposite nitrogen atoms.
NASA Astrophysics Data System (ADS)
Chiu, Y.; Nishikawa, T.
2013-12-01
With the increasing complexity of parameter-structure identification (PSI) in groundwater modeling, there is a need for robust, fast, and accurate optimizers in the groundwater-hydrology field. For this work, PSI is defined as identifying parameter dimension, structure, and value. In this study, Voronoi tessellation and differential evolution (DE) are used to solve the optimal PSI problem. Voronoi tessellation is used for automatic parameterization, whereby stepwise regression and the error covariance matrix are used to determine the optimal parameter dimension. DE is a novel global optimizer that can be used to solve nonlinear, nondifferentiable, and multimodal optimization problems. It can be viewed as an improved version of genetic algorithms and employs a simple cycle of mutation, crossover, and selection operations. DE is used to estimate the optimal parameter structure and its associated values. A synthetic numerical experiment of continuous hydraulic conductivity distribution was conducted to demonstrate the proposed methodology. The results indicate that DE can identify the global optimum effectively and efficiently. A sensitivity analysis of the control parameters (i.e., the population size, mutation scaling factor, crossover rate, and mutation schemes) was performed to examine their influence on the objective function. The proposed DE was then applied to solve a complex parameter-estimation problem for a small desert groundwater basin in Southern California. Hydraulic conductivity, specific yield, specific storage, fault conductance, and recharge components were estimated simultaneously. Comparison of DE and a traditional gradient-based approach (PEST) shows DE to be more robust and efficient. The results of this work not only provide an alternative for PSI in groundwater models, but also extend DE applications towards solving complex, regional-scale water management optimization problems.
Doherty, Kimberly R. Talbert, Dominique R.; Trusk, Patricia B.; Moran, Diarmuid M.; Shell, Scott A.; Bacus, Sarah
2015-05-15
Safety pharmacology studies that evaluate new drug entities for potential cardiac liability remain a critical component of drug development. Current studies have shown that in vitro tests utilizing human induced pluripotent stem cell-derived cardiomyocytes (hiPS-CM) may be beneficial for preclinical risk evaluation. We recently demonstrated that an in vitro multi-parameter test panel assessing overall cardiac health and function could accurately reflect the associated clinical cardiotoxicity of 4 FDA-approved targeted oncology agents using hiPS-CM. The present studies expand upon this initial observation to assess whether this in vitro screen could detect cardiotoxicity across multiple drug classes with known clinical cardiac risks. Thus, 24 drugs were examined for their effect on both structural (viability, reactive oxygen species generation, lipid formation, troponin secretion) and functional (beating activity) endpoints in hiPS-CM. Using this screen, the cardiac-safe drugs showed no effects on any of the tests in our panel. However, 16 of 18 compounds with known clinical cardiac risk showed drug-induced changes in hiPS-CM by at least one method. Moreover, when taking into account the Cmax values, these 16 compounds could be further classified depending on whether the effects were structural, functional, or both. Overall, the most sensitive test assessed cardiac beating using the xCELLigence platform (88.9%) while the structural endpoints provided additional insight into the mechanism of cardiotoxicity for several drugs. These studies show that a multi-parameter approach examining both cardiac cell health and function in hiPS-CM provides a comprehensive and robust assessment that can aid in the determination of potential cardiac liability. - Highlights: • 24 drugs were tested for cardiac liability using an in vitro multi-parameter screen. • Changes in beating activity were the most sensitive in predicting cardiac risk. • Structural effects add in
Topology Synthesis of Structures Using Parameter Relaxation and Geometric Refinement
NASA Technical Reports Server (NTRS)
Hull, P. V.; Tinker, M. L.
2007-01-01
Typically, structural topology optimization problems undergo relaxation of certain design parameters to allow the existence of intermediate variable optimum topologies. Relaxation permits the use of a variety of gradient-based search techniques and has been shown to guarantee the existence of optimal solutions and eliminate mesh dependencies. This Technical Publication (TP) will demonstrate the application of relaxation to a control point discretization of the design workspace for the structural topology optimization process. The control point parameterization with subdivision has been offered as an alternative to the traditional method of discretized finite element design domain. The principle of relaxation demonstrates the increased utility of the control point parameterization. One of the significant results of the relaxation process offered in this TP is that direct manufacturability of the optimized design will be maintained without the need for designer intervention or translation. In addition, it will be shown that relaxation of certain parameters may extend the range of problems that can be addressed; e.g., in permitting limited out-of-plane motion to be included in a path generation problem.
THE LUMINOSITY PROFILE AND STRUCTURAL PARAMETERS OF THE ANDROMEDA GALAXY
Courteau, Stephane; Widrow, Lawrence M.; McDonald, Michael; Guhathakurta, Puragra; Zhu Yucong
2011-09-20
We have constructed an extended composite luminosity profile for the Andromeda galaxy, M31, and have decomposed it into three basic luminous structural components: a bulge, a disk, and a halo. The dust-free Spitzer/Infrared Array Camera (IRAC) imaging and extended spatial coverage of ground-based optical imaging and deep star counts allow us to map M31's structure from its center to 22 kpc along the major axis. We apply, and address the limitations of, different decomposition methods for the one-dimensional luminosity profiles and two-dimensional images. These methods include nonlinear least-squares and Bayesian Monte Carlo Markov chain analyses. The basic photometric model for M31 has a Sersic bulge with shape index n {approx_equal} 2.2 {+-} .3 and effective radius R{sub e} = 1.0 {+-} 0.2 kpc, and a dust-free exponential disk of scale length R{sub d} = 5.3 {+-} .5 kpc; the parameter errors reflect the range between different decomposition methods. Despite model covariances, the convergence of solutions based on different methods and current data suggests a stable set of structural parameters. The ellipticities ({epsilon} = 1 - b/a) of the bulge and the disk from the IRAC image are 0.37 {+-} 0.03 and 0.73 {+-} 0.03, respectively. The bulge parameter n is rather insensitive to bandpass effects and its value (2.2) suggests a first rapid formation via mergers followed by secular growth from the disk. The M31 halo has a two-dimensional power-law index {approx_equal} - 2.5 {+-} 0.2 (or -3.5 in three-dimensional), comparable to that of the Milky Way. We find that the M31 bulge light is mostly dominant over the range R{sub min} {approx}< 1.2 kpc. The disk takes over in the range 1.2 kpc {approx}< R{sub min} {approx}< 9 kpc, whereas the halo dominates at R{sub min} {approx}> 9 kpc. The stellar nucleus, bulge, disk, and halo components each contribute roughly 0.05%, 23%, 73%, and 4% of the total light of M31 out to 200 kpc along the minor axis. Nominal errors for the
Dunn, Nicholas J. H.; Noid, W. G.
2015-12-28
The present work investigates the capability of bottom-up coarse-graining (CG) methods for accurately modeling both structural and thermodynamic properties of all-atom (AA) models for molecular liquids. In particular, we consider 1, 2, and 3-site CG models for heptane, as well as 1 and 3-site CG models for toluene. For each model, we employ the multiscale coarse-graining method to determine interaction potentials that optimally approximate the configuration dependence of the many-body potential of mean force (PMF). We employ a previously developed “pressure-matching” variational principle to determine a volume-dependent contribution to the potential, U{sub V}(V), that approximates the volume-dependence of the PMF. We demonstrate that the resulting CG models describe AA density fluctuations with qualitative, but not quantitative, accuracy. Accordingly, we develop a self-consistent approach for further optimizing U{sub V}, such that the CG models accurately reproduce the equilibrium density, compressibility, and average pressure of the AA models, although the CG models still significantly underestimate the atomic pressure fluctuations. Additionally, by comparing this array of models that accurately describe the structure and thermodynamic pressure of heptane and toluene at a range of different resolutions, we investigate the impact of bottom-up coarse-graining upon thermodynamic properties. In particular, we demonstrate that U{sub V} accounts for the reduced cohesion in the CG models. Finally, we observe that bottom-up coarse-graining introduces subtle correlations between the resolution, the cohesive energy density, and the “simplicity” of the model.
Son, Sang-Kil
2011-03-01
We introduce a new numerical grid-based method on unstructured grids in the three-dimensional real-space to investigate the electronic structure of polyatomic molecules. The Voronoi-cell finite difference (VFD) method realizes a discrete Laplacian operator based on Voronoi cells and their natural neighbors, featuring high adaptivity and simplicity. To resolve multicenter Coulomb singularity in all-electron calculations of polyatomic molecules, this method utilizes highly adaptive molecular grids which consist of spherical atomic grids. It provides accurate and efficient solutions for the Schroedinger equation and the Poisson equation with the all-electron Coulomb potentials regardless of the coordinate system and the molecular symmetry. For numerical examples, we assess accuracy of the VFD method for electronic structures of one-electron polyatomic systems, and apply the method to the density-functional theory for many-electron polyatomic molecules.
Variation in the modal parameters of space structures
NASA Technical Reports Server (NTRS)
Crawley, Edward F.; Barlow, Mark S.; Van Schoor, Marthinus C.; Bicos, Andrew S.
1992-01-01
An analytic and experimental study of gravity and suspension influences on space structural test articles is presented. A modular test article including deployable, erectable, and rotary modules was assembled in three one- and two-dimensional structures. The two deployable modules utilized cable diagonal bracing rather than rigid cross members; within a bay of one of the deployable modules, the cable preload was adjustable. A friction lock was used on the alpha joint to either allow or prohibit rotary motion. Suspension systems with plunge fundamentals of 1, 2, and 5 Hz were used for ground testing to evaluate the influences of suspension stiffness. Assembly and reassembly testing was performed, as was testing on two separate shipsets at two test sites. Trends and statistical variances in modal parameters are presented as a function of force amplitude, joint preload, reassembly, shipset and suspension. Linear finite element modeling of each structure provided analytical results for 0-g unsuspended and 1-g suspended models, which are correlated with the analytical model.
Puzzarini, Cristina; Biczysko, Malgorzata
2015-05-28
State-of-the-art quantum-chemical computations have been employed to accurately determine the equilibrium structure and interaction energy of the 2-thiouracil-water complex, thus extending available reference data for biomolecule solvation patterns. The coupled-cluster level of theory in conjunction with a triple-ζ basis set has been considered together with extrapolation to the basis set limit, performed by employing second-order Møller-Plesset perturbation theory, and inclusion of core-correlation and diffuse-function corrections. On the basis of the comparison of experiment and theory for 2-thiouracil [ Puzzarini et al. Phys. Chem. Chem. Phys. 2013 , 15 , 16965 - 16975 ], structural changes due to water complexation have been pointed out. Molecular and spectroscopic properties of the 2-thiouracil-water complex have then been studied by means of the composite computational approach introduced for the molecular structure evaluation. Among the results achieved, we mention the accurate determination of the molecular dipole moment and of the spectroscopic parameters required for predicting the rotational spectrum. PMID:25474644
STRUCTURAL PARAMETERS FOR GLOBULAR CLUSTERS IN M31
Wang Song; Ma Jun
2013-08-01
In this paper, we present surface brightness profiles for 79 globular clusters in M31, using images observed with the Hubble Space Telescope, some of which are from new observations. The structural and dynamical parameters are derived from fitting the profiles to several different models for the first time. The results show that in the majority of cases, King models fit the M31 clusters just as well as Wilson models and better than Sersic models. However, there are 11 clusters best fitted by Sersic models with the Sersic index n > 2, meaning that they have cuspy central density profiles. These clusters may be the well-known core-collapsed candidates. There is a bimodality in the size distribution of M31 clusters at large radii, which is different from their Galactic counterparts. In general, the properties of clusters in M31 and the Milky Way fall in the same regions of parameter spaces. The tight correlations of cluster properties indicate a ''fundamental plane'' for clusters, which reflects some universal physical conditions and processes operating at the epoch of cluster formation.
Effect of mechanical parameters on dielectric elastomer minimum energy structures
NASA Astrophysics Data System (ADS)
Shintake, Jun; Rosset, Samuel; Floreano, Dario; Shea, Herbert R.
2013-04-01
Soft robotics may provide many advantages compared to traditional robotics approaches based on rigid materials, such as intrinsically safe physical human-robot interaction, efficient/stable locomotion, adaptive morphology, etc. The objective of this study is to develop a compliant structural actuator for soft a soft robot using dielectric elastomer minimum energy structures (DEMES). DEMES consist of a pre-stretched dielectric elastomer actuator (DEA) bonded to an initially planar flexible frame, which deforms into an out-of-plane shape which allows for large actuation stroke. Our initial goal is a one-dimensional bending actuator with 90 degree stroke. Along with frame shape, the actuation performance of DEMES depends on mechanical parameters such as thickness of the materials and pre-stretch of the elastomer membrane. We report here the characterization results on the effect of mechanical parameters on the actuator performance. The tested devices use a cm-size flexible-PCB (polyimide, 50 μm thickness) as the frame-material. For the DEA, PDMS (approximately 50 μm thickness) and carbon black mixed with silicone were used as membrane and electrode, respectively. The actuators were characterized by measuring the tip angle and the blocking force as functions of applied voltage. Different pre-stretch methods (uniaxial, biaxial and their ratio), and frame geometries (rectangular with different width, triangular and circular) were used. In order to compare actuators with different geometries, the same electrode area was used in all the devices. The results showed that the initial tip angle scales inversely with the frame width, the actuation stroke and the blocking force are inversely related (leading to an interesting design trade-off), using anisotropic pre-stretch increased the actuation stroke and the initial bending angle, and the circular frame shape exhibited the highest actuation performance.
Structural Design Parameters for Highly Birefringent Coordination Polymers.
Thompson, John R; Katz, Michael J; Williams, Vance E; Leznoff, Daniel B
2015-07-01
A series of coordination polymer materials incorporating the highly anisotropic 2-(2-pyridyl)-1,10-phenanthroline (phenpy) building block have been synthesized and structurally characterized. M(phenpy)[Au(CN)2]2 (M = Cd, Mn) are isostructural and form a 1-D chain through bridging [Au(CN)2](-) units and extend into a 2-D sheet through aurophilic interactions. M(phenpy)(H2O)[Au(CN)2]2·2H2O (M = Cd, Mn, and Zn) are also isostructural but differ from the first set via the inclusion of a water molecule into the coordination sphere, resulting in a 1-D topology through aurophilic interactions. In(phenpy)(Cl)2[Au(CN)2]·0.5H2O forms a dimer through bridging chlorides and contains a free [Au(CN)2](-) unit. In the plane of the primary crystal growth direction, the birefringence values (Δn) of 0.37(2) (Cd(phenpy)[Au(CN)2]2), 0.50(3) (In(phenpy)(Cl)2[Au(CN)2]·0.5H2O), 0.56(3) and 0.59(6) (M(phenpy)(H2O)[Au(CN)2]2·2H2O M = Cd and Zn, respectively) were determined. β, a structural parameter defined by phenpy units rotated in the A-C plane relative to the light propagation (C) direction, was found to correlate to Δn magnitudes. The addition of a carbon-carbon double bond to terpy has increased the molecular polarizability anisotropy of the building block, and all structures have reduced deviation from planarity in comparison to terpy and terpy derivative structures, leading to these higher Δn values, which are among the highest reported for crystalline solids. PMID:26098267
NASA Astrophysics Data System (ADS)
Wu, Shunguang; Hong, Lang
2008-04-01
A framework of simultaneously estimating the motion and structure parameters of a 3D object by using high range resolution (HRR) and ground moving target indicator (GMTI) measurements with template information is given. By decoupling the motion and structure information and employing rigid-body constraints, we have developed the kinematic and measurement equations of the problem. Since the kinematic system is unobservable by using only one scan HRR and GMTI measurements, we designed an architecture to run the motion and structure filters in parallel by using multi-scan measurements. Moreover, to improve the estimation accuracy in large noise and/or false alarm environments, an interacting multi-template joint tracking (IMTJT) algorithm is proposed. Simulation results have shown that the averaged root mean square errors for both motion and structure state vectors have been significantly reduced by using the template information.
Spatial dependence clusters in the estimation of forest structural parameters
NASA Astrophysics Data System (ADS)
Wulder, Michael Albert
1999-12-01
In this thesis we provide a summary of the methods by which remote sensing may be applied in forestry, while also acknowledging the various limitations which are faced. The application of spatial statistics to high spatial resolution imagery is explored as a means of increasing the information which may be extracted from digital images. A number of high spatial resolution optical remote sensing satellites that are soon to be launched will increase the availability of imagery for the monitoring of forest structure. This technological advancement is timely as current forest management practices have been altered to reflect the need for sustainable ecosystem level management. The low accuracy level at which forest structural parameters have been estimated in the past is partly due to low image spatial resolution. A large pixel is often composed of a number of surface features, resulting in a spectral value which is due to the reflectance characteristics of all surface features within that pixel. In the case of small pixels, a portion of a surface feature may be represented by a single pixel. When a single pixel represents a portion of a surface object, the potential to isolate distinct surface features exists. Spatial statistics, such as the Gets statistic, provide for an image processing method to isolate distinct surface features. In this thesis, high spatial resolution imagery sensed over a forested landscape is processed with spatial statistics to combine distinct image objects into clusters, representing individual or groups of trees. Tree clusters are a means to deal with the inevitable foliage overlap which occurs within complex mixed and deciduous forest stands. The generation of image objects, that is, clusters, is necessary to deal with the presence of spectrally mixed pixels. The ability to estimate forest inventory and biophysical parameters from image clusters generated from spatially dependent image features is tested in this thesis. The inventory
Effects of the projectile electronic structure on stopping parameters for nickel
NASA Astrophysics Data System (ADS)
Damache, S.; Ouichaoui, S.; Moussa, D.; Dib, A.
2006-08-01
The stopping powers of nickel foils for 1H+, 2H+ and 4He2+ ions have been accurately measured over the energy range E ≈ (0.166-2.725) MeV/amu. The data were compared to those reported in the literature and to values derived by the SRIM-2003 code. They were analysed in the framework of the modified Bethe-Bloch theory to extract the target mean excitation and ionisation potential, I, and the Barkas effect parameter, b. The values derived from the proton and the alpha particle data are {I = (311.11 ± 2.61) eV, b = 1.38 ± 0.01} and {I = (276.12 ± 2.84) eV, b = 1.58 ± 0.01}, respectively. They thus show opposite variations of the I and b parameters as the projectile charge increases, presumably due only to the incident ions electronic structure. Besides, the reduced stopping power data for incident deuterons show to be fully consistent with those for protons indicating no effect of the projectile isotopic structure on the target stopping parameters. These results are discussed in comparison with ones reported in the literature.
Dynamic parameters in models of atmospheric vortex structures
NASA Astrophysics Data System (ADS)
Dobryshman, E. M.; Kochina, V. G.; Letunova, T. A.
2013-09-01
Vortex simulation and the computation of fields of dynamic parameters of vortex structures (velocity, rotor velocity, and helicity) are carried out with the use of exact hydrodynamic equations in a cylindrical coordinate system. Components of centripetal and Coriolis accelerations are taken into account in the initial equations. Internal and external solutions are defined. Internal solutions ignore the disturbances of the pressure field, but they are considered in external solutions. The simulation is carried out so that the effect of accounting for spatial coordinates on the structure of the above fields is pronounced. It is shown that the initial kinetic energy of rotating motion transforms into the kinetic energy of radial and vertical velocity components in models with centripetal acceleration. In models with Coriolis acceleration, the Rossby effect is clearly pronounced. The method of an "inverse problem" is used for finding external solutions, i.e., reconstruction of the pressure field at specified velocity components. Computations have shown that tangential components mainly contribute to the velocity and helicity vortex moduli at the initial stage.
Modal parameter extraction from large operating structures using ambient excitation
James, G.H. III; Carne, T.G.; Mayes, R.L.
1995-12-31
A technique called the Natural Excitation Technique or has been developed to response extract response parameters from large operational structure when subjected to random and unmeasured forces such as wind, road noise, aerodynamics, or waves. Six applications of NExT to ambient excitation testing and NExT analysis are surveyed in this paper with a minimum of technical detail. In the first application, NExT was applied to a controlled-yaw Horizontal-Axis Wind Turbine (HAWT). By controlling the yaw degree of freedom an important class of rotating coordinate system effects are reduced. A new shape extraction procedure was applied to this data set with good results. The second application was to a free-yaw HAWT. The complexity of the response has prompted further analytical studies and the development of a specialized visualization package. The third application of NExT was to a parked three-bladed Vertical-Axis Wind Turbine (VAWT) in which traditional modal testing could not excite all modes of interest. The shape extraction process used cross-correlation functions directly in a time-domain shape-fitting routine. The fourth application was to ground transportation systems. Ongoing work to improve driver and passenger comfort in tractor-trailer vehicles and to refine automobile body and tire models will use NExT. NExT has been used to process ambient vibration data for Finite Element Model correlation and is being used to study Structural Health Monitoring with ambient excitation. Shape fitting was performed using amplitude and phase information taken directly from the cross-spectra. The final application is to an offshore structure. This work is on-going, however initial studies have found a high-modal density, high noise content, and sparse data set.
Srajer, V; Crosson, S; Schmidt, M; Key, J; Schotte, F; Anderson, S; Perman, B; Ren, Z; Teng, T Y; Bourgeois, D; Wulff, M; Moffat, K
2000-07-01
Wavelength normalization is an essential part of processing of Laue X-ray diffraction data and is critically important for deriving accurate structure-factor amplitudes. The results of wavelength normalization for Laue data obtained in nanosecond time-resolved experiments at the ID09 beamline at the European Synchrotron Radiation Facility, Grenoble, France, are presented. Several wiggler and undulator insertion devices with complex spectra were used. The results show that even in the most challenging cases, such as wiggler/undulator tandems or single-line undulators, accurate wavelength normalization does not require unusually redundant Laue data and can be accomplished using typical Laue data sets. Single-line undulator spectra derived from Laue data compare well with the measured incident X-ray spectra. Successful wavelength normalization of the undulator data was also confirmed by the observed signal in nanosecond time-resolved experiments. Single-line undulators, which are attractive for time-resolved experiments due to their high peak intensity and low polychromatic background, are compared with wigglers, based on data obtained on the same crystal. PMID:16609201
NASA Astrophysics Data System (ADS)
Maeda, Yuta; Yamaoka, Koshun; Miyamachi, Hiroki; Watanabe, Toshiki; Kunitomo, Takahiro; Ikuta, Ryoya; Yakiwara, Hiroshi; Iguchi, Masato
2015-07-01
Temporal variations of Green functions associated with the eruptive activity at Sakurajima Volcano, Japan, were estimated using an accurately controlled routinely operated signal system (ACROSS). We deconvolved 400 s waveforms of the ACROSS signal at nearby stations by a known source time function and stacked the results based on the time relative to individual eruptions and the eruption intervals; the quantities obtained by this procedure are Green functions corresponding to various stages of the eruptive activity. We found an energy decrease in the later phase of the Green functions in active eruptive periods. This energy decrease, localized in the 2-6 s window of the Green functions, is difficult to explain by contamination from volcanic earthquakes and tremors. The decrease could be more reasonably attributed to a subsurface structure change caused by the volcanic activity.
Barone, Vincenzo; Biczysko, Malgorzata Bloino, Julien; Puzzarini, Cristina
2014-07-21
Oxirane derivatives are the most used benchmarks for chiroptical spectroscopies in view of their small size and relative rigidity. The molecular structure, vibrational harmonic and anharmonic frequencies, and infrared intensities of the ground electronic states are analyzed in this paper. Equilibrium structure and harmonic force fields have been evaluated by means of high-level quantum-chemical calculations at the coupled-cluster level including single and double excitations together with a perturbative treatment of triples (CCSD(T)). Extrapolation to the complete basis-set limit as well as core-correlation effects have also been taken into account. Anharmonic contributions have been computed at the CCSD(T)/cc-pVTZ level for trans-2,3-dideuterooxirane. These data can serve as references to evaluate the accuracy of less expensive computational approaches rooted in the density functional theory (DFT). The latter have been used within hybrid CC/DFT approaches, which have been applied to simulate fully anharmonic infrared (IR) spectra. Finally, the best theoretical estimates of the equilibrium structures and vibrational wavenumbers are compared to the most accurate experimental data and show in all cases very good agreement, i.e., within 0.001 Å, 0.1 deg, 10 cm{sup −1}, and 0.5 km mol{sup −1}, for bond lengths, angles, wavenumbers, and IR intensities, respectively.
THE STAR CLUSTERS OF THE LARGE MAGELLANIC CLOUD: STRUCTURAL PARAMETERS
Werchan, Felicia; Zaritsky, Dennis E-mail: dzaritsky@as.arizona.edu
2011-08-15
We present and analyze the radial luminosity profiles of a sample of 1066 stellar clusters in the Large Magellanic Cloud (LMC). By design, this study closely follows the compilation by Hill and Zaritsky of the structural parameters of stellar clusters in the Small Magellanic Cloud (SMC). Both King and Elson-Fall-Freeman model profiles are fit to V-band surface brightness profiles measured from the Magellanic Cloud Photometric Survey images. We tabulate the concentration, central surface brightness, tidal radii, 90% enclosed luminosity radii (r{sub 90}), and local background luminosity density. Over two-thirds of the clusters in the sample are adequately fit by one or both of these models. One notable and systematic exception, as in the SMC, is those clusters that lack a central brightness concentration, the 'ring' clusters. While the bulk properties of the clusters are similar between the LMC and SMC populations, we find that the LMC lacks clusters that are as large, either in terms of core radii or r{sub 90}, as the largest in the SMC, perhaps a signature of larger tidal stresses in the LMC.
Analysing DNA structural parameters using a mesoscopic model
NASA Astrophysics Data System (ADS)
Amarante, Tauanne D.; Weber, Gerald
2014-03-01
The Peyrard-Bishop model is a mesoscopic approximation to model DNA and RNA molecules. Several variants of this model exists, from 3D Hamiltonians, including torsional angles, to simpler 2D versions. Currently, we are able to parametrize the 2D variants of the model which allows us to extract important information about the molecule. For example, with this technique we were able recently to obtain the hydrogen bonds of RNA from melting temperatures, which previously were obtainable only from NMR measurements. Here, we take the 3D torsional Hamiltonian and set the angles to zero. Curiously, in doing this we do not recover the traditional 2D Hamiltonians. Instead, we obtain a different 2D Hamiltonian which now includes a base pair step distance, commonly known as rise. A detailed knowledge of the rise distance is important as it determines the overall length of the DNA molecule. This 2D Hamiltonian provides us with the exciting prospect of obtaining DNA structural parameters from melting temperatures. Our results of the rise distance at low salt concentration are in good qualitative agreement with those from several published x-ray measurements. We also found an important dependence of the rise distance with salt concentration. In contrast to our previous calculations, the elastic constants now show little dependence with salt concentrations which appears to be closer to what is seen experimentally in DNA flexibility experiments.
Extracting galactic structure parameters from multivariated density estimation
NASA Technical Reports Server (NTRS)
Chen, B.; Creze, M.; Robin, A.; Bienayme, O.
1992-01-01
Multivariate statistical analysis, including includes cluster analysis (unsupervised classification), discriminant analysis (supervised classification) and principle component analysis (dimensionlity reduction method), and nonparameter density estimation have been successfully used to search for meaningful associations in the 5-dimensional space of observables between observed points and the sets of simulated points generated from a synthetic approach of galaxy modelling. These methodologies can be applied as the new tools to obtain information about hidden structure otherwise unrecognizable, and place important constraints on the space distribution of various stellar populations in the Milky Way. In this paper, we concentrate on illustrating how to use nonparameter density estimation to substitute for the true densities in both of the simulating sample and real sample in the five-dimensional space. In order to fit model predicted densities to reality, we derive a set of equations which include n lines (where n is the total number of observed points) and m (where m: the numbers of predefined groups) unknown parameters. A least-square estimation will allow us to determine the density law of different groups and components in the Galaxy. The output from our software, which can be used in many research fields, will also give out the systematic error between the model and the observation by a Bayes rule.
Extracting Structure Parameters of Dimers for Molecular Tunneling Ionization Model
NASA Astrophysics Data System (ADS)
Song-Feng, Zhao; Fang, Huang; Guo-Li, Wang; Xiao-Xin, Zhou
2016-03-01
We determine structure parameters of the highest occupied molecular orbital (HOMO) of 27 dimers for the molecular tunneling ionization (so called MO-ADK) model of Tong et al. [Phys. Rev. A 66 (2002) 033402]. The molecular wave functions with correct asymptotic behavior are obtained by solving the time-independent Schrödinger equation with B-spline functions and molecular potentials which are numerically created using the density functional theory. We examine the alignment-dependent tunneling ionization probabilities from MO-ADK model for several molecules by comparing with the molecular strong-field approximation (MO-SFA) calculations. We show the molecular Perelomov–Popov–Terent'ev (MO-PPT) can successfully give the laser wavelength dependence of ionization rates (or probabilities). Based on the MO-PPT model, two diatomic molecules having valence orbital with antibonding systems (i.e., Cl2, Ne2) show strong ionization suppression when compared with their corresponding closest companion atoms. Supported by National Natural Science Foundation of China under Grant Nos. 11164025, 11264036, 11465016, 11364038, the Specialized Research Fund for the Doctoral Program of Higher Education of China under Grant No. 20116203120001, and the Basic Scientific Research Foundation for Institution of Higher Learning of Gansu Province
Extracting Structure Parameters of Dimers for Molecular Tunneling Ionization Model
NASA Astrophysics Data System (ADS)
Zhao, Song-Feng; Huang, Fang; Wang, Guo-Li; Zhou, Xiao-Xin
2016-03-01
We determine structure parameters of the highest occupied molecular orbital (HOMO) of 27 dimers for the molecular tunneling ionization (so called MO-ADK) model of Tong et al. [Phys. Rev. A 66 (2002) 033402]. The molecular wave functions with correct asymptotic behavior are obtained by solving the time-independent Schrödinger equation with B-spline functions and molecular potentials which are numerically created using the density functional theory. We examine the alignment-dependent tunneling ionization probabilities from MO-ADK model for several molecules by comparing with the molecular strong-field approximation (MO-SFA) calculations. We show the molecular Perelomov-Popov-Terent'ev (MO-PPT) can successfully give the laser wavelength dependence of ionization rates (or probabilities). Based on the MO-PPT model, two diatomic molecules having valence orbital with antibonding systems (i.e., Cl2, Ne2) show strong ionization suppression when compared with their corresponding closest companion atoms. Supported by National Natural Science Foundation of China under Grant Nos. 11164025, 11264036, 11465016, 11364038, the Specialized Research Fund for the Doctoral Program of Higher Education of China under Grant No. 20116203120001, and the Basic Scientific Research Foundation for Institution of Higher Learning of Gansu Province
NASA Technical Reports Server (NTRS)
Cantrell, John H.; Yost, William T.
1990-01-01
The effects of material structure on the nonlinearity parameters are reviewed. Problems discussed include definition of nonlinearity parameters, square-law nonlinearity and collinear beam-mixing, structure dependence of the nonlinearity parameters, negative nonlinearity parameters, and implications for materials characterization.
Moran, Robert F; McKay, David; Pickard, Chris J; Berry, Andrew J; Griffin, John M; Ashbrook, Sharon E
2016-04-21
The structural chemistry of materials containing low levels of nonstoichiometric hydrogen is difficult to determine, and producing structural models is challenging where hydrogen has no fixed crystallographic site. Here we demonstrate a computational approach employing ab initio random structure searching (AIRSS) to generate a series of candidate structures for hydrous wadsleyite (β-Mg2SiO4 with 1.6 wt% H2O), a high-pressure mineral proposed as a repository for water in the Earth's transition zone. Aligning with previous experimental work, we solely consider models with Mg3 (over Mg1, Mg2 or Si) vacancies. We adapt the AIRSS method by starting with anhydrous wadsleyite, removing a single Mg(2+) and randomly placing two H(+) in a unit cell model, generating 819 candidate structures. 103 geometries were then subjected to more accurate optimisation under periodic DFT. Using this approach, we find the most favourable hydration mechanism involves protonation of two O1 sites around the Mg3 vacancy. The formation of silanol groups on O3 or O4 sites (with loss of stable O1-H hydroxyls) coincides with an increase in total enthalpy. Importantly, the approach we employ allows observables such as NMR parameters to be computed for each structure. We consider hydrous wadsleyite (∼1.6 wt%) to be dominated by protonated O1 sites, with O3/O4-H silanol groups present as defects, a model that maps well onto experimental studies at higher levels of hydration (J. M. Griffin et al., Chem. Sci., 2013, 4, 1523). The AIRSS approach adopted herein provides the crucial link between atomic-scale structure and experimental studies. PMID:27020937
NASA Astrophysics Data System (ADS)
Hsieh, H. P.; Sung, K. B.; Hsu, F. W.
2014-05-01
Diffuse reflectance spectroscopy has been applied as a non-invasive method to measure tissue optical properties, which are associate with anatomical information. The algorithm widely used to extract, optical parameters from reflectance spectra is the regression method, which is time-consuming and frequently converge to local maxima. In this study, the effects of parameters changes on spectra are analyzed in different fiber geometries, source-detector separations and wavelengths. In the end of this paper, a new fitting algorithm is proposed base on parameters features found. The new algorithm is expected to enhance the accuracy of parameters extracted and save 75% of the process time.
Projection filters for modal parameter estimate for flexible structures
NASA Technical Reports Server (NTRS)
Huang, Jen-Kuang; Chen, Chung-Wen
1987-01-01
Single-mode projection filters are developed for eigensystem parameter estimates from both analytical results and test data. Explicit formulations of these projection filters are derived using the pseudoinverse matrices of the controllability and observability matrices in general use. A global minimum optimization algorithm is developed to update the filter parameters by using interval analysis method. Modal parameters can be attracted and updated in the global sense within a specific region by passing the experimental data through the projection filters. For illustration of this method, a numerical example is shown by using a one-dimensional global optimization algorithm to estimate model frequencies and dampings.
Parameter Optimization and Field Validation of the Functional–Structural Model GREENLAB for Maize
GUO, YAN; MA, YUNTAO; ZHAN, ZHIGANG; LI, BAOGUO; DINGKUHN, MICHAEL; LUQUET, DELPHINE; DE REFFYE, PHILIPPE
2006-01-01
• Background and Aims There are three reasons for the increasing demand for crop models that build the plant on the basis of architectural principles and organogenetic processes: (1) realistic concepts for developing new crops need to be guided by such models; (2) there is an increasing interest in crop phenotypic plasticity, based on variable architecture and morphology; and (3) engineering of mechanized cropping systems requires information on crop architecture. The functional–structural model GREENLAB was recently presented that simulates resource-dependent plasticity of plant architecture. This study introduces a new methodology for crop parameter optimization against measured data called multi-fitting, validates the calibrated model for maize with independent field data, and describes a technique for 3D visualization of outputs. • Methods Maize was grown near Beijing during the 2000, 2001 and 2003 (two sowing dates) summer seasons in a block design with four to five replications. Detailed morphological and topological observations were made on the plant architecture throughout the development of the four crops. Data obtained in 2000 was used to establish target files for parameter optimization using the generalized least square method, and parameter accuracy was evaluated by coefficient of variance. In situ plant digitization was used to establish 3D symbol files for organs that were then used to translate model outputs directly into 3D representations for each time step of model execution. •Key Results and Conclusions Multi-fitting against several target files obtained at different growth stages gave better parameter accuracy than single fitting at maturity only, and permitted extracting generic organ expansion kinetics from the static observations. The 2000 model gave excellent predictions of plant architecture and vegetative growth for the other three seasons having different temperature regimes, but predictions of inter-seasonal variability of
Microwave spectra and gas phase structural parameters for N-hydroxypyridine-2(1H)-thione.
Daly, Adam M; Mitchell, Erik G; Sanchez, Daniel A; Block, Eric; Kukolich, Stephen G
2011-12-22
The microwave spectrum for N-hydroxypyridine-2(1H)-thione (pyrithione) was measured in the frequency range 6-18 GHz, providing accurate rotational constants and nitrogen quadrupole coupling strengths for three isotopologues, C(5)H(4)(32)S(14)NOH, C(5)H(4)(32)S(14)NOD, and C(5)H(4)(34)S(14)NOH. Pyrithione was found to be in a higher concentration in the gas phase than the other tautomer, 2-mercaptopyridine-N-oxide (MPO). Microwave spectroscopy is best suited to determine which structure predominates in the gas phase. The measured rotational constants were used to accurately determine the coordinates of the substituted atoms and provided sufficient data to determine some of the important structural parameters for pyrithione, the only tautomer observed in the present work. The spectra were obtained using a pulsed-beam Fourier transform microwave spectrometer, with sufficient resolution to allow accurate measurements of the (14)N nuclear quadrupole hyperfine interactions. Ab initio calculations provided structural parameters and quadrupole coupling strengths that are in very good agreement with measured values. The experimental rotational constants for the parent compound are A = 3212.10(1), B = 1609.328(7), and C = 1072.208(6) MHz, yielding the inertial defect Δ(0) = -0.023 amu·Å(2) for the C(5)H(4)(32)S(14)NOH isotopologue. The observed near zero inertial defect clearly indicates a planar structure. The least-squares fit structural analysis yielded the experimental bond lengths R(O-H) = 0.93(2) Å, R(C-S) = 1.66(2) Å, and angle (N-O-H) = 105(4)° for the ground state structure. PMID:22070758
Constructing Approximate Confidence Intervals for Parameters with Structural Equation Models
ERIC Educational Resources Information Center
Cheung, Mike W. -L.
2009-01-01
Confidence intervals (CIs) for parameters are usually constructed based on the estimated standard errors. These are known as Wald CIs. This article argues that likelihood-based CIs (CIs based on likelihood ratio statistics) are often preferred to Wald CIs. It shows how the likelihood-based CIs and the Wald CIs for many statistics and psychometric…
Identification of random variation in structures and their parameter estimates.
Farrar, C. R.; Aumann, R. J.; McCarty, A. A.; Olson, C. C.
2002-01-01
Structures that are members of an ensemble of nominally identical systems actually differ due to variations in details among individuals. Furthermore, there are variations in the system response of an individual structure that can be attributed to unmeasured conditions (such as temperature and humidity) that are present during experiments. Finally, noise is present in all measurements of structural excitations and responses. For these reasons, there is always random variation associated with the characterizations of structural dynamic systems, and descriptions of results must be in statistical or probabilistic terms;. This study identifies and assesses the sources and the degrees of randomness in a metric of structural dynamics of a given system through experiments and analysis.
Derivation of tree stem structural parameters from static terrestrial laser scanning data
NASA Astrophysics Data System (ADS)
Tian, Wei; Lin, Yi; Liu, Yajing; Niu, Zheng
2014-11-01
Accurate tree-level characteristic information is increasingly demanded for forest management and environment protection. The cutting-edge remote sensing technique of terrestrial laser scanning (TLS) shows the potential of filling this gap. This study focuses on exploring the methods for deriving various tree stem structural parameters, such as stem position, diameter at breast height (DBH), the degree of stem shrinkage, and the elevation angle and azimuth angle of stem inclination. The data for test was collected with a Leica HDS6100 TLS system in Seurasaari, Southern Finland in September 2010. In the field, the reference positions and DBHs of 100 trees were measured manually. The isolation of individual trees is based on interactive segmentation of point clouds. The estimation of stem position and DBH is based on the schematic of layering and then least-square-based circle fitting in each layer. The slope of robust fit line between the height of each layer and DBH is used to characterize the stem shrinkage. The elevation angle of stem inclination is described by the angle between the ground plane and the fitted stem axis. The angle between the north direction and the fitted stem axis gives the azimuth angle of stem inclination. The estimation of the DBHs performed with R square (R2) of 0.93 and root mean square error (RMSE) of 0.038m.The average angle corresponding to stem shrinkage is -1.86°. The elevation angles of stem inclinations are ranged from 31° to 88.3°. The results have basically validated TLS for deriving multiple structural parameters of stem, which help better grasp tree specialties.
Direct structural parameter identification by modal test results
NASA Technical Reports Server (NTRS)
Chen, J.-C.; Kuo, C.-P.; Garba, J. A.
1983-01-01
A direct identification procedure is proposed to obtain the mass and stiffness matrices based on the test measured eigenvalues and eigenvectors. The method is based on the theory of matrix perturbation in which the correct mass and stiffness matrices are expanded in terms of analytical values plus a modification matrix. The simplicity of the procedure enables real time operation during the structural testing.
Influence of solvent parameters on structure of polyhydroxybutyrate films
NASA Astrophysics Data System (ADS)
Olkhov, A. A.; Iordanskii, A. L.
2012-07-01
The polarity of dissolvent increases the perfection of a crystal structure and decrease the amount of amorphous phase of polyhydroxybutyrate (PHB). It is shown, that the amount of a defective crystal phase in films PHB is directly proportional to magnitude of cohesive energy of dissolvent.
Cosmological parameter estimation with large scale structure observations
NASA Astrophysics Data System (ADS)
Di Dio, Enea; Montanari, Francesco; Durrer, Ruth; Lesgourgues, Julien
2014-01-01
We estimate the sensitivity of future galaxy surveys to cosmological parameters, using the redshift dependent angular power spectra of galaxy number counts, Cl(z1,z2), calculated with all relativistic corrections at first order in perturbation theory. We pay special attention to the redshift dependence of the non-linearity scale and present Fisher matrix forecasts for Euclid-like and DES-like galaxy surveys. We compare the standard P(k) analysis with the new Cl(z1,z2) method. We show that for surveys with photometric redshifts the new analysis performs significantly better than the P(k) analysis. For spectroscopic redshifts, however, the large number of redshift bins which would be needed to fully profit from the redshift information, is severely limited by shot noise. We also identify surveys which can measure the lensing contribution and we study the monopole, C0(z1,z2).
Structure parameters in rotating Couette-Poiseuille channel flow
NASA Technical Reports Server (NTRS)
Knightly, George H.; Sather, D.
1986-01-01
It is well-known that a number of steady state problems in fluid mechanics involving systems of nonlinear partial differential equations can be reduced to the problem of solving a single operator equation of the form: v + lambda Av + lambda B(v) = 0, v is the summation of H, lambda is the summation of one-dimensional Euclid space, where H is an appropriate (real or complex) Hilbert space. Here lambda is a typical load parameter, e.g., the Reynolds number, A is a linear operator, and B is a quadratic operator generated by a bilinear form. In this setting many bifurcation and stability results for problems were obtained. A rotating Couette-Poiseuille channel flow was studied, and it showed that, in general, the superposition of a Poiseuille flow on a rotating Couette channel flow is destabilizing.
NASA Astrophysics Data System (ADS)
Bruntt, H.
2009-10-01
Context: The CoRoT satellite has provided high-quality light curves of several solar-like stars. Analysis of these light curves provides oscillation frequencies that make it possible to probe the interior of the stars. However, additional constraints on the fundamental parameters of the stars are important for the theoretical modelling to be successful. Aims: We estimate the fundamental parameters (mass, radius, and luminosity) of the first four solar-like targets to be observed in the asteroseismic field. In addition, we determine their effective temperature, metallicity, and detailed abundance patterns. Methods: To constrain the stellar mass, radius and age we used the shotgun software, which compares the location of the stars in the Hertzsprung-Russell diagram with theoretical evolution models. This method takes the uncertainties of the observed parameters into account, including the large separation determined from the solar-like oscillations. We determined the effective temperatures and abundance patterns in the stars from the analysis of high-resolution spectra obtained with the HARPS, NARVAL, ELODIE and FEROS spectrographs. Results: We determined the mass, radius, and luminosity of the four CoRoT targets to within 5{-}10%, 2{-}4% and 5{-}13%, respectively. The quality of the stellar spectra determines how well we can constrain the effective temperature. For the two best spectra we get 1-σ uncertainties below 60 K and 100{-}150 K for the other two. The uncertainty on the surface gravity is less than 0.08 dex for three stars, while it is 0.15 dex for HD 181906. The reason for the larger uncertainty is that the spectrum has two components with a luminosity ratio of L_p/Ls = 0.50±0.15. While Hipparcos astrometric data strongly suggest it is a binary star, we find evidence that the fainter star may be a background star, since it is less luminous but hotter.
Quantitative structure parameters from the NMR spectroscopy of quadrupolar nuclei
Perras, Frédéric A.
2016-01-01
Nuclear magnetic resonance (NMR) spectroscopy is one of the most important characterization tools in chemistry, however, 3/4 of the NMR active nuclei are underutilized due to their quadrupolar nature. This short review centers on the development of methods that use solid-state NMR of quadrupolar nuclei for obtaining quantitative structural information. Namely, techniques using dipolar recoupling as well as the resolution afforded by double-rotation are presented for the measurement of spin–spin coupling between quadrupoles, enabling the measurement of internuclear distances and connectivities. Two-dimensional
Quantitative structure parameters from the NMR spectroscopy of quadrupolar nuclei
Perras, Frederic A.
2015-12-15
Here, nuclear magnetic resonance (NMR) spectroscopy is one of the most important characterization tools in chemistry, however, 3/4 of the NMR active nuclei are underutilized due to their quadrupolar nature. This short review centers on the development of methods that use solid-state NMR of quadrupolar nuclei for obtaining quantitative structural information. Namely, techniques using dipolar recoupling as well as the resolution afforded by double-rotation are presented for the measurement of spin–spin coupling between quadrupoles, enabling the measurement of internuclear distances and connectivities.
On the atmospheric structure and fundamental parameters of red supergiants
NASA Astrophysics Data System (ADS)
Wittkowski, M.; Arroyo-Torres, B.; Marcaide, J. M.; Abellan, F. J.; Chiavassa, A.; Freytag, B.; Scholz, M.; Wood, P. R.; Hauschildt, P. H.
2015-01-01
We present near-infrared spectro-interferometric studies of red supergiant (RSG) stars using the VLTI/AMBER instrument, which are compared to previously obtained similar observations of AGB stars. Our observations indicate spatially extended atmospheric molecular layers of water vapor and CO, similar as previously observed for Mira stars. Data of VY~CMa indicate that the molecular layers are asymmetric, possibly clumpy. Thanks to the spectro-interferometric capabilities of the VLTI/AMBER instrument, we can isolate continuum bandpasses, estimate fundamental parameters of our sources, locate them in the HR diagram, and compare their positions to recent evolutionary tracks. For the example of VY CMa, this puts it close to evolutionary tracks of initial mass 25-32 M ⊙. Comparisons of our data to hydrostatic model atmospheres, 3d simulations of convection, and 1d dynamic model atmospheres based on self-excited pulsation models indicate that none of these models can presently explain the observed atmospheric extensions for RSGs. The mechanism that levitates the atmospheres of red supergiant is thus a currently unsolved problem.
Cosmological parameter estimation with large scale structure observations
Dio, Enea Di; Montanari, Francesco; Durrer, Ruth; Lesgourgues, Julien E-mail: Francesco.Montanari@unige.ch E-mail: Julien.Lesgourgues@cern.ch
2014-01-01
We estimate the sensitivity of future galaxy surveys to cosmological parameters, using the redshift dependent angular power spectra of galaxy number counts, C{sub ℓ}(z{sub 1},z{sub 2}), calculated with all relativistic corrections at first order in perturbation theory. We pay special attention to the redshift dependence of the non-linearity scale and present Fisher matrix forecasts for Euclid-like and DES-like galaxy surveys. We compare the standard P(k) analysis with the new C{sub ℓ}(z{sub 1},z{sub 2}) method. We show that for surveys with photometric redshifts the new analysis performs significantly better than the P(k) analysis. For spectroscopic redshifts, however, the large number of redshift bins which would be needed to fully profit from the redshift information, is severely limited by shot noise. We also identify surveys which can measure the lensing contribution and we study the monopole, C{sub 0}(z{sub 1},z{sub 2})
Estimates of forest structure parameters from GLAS data and multi-angle imaging spectrometer data
NASA Astrophysics Data System (ADS)
Yu, Ying; Yang, Xiguang; Fan, Wenyi
2015-06-01
Quantitative estimates of forest vertical and spatial distribution using remote sensing technology play an important role in better understanding forest ecosystem function, forest carbon storage and the global carbon cycle. Although most remote sensing systems can provide horizontal distribution of canopies, information concerning the vertical distribution of canopies cannot be detected. Fortunately, laser radars have become available, such as GLAS (Geoscience laser altimeter system). Because laser radar can penetrate foliage, it is superior to other remote sensing technologies for detecting vertical forest structure and has higher accuracy. GLAS waveform data were used in this study to retrieve average tree height and biomass in a GLAS footprint area in Heilongjiang Province. However, GLAS data are not spatially continuous. To fill the gaps, MISR (multi- angle imaging spectrometer) spectral radiance was chosen to predict the regional continuous tree height by developing a multivariate linear regression model. We compared tree height estimated by the regression model and GLAS data. The results confirmed that estimates of tree height and biomass based on GLAS data are considerably more accurate than estimates based on traditional methods. The accuracy is approximately 90%. MISR can be used to estimate tree height in continuous areas with a robust regression model. The R2, precision and root mean square error of the regression model were 0.8, 83% and 1 m, respectively. This study provides an important reference for mapping forest vertical parameters.
Further Studies of Forest Structure Parameter Retrievals Using the Echidna® Ground-Based Lidar
NASA Astrophysics Data System (ADS)
Strahler, A. H.; Yao, T.; Zhao, F.; Yang, X.; Schaaf, C.; Wang, Z.; Li, Z.; Woodcock, C. E.; Culvenor, D.; Jupp, D.; Newnham, G.; Lovell, J.
2012-12-01
Ongoing work with the Echidna® Validation Instrument (EVI), a full-waveform, ground-based scanning lidar (1064 nm) developed by Australia's CSIRO and deployed by Boston University in California conifers (2008) and New England hardwood and softwood (conifer) stands (2007, 2009, 2010), confirms the importance of slope correction in forest structural parameter retrieval; detects growth and disturbance over periods of 2-3 years; provides a new way to measure the between-crown clumping factor in leaf area index retrieval using lidar range; and retrieves foliage profiles with more lower-canopy detail than a large-footprint aircraft scanner (LVIS), while simulating LVIS foliage profiles accurately from a nadir viewpoint using a 3-D point cloud. Slope correction is important for accurate retrieval of forest canopy structural parameters, such as mean diameter at breast height (DBH), stem count density, basal area, and above-ground biomass. Topographic slope can induce errors in parameter retrievals because the horizontal plane of the instrument scan, which is used to identify, measure, and count tree trunks, will intersect trunks below breast height in the uphill direction and above breast height in the downhill direction. A test of three methods at southern Sierra Nevada conifer sites improved the range of correlations of these EVI-retrieved parameters with field measurements from 0.53-0.68 to 0.85-0.93 for the best method. EVI scans can detect change, including both growth and disturbance, in periods of two to three years. We revisited three New England forest sites scanned in 2007-2009 or 2007-2010. A shelterwood stand at the Howland Experimental Forest, Howland, Maine, showed increased mean DBH, above-ground biomass and leaf area index between 2007 and 2009. Two stands at the Harvard Forest, Petersham, Massachusetts, suffered reduced leaf area index and reduced stem count density as the result of an ice storm that damaged the stands. At one stand, broken tops were
Shi, Deheng; Liu, Qionglan; Sun, Jinfeng; Zhu, Zunlue
2014-03-25
The potential energy curves (PECs) of 28 Ω states generated from the 12 states (X(4)Σ(-), 1(2)Π, 1(2)Σ(-), 1(2)Δ, 1(2)Σ(+), 2(2)Π, A(4)Π, B(4)Σ(-), 3(2)Π, 1(6)Σ(-), 2(2)Σ(-) and 1(6)Π) of the BN(+) cation are studied for the first time for internuclear separations from about 0.1 to 1.0 nm using an ab initio quantum chemical method. All the Λ-S states correlate to the first four dissociation channels. The 1(6)Σ(-), 3(2)Π and A(4)Π states are found to be the inverted ones. The 1(2)Σ(+), 2(2)Π, 3(2)Π and 2(2)Σ(-) states are found to possess the double well. The PECs are calculated by the complete active space self-consistent field method, which is followed by the internally contracted multireference configuration interaction approach with the Davidson correction. Core-valence correlation correction is included by a cc-pCV5Z basis set. Scalar relativistic correction is calculated by the third-order Douglas-Kroll Hamiltonian approximation at the level of a cc-pV5Z basis set. The convergent behavior of present calculations is discussed with respect to the basis set and level of theory. The spin-orbit coupling is accounted for by the state interaction approach with the Breit-Pauli Hamiltonian using the all-electron cc-pCV5Z basis set. All the PECs are extrapolated to the complete basis set limit. The spectroscopic parameters are obtained, and the vibrational properties of 1(2)Σ(+), 2(2)Π, 3(2)Π and 2(2)Σ(-) states are evaluated. Analyses demonstrate that the spectroscopic parameters reported here can be expected to be reliably predicted ones. The conclusion is gained that the effect of spin-orbit coupling on the spectroscopic parameters are not obvious almost for all the Λ-S states involved in the present paper. PMID:24334021
Shi, Deheng; Li, Peiling; Sun, Jinfeng; Zhu, Zunlue
2014-01-01
The potential energy curves (PECs) of 28 Ω states generated from 9 Λ-S states (X(2)Π, 1(4)Π, 1(6)Π, 1(2)Σ(+), 1(4)Σ(+), 1(6)Σ(+), 1(4)Σ(-), 2(4)Π and 1(4)Δ) are studied for the first time using an ab initio quantum chemical method. All the 9 Λ-S states correlate to the first two dissociation limits, N((4)Su)+Se((3)Pg) and N((4)Su)+Se((3)Dg), of NSe radical. Of these Λ-S states, the 1(6)Σ(+), 1(4)Σ(+), 1(6)Π, 2(4)Π and 1(4)Δ are found to be rather weakly bound states. The 1(2)Σ(+) is found to be unstable and has double wells. And the 1(6)Σ(+), 1(4)Σ(+), 1(4)Π and 1(6)Π are found to be the inverted ones with the SO coupling included. The PEC calculations are made by the complete active space self-consistent field method, which is followed by the internally contracted multireference configuration interaction approach with the Davidson modification. The spin-orbit coupling is accounted for by the state interaction approach with the Breit-Pauli Hamiltonian. The convergence of the present calculations is discussed with respect to the basis set and the level of theory. Core-valence correlation corrections are included with a cc-pCVTZ basis set. Scalar relativistic corrections are calculated by the third-order Douglas-Kroll Hamiltonian approximation at the level of a cc-pV5Z basis set. All the PECs are extrapolated to the complete basis set limit. The variation with internuclear separation of spin-orbit coupling constants is discussed in brief for some Λ-S states with one shallow well on each PEC. The spectroscopic parameters of 9 Λ-S and 28 Ω states are determined by fitting the first ten vibrational levels whenever available, which are calculated by solving the rovibrational Schrödinger equation with Numerov's method. The splitting energy in the X(2)Π Λ-S state is determined to be about 864.92 cm(-1), which agrees favorably with the measurements of 891.80 cm(-1). Moreover, other spectroscopic parameters of Λ-S and Ω states involved here are
Accurate wavelength calibration method for flat-field grating spectrometers.
Du, Xuewei; Li, Chaoyang; Xu, Zhe; Wang, Qiuping
2011-09-01
A portable spectrometer prototype is built to study wavelength calibration for flat-field grating spectrometers. An accurate calibration method called parameter fitting is presented. Both optical and structural parameters of the spectrometer are included in the wavelength calibration model, which accurately describes the relationship between wavelength and pixel position. Along with higher calibration accuracy, the proposed calibration method can provide information about errors in the installation of the optical components, which will be helpful for spectrometer alignment. PMID:21929865
Badran, Yasser Ali; Abdelaziz, Alsayed Saad; Shehab, Mohamed Ahmed; Mohamed, Hazem Abdelsabour Dief; Emara, Absel-Aziz Ali; Elnabtity, Ali Mohamed Ali; Ghanem, Maged Mohammed; ELHelaly, Hesham Abdel Azim
2016-01-01
Objective: The objective was to determine the predicting success of shock wave lithotripsy (SWL) using a combination of computed tomography based metric parameters to improve the treatment plan. Patients and Methods: Consecutive 180 patients with symptomatic upper urinary tract calculi 20 mm or less were enrolled in our study underwent extracorporeal SWL were divided into two main groups, according to the stone size, Group A (92 patients with stone ≤10 mm) and Group B (88 patients with stone >10 mm). Both groups were evaluated, according to the skin to stone distance (SSD) and Hounsfield units (≤500, 500–1000 and >1000 HU). Results: Both groups were comparable in baseline data and stone characteristics. About 92.3% of Group A rendered stone-free, whereas 77.2% were stone-free in Group B (P = 0.001). Furthermore, in both group SWL success rates was a significantly higher for stones with lower attenuation <830 HU than with stones >830 HU (P < 0.034). SSD were statistically differences in SWL outcome (P < 0.02). Simultaneous consideration of three parameters stone size, stone attenuation value, and SSD; we found that stone-free rate (SFR) was 100% for stone attenuation value <830 HU for stone <10 mm or >10 mm but total number SWL sessions and shock waves required for the larger stone group were higher than in the smaller group (P < 0.01). Furthermore, SFR was 83.3% and 37.5% for stone <10 mm, mean HU >830, SSD 90 mm and SSD >120 mm, respectively. On the other hand, SFR was 52.6% and 28.57% for stone >10 mm, mean HU >830, SSD <90 mm and SSD >120 mm, respectively. Conclusion: Stone size, stone density (HU), and SSD is simple to calculate and can be reported by radiologists to applying combined score help to augment predictive power of SWL, reduce cost, and improving of treatment strategies. PMID:27141192
NASA Astrophysics Data System (ADS)
Suchomska, K.; Graczyk, D.; Smolec, R.; Pietrzyński, G.; Gieren, W.; Stȩpień, K.; Konorski, P.; Pilecki, B.; Villanova, S.; Thompson, I. B.; Górski, M.; Karczmarek, P.; Wielgórski, P.; Anderson, R. I.
2015-07-01
We have analyzed the double-lined eclipsing binary system ASAS J180057-2333.8 from the All Sky Automated Survey (ASAS) catalogue. We measure absolute physical and orbital parameters for this system based on archival V-band and I-band ASAS photometry, as well as on high-resolution spectroscopic data obtained with ESO 3.6 m/HARPS and CORALIE spectrographs. The physical and orbital parameters of the system were derived with an accuracy of about 0.5-3 per cent. The system is a very rare configuration of two bright well-detached giants of spectral types K1 and K4 and luminosity class II. The radii of the stars are R1 = 52.12 ± 1.38 and R2 = 67.63 ± 1.40 R⊙ and their masses are M1 = 4.914 ± 0.021 and M2 = 4.875 ± 0.021 M⊙. The exquisite accuracy of 0.5 per cent obtained for the masses of the components is one of the best mass determinations for giants. We derived a precise distance to the system of 2.14 ± 0.06 kpc (stat.) ± 0.05 (syst.) which places the star in the Sagittarius-Carina arm. The Galactic rotational velocity of the star is Θs = 258 ± 26 km s-1 assuming Θ0 = 238 km s-1. A comparison with PARSEC isochrones places the system at the early phase of core helium burning with an age of slightly larger than 100 million years. The effect of overshooting on stellar evolutionary tracks was explored using the MESA star code.
Drew, Simon C; Young, Charles G; Hanson, Graeme R
2007-04-01
The electron paramagnetic resonance spin Hamiltonian parameters of mononuclear thiomolybdenyl complexes based upon the tris(pyrazolyl)borate ligand, together with their molybdenyl analogues, are calculated using density functional theory. The electronic g matrix and 95Mo hyperfine matrix are calculated as second-order response properties from the coupled-perturbed Kohn-Sham equations. The scalar relativistic zero-order regular approximation (ZORA) is used with an all-electron basis and an accurate mean-field spin-orbit operator which includes all one- and two-electron terms. The principal values and relative orientations of the g and A interaction matrices obtained from the experimental spectra in a previous EPR study are compared with those obtained from unrestricted Kohn-Sham calculations at the BP86 and B3LYP level, and the latter are found to be in good quantitative agreement. A quasi-restricted approach is used to analyze the influence of the various molecular orbitals on g and A. In all complexes the ground state magnetic orbital is dX2-Y2-based and the orientation of the A matrix is directly related to the orientation of this orbital. The largest single contribution to the orientation of the g matrix arises from the spin-orbit coupling of the dYZ-based lowest-unoccupied molecular orbital into the ground state. A number of smaller, cumulative charge-transfer contributions augment the d-d contributions. A comparison of the theoretical EPR parameters obtained using both crystallographic and gas-phase geometry-optimized structures of Tp*MoO(bdt) (Tp* = hydrotris(3,5-dimethylpyrazol-1-yl)borate, bdt = 1,2-benzenedithiolate) suggests a correspondence between the metal-dithiolate fold angle and the angle of noncoincidence between g and A. PMID:17305330
Navratil, P; Caurier, E
2003-10-14
The authors calculate properties of A = 6 system using the accurate charge-dependent nucleon-nucleon (NN) potential at fourth order of chiral perturbation theory. By application of the ab initio no-core shell model (NCSM) and a variational calculation in the harmonic oscillator basis with basis size up to 16 {h_bar}{Omega} they obtain the {sup 6}Li binding energy of 28.5(5) MeV and a converged excitation spectrum. Also, they calculate properties of {sup 10}B using the same NN potential in a basis space of up to 8 {h_bar}{Omega}. The results are consistent with results obtained by standard accurate NN potentials and demonstrate a deficiency of Hamiltonians consisting of only two-body terms. At this order of chiral perturbation theory three-body terms appear. It is expected that inclusion of such terms in the Hamiltonian will improve agreement with experiment.
NASA Astrophysics Data System (ADS)
Orkin, V. L.; Khamaganov, V. G.; Martynova, L. E.; Kurylo, M. J.
2012-12-01
The emissions of halogenated (Cl, Br containing) organics of both natural and anthropogenic origin contribute to the balance of and changes in the stratospheric ozone concentration. The associated chemical cycles are initiated by the photochemical decomposition of the portion of source gases that reaches the stratosphere. Reactions with hydroxyl radicals and photolysis are the main processes dictating the compound lifetime in the troposphere and release of active halogen in the stratosphere for a majority of halogen source gases. Therefore, the accuracy of photochemical data is of primary importance for the purpose of comprehensive atmospheric modeling and for simplified kinetic estimations of global impacts on the atmosphere, such as in ozone depletion (i.e., the Ozone Depletion Potential, ODP) and climate change (i.e., the Global Warming Potential, GWP). The sources of critically evaluated photochemical data for atmospheric modeling, NASA/JPL Publications and IUPAC Publications, recommend uncertainties within 10%-60% for the majority of OH reaction rate constants with only a few cases where uncertainties lie at the low end of this range. These uncertainties can be somewhat conservative because evaluations are based on the data from various laboratories obtained during the last few decades. Nevertheless, even the authors of the original experimental works rarely estimate the total combined uncertainties of the published OH reaction rate constants to be less than ca. 10%. Thus, uncertainties in the photochemical properties of potential and current atmospheric trace gases obtained under controlled laboratory conditions still may constitute a major source of uncertainty in estimating the compound's environmental impact. One of the purposes of the presentation is to illustrate the potential for obtaining accurate laboratory measurements of the OH reaction rate constant over the temperature range of atmospheric interest. A detailed inventory of accountable sources of
Liu, Hui; Shi, Deheng; Sun, Jinfeng; Zhu, Zunlue; Shulin, Zhang
2014-04-24
The potential energy curves (PECs) of 54 spin-orbit states generated from the 22 electronic states of O2 molecule are investigated for the first time for internuclear separations from about 0.1 to 1.0nm. Of the 22 electronic states, the X(3)Σg(-), A(')(3)Δu, A(3)Σu(+), B(3)Σu(-), C(3)Πg, a(1)Δg, b(1)Σg(+), c(1)Σu(-), d(1)Πg, f(1)Σu(+), 1(5)Πg, 1(3)Πu, 2(3)Σg(-), 1(5)Σu(-), 2(1)Σu(-) and 2(1)Δg are found to be bound, whereas the 1(5)Σg(+), 2(5)Σg(+), 1(1)Πu, 1(5)Δg, 1(5)Πu and 2(1)Πu are found to be repulsive ones. The B(3)Σu(-) and d(1)Πg states possess the double well. And the 1(3)Πu, C(3)Πg, A'(3)Δu, 1(5)Δg and 2(5)Σg(+) states are the inverted ones when the spin-orbit coupling is included. The PEC calculations are done by the complete active space self-consistent field (CASSCF) method, which is followed by the internally contracted multireference configuration interaction (icMRCI) approach with the Davidson correction. Core-valence correlation and scalar relativistic corrections are taken into account. The convergence of present calculations is evaluated with respect to the basis set and level of theory. The vibrational properties are discussed for the 1(5)Πg, 1(3)Πu, d(1)Πg and 1(5)Σu(-) states and for the second well of the B(3)Σu(-) state. The spin-orbit coupling effect is accounted for by the state interaction method with the Breit-Pauli Hamiltonian. The PECs of all the electronic states and spin-orbit states are extrapolated to the complete basis set limit. The spectroscopic parameters are obtained, and compared with available experimental and other theoretical results. Analyses demonstrate that the spectroscopic parameters reported here can be expected to be reliably predicted ones. The conclusion is obtained that the effect of spin-orbit coupling on the spectroscopic parameters are small almost for all the electronic states involved in this paper except for the 1(5)Σu(-), 1(5)Πg and 1(3)Πu. PMID:24486866
NASA Astrophysics Data System (ADS)
Shi, De-Heng; Liu, Qionglan; Yu, Wei; Sun, Jinfeng; Zhu, Zunlue
2014-05-01
The potential energy curves (PECs) of 23 Ω states generated from the 12 electronic states (X1 Σ +, 21 Σ +, 11 Σ -, 11 Π, 21 Π, 11 Δ, 13 Σ +, 23 Σ +, 13 Σ -, a3 Π, 23 Π and 13 Δ) are studied for the first time. All the states correlate to the first dissociation channel of the SiBr+ cation. Of these electronic states, the 23 Σ + is the repulsive one without the spin-orbit coupling, whereas it becomes the bound one with the spin-orbit coupling added. On the one hand, without the spin-orbit coupling, the 11 Π, 21 Π and 23 Π are the rather weakly bound states, and only the 11 Π state possesses the double well; on the other hand, with the spin-orbit coupling included, the a3 Π and 11 Π states possess the double well, and the 13 Σ + and 13 Σ - are the inverted states. The PECs are calculated by the CASSCF method, which is followed by the internally contracted MRCI approach with the Davidson modification. Scalar relativistic correction is calculated by the third-order Douglas-Kroll Hamiltonian approximation with a cc-pVTZ-DK basis set. Core-valence correlation correction is included with a cc-pCVTZ basis set. The spin-orbit coupling is accounted for by the state interaction method with the Breit-Pauli Hamiltonian using the all-electron aug-cc-pCVTZ basis set. All the PECs are extrapolated to the complete basis set limit. The variation with internuclear separation of the spin-orbit coupling constant is discussed in brief. The spectroscopic parameters are evaluated for the 11 bound electronic states and the 23 bound Ω states, and are compared with available measurements. Excellent agreement has been found between the present results and the experimental data. It demonstrates that the spectroscopic parameters reported here can be expected to be reliably predicted ones. The Franck-Condon factors and radiative lifetimes of the transitions from the a3 Π 0 + and a3 Π 1 states to the X1 Σ + 0+ state are calculated for several low vibrational levels, and
NASA Astrophysics Data System (ADS)
Long, Xu; Ge, Fei; Wang, Lei; Hong, Youshi
2009-06-01
This paper investigates the effects of structure parameters on dynamic responses of submerged floating tunnel (SFT) under hydrodynamic loads. The structure parameters includes buoyancy-weight ratio (BWR), stiffness coefficients of the cable systems, tunnel net buoyancy and tunnel length. First, the importance of structural damp in relation to the dynamic responses of SFT is demonstrated and the mechanism of structural damp effect is discussed. Thereafter, the fundamental structure parameters are investigated through the analysis of SFT dynamic responses under hydrodynamic loads. The results indicate that the BWR of SFT is a key structure parameter. When BWR is 1.2, there is a remarkable trend change in the vertical dynamic response of SFT under hydrodynamic loads. The results also indicate that the ratio of the tunnel net buoyancy to the cable stiffness coefficient is not a characteristic factor affecting the dynamic responses of SFT under hydrodynamic loads.
Effective parameters in beam acoustic metamaterials based on energy band structures
NASA Astrophysics Data System (ADS)
Jing, Li; Wu, Jiu Hui; Guan, Dong; Hou, Mingming; Kuan, Lu; Shen, Li
2016-07-01
We present a method to calculate the effective material parameters of beam acoustic metamaterials. The effective material parameters of a periodic beam are calculated as an example. The dispersion relations and energy band structures of this beam are calculated. Subsequently, the effective material parameters of the beam are investigated by using the energy band structures. Then, the modal analysis and transmission properties of the beams with finite cells are simulated in order to confirm the correctness of effective approximation. The results show that the periodic beam can be equivalent to the homogeneous beam with dynamic effective material parameters in passband.
Sharapov, Vladimir A; Mandelshtam, Vladimir A
2007-10-18
We consider systems undergoing very-low-temperature solid-solid transitions associated with minima of similar energy but different symmetry, and separated by a high potential barrier. In such cases the well-known "broken-ergodicity" problem is often difficult to overcome, even using the most advanced Monte Carlo (MC) techniques, including the replica exchange method (REM). The methodology that we develop in this paper is suitable for the above specified cases and is numerically accurate and efficient. It is based on a new MC move implemented within the REM framework, in which trial points are generated analytically using an auxiliary harmonic superposition system that mimics well the true system at low temperatures. Due to the new move, the low-temperature random walks are able to frequently switch the relevant potential energy funnels leading to an efficient sampling. Numerically accurate results are obtained for a number of Lennard-Jones clusters, including those that have so far been treated only by the harmonic superposition approximation (HSA). The latter is believed to provide good estimates for low-temperature equilibrium properties but is manifestly uncontrollable and is difficult to validate. The present results provide a good test for the HSA and demonstrate its reliability, particularly for estimation of the solid-solid transition temperatures in most cases considered. PMID:17685597
Li, Liqi; Cui, Xiang; Yu, Sanjiu; Zhang, Yuan; Luo, Zhong; Yang, Hua; Zhou, Yue; Zheng, Xiaoqi
2014-01-01
Protein structure prediction is critical to functional annotation of the massively accumulated biological sequences, which prompts an imperative need for the development of high-throughput technologies. As a first and key step in protein structure prediction, protein structural class prediction becomes an increasingly challenging task. Amongst most homological-based approaches, the accuracies of protein structural class prediction are sufficiently high for high similarity datasets, but still far from being satisfactory for low similarity datasets, i.e., below 40% in pairwise sequence similarity. Therefore, we present a novel method for accurate and reliable protein structural class prediction for both high and low similarity datasets. This method is based on Support Vector Machine (SVM) in conjunction with integrated features from position-specific score matrix (PSSM), PROFEAT and Gene Ontology (GO). A feature selection approach, SVM-RFE, is also used to rank the integrated feature vectors through recursively removing the feature with the lowest ranking score. The definitive top features selected by SVM-RFE are input into the SVM engines to predict the structural class of a query protein. To validate our method, jackknife tests were applied to seven widely used benchmark datasets, reaching overall accuracies between 84.61% and 99.79%, which are significantly higher than those achieved by state-of-the-art tools. These results suggest that our method could serve as an accurate and cost-effective alternative to existing methods in protein structural classification, especially for low similarity datasets. PMID:24675610
Impact of Post Hoc Measurement Model Overspecification on Structural Parameter Integrity
ERIC Educational Resources Information Center
Fan, Weihua; Hancock, Gregory R.
2006-01-01
In the common two-step structural equation modeling process, modifications are routinely made to the measurement portion of the model prior to assessing structural relations. The effect of such measurement model modifications on the structural parameter estimates, however, is not well known and is the subject of the current investigation. For a…
Effect of UV light on different structural and transport parameters of cellophane membranes
Benavente, J.; Vazquez, M.I.; De Abajo, J.
1996-01-01
A comparative study of UV light influence on structural and transport parameters of cellophane membranes was made. Changes in the chemical structure and electrical behavior of cellophane membranes were considered by determining the hydraulic permeability, salt diffusion coefficient, and resistance values, as well as some geometrical parameters, for an untreated membrane and two differently UV-treated cellophane membranes. Differences in the characteristic parameters for the three samples showed that radiation mainly affected the membrane structure, while only small changes in membrane electrical behavior were determined.
Nägele, Edgar; Moritz, Ralf
2005-10-01
Today, it is necessary to identify relevant compounds appearing in discovery and development of new drug substances in the pharmaceutical industry. For that purpose, the measurement of accurate molecular mass and empirical formula calculation is very important for structure elucidation in addition to other available analytical methods. In this work, the identification and confirmation of degradation products in a finished dosage form of the antibiotic drug amoxicillin obtained under stress conditions will be demonstrated. Structure elucidation is performed utilizing liquid chromatography (LC) ion trap MS/MS and MS3 together with accurate mass measurement of the molecular ions and of the collision induced dissociation (CID) fragments by liquid chromatography electro spray ionization time-of-flight mass spectrometry (LC/ESI-TOF). PMID:16099170
Orenstein, Yaron; Wang, Yuhao; Berger, Bonnie
2016-01-01
Motivation: Protein–RNA interactions, which play vital roles in many processes, are mediated through both RNA sequence and structure. CLIP-based methods, which measure protein–RNA binding in vivo, suffer from experimental noise and systematic biases, whereas in vitro experiments capture a clearer signal of protein RNA-binding. Among them, RNAcompete provides binding affinities of a specific protein to more than 240 000 unstructured RNA probes in one experiment. The computational challenge is to infer RNA structure- and sequence-based binding models from these data. The state-of-the-art in sequence models, Deepbind, does not model structural preferences. RNAcontext models both sequence and structure preferences, but is outperformed by GraphProt. Unfortunately, GraphProt cannot detect structural preferences from RNAcompete data due to the unstructured nature of the data, as noted by its developers, nor can it be tractably run on the full RNACompete dataset. Results: We develop RCK, an efficient, scalable algorithm that infers both sequence and structure preferences based on a new k-mer based model. Remarkably, even though RNAcompete data is designed to be unstructured, RCK can still learn structural preferences from it. RCK significantly outperforms both RNAcontext and Deepbind in in vitro binding prediction for 244 RNAcompete experiments. Moreover, RCK is also faster and uses less memory, which enables scalability. While currently on par with existing methods in in vivo binding prediction on a small scale test, we demonstrate that RCK will increasingly benefit from experimentally measured RNA structure profiles as compared to computationally predicted ones. By running RCK on the entire RNAcompete dataset, we generate and provide as a resource a set of protein–RNA structure-based models on an unprecedented scale. Availability and Implementation: Software and models are freely available at http://rck.csail.mit.edu/ Contact: bab@mit.edu Supplementary information
Bitsie, Fernando; Jensen, Brian D.; de Boer, Maarten
1999-07-15
We have designed, fabricated, tested and modeled a first generation small area test structure for MEMS fracture studies by electrostatic rather than mechanical probing. Because of its small area, this device has potential applications as a lot monitor of strength or fatigue of the MEMS structural material. By matching deflection versus applied voltage data to a 3-D model of the test structure, we develop high confidence that the local stresses achieved in the gage section are greater than 1 GPa. Brittle failure of the polycrystalline silicon was observed.
NASA Technical Reports Server (NTRS)
Hyland, D. C.
1983-01-01
A stochastic structural control model is described. In contrast to the customary deterministic model, the stochastic minimum data/maximum entropy model directly incorporates the least possible a priori parameter information. The approach is to adopt this model as the basic design model, thus incorporating the effects of parameter uncertainty at a fundamental level, and design mean-square optimal controls (that is, choose the control law to minimize the average of a quadratic performance index over the parameter ensemble).
The computer program SPARC (SPARC Performs Automated Reasoning in Chemistry) has been under development for several years to estimate physical properties and chemical reactivity parameters of organic compounds strictly from molecular structure. SPARC uses computational algorithms...
Parameter identification methods for improving structural dynamic models. Ph.D. Thesis
NASA Technical Reports Server (NTRS)
Lawrence, Charles
1988-01-01
There is an increasing need to develop Parameter Identification methods for improving structural dynamic models, based on the inability of engineers to produce mathematical models which correlate with experimental data. This research explores the efficiency of combining Component Mode Synthesis (substructuring) methods with Parameter Identification procedures in order to improve analytical modeling of structural components and their connections. Improvements are computed in terms of physical stiffness and damping parameters in order that the physical characteristics of the model can be better understood. Connections involving both viscous and friction damping are investigated. Substructuring methods are utilized to reduce the complexity of the identification problem. Component and inter-component structural connection properties are evaluated and identified independently, thus simplifying the identification problem. It is shown that modal test data is effective for identifying modeling problems associated with structural components, and for determining the stiffness and damping properties of intercomponent connections. In general, Parameter Identification is improved when greater quantities of experimental data are available.
Distributions of electric parameters in MOS structures on 3C-SiC substrate
NASA Astrophysics Data System (ADS)
Piskorski, Krzysztof; Przewlocki, Henryk M.; Esteve, Romain; Bakowski, Mietek
2013-02-01
In this work studies of some electrical parameters of the MOS structure based on 3C-SiC substrate are presented. The effective contact potential difference ϕ MS , the barrier height at the gate-dielectric interface E BG and the flat-band in semiconductor voltage V FB were measured using several electric and photoelectric techniques. Values of these parameters obtained on structures with different gate areas decrease monotonically with increasing parameter R, defined as the ratio of the gate perimeter to the gate area. Such behavior confirmed results obtained on MOS structures on silicon substrate and also supported our hypothesis that the mechanical stress in the dielectric layer under the metal gate causes non uniform distribution of some parameters over the gate area of MOS structure.
Distributions of electric parameters in MOS structures on 3C-SiC substrate
NASA Astrophysics Data System (ADS)
Piskorski, Krzysztof; Przewlocki, Henryk; Esteve, Romain; Bakowski, Mietek
2013-02-01
In this work studies of some electrical parameters of the MOS structure based on 3C-SiC substrate are presented. The effective contact potential difference ϕMS, the barrier height at the gate-dielectric interface E BG and the flat-band in semiconductor voltage V FB were measured using several electric and photoelectric techniques. Values of these parameters obtained on structures with different gate areas decrease monotonically with increasing parameter R, defined as the ratio of the gate perimeter to the gate area. Such behavior confirmed results obtained on MOS structures on silicon substrate and also supported our hypothesis that the mechanical stress in the dielectric layer under the metal gate causes non uniform distribution of some parameters over the gate area of MOS structure.
Standard Errors of Estimated Latent Variable Scores with Estimated Structural Parameters
ERIC Educational Resources Information Center
Hoshino, Takahiro; Shigemasu, Kazuo
2008-01-01
The authors propose a concise formula to evaluate the standard error of the estimated latent variable score when the true values of the structural parameters are not known and must be estimated. The formula can be applied to factor scores in factor analysis or ability parameters in item response theory, without bootstrap or Markov chain Monte…
Bias-Corrected Estimation of Noncentrality Parameters of Covariance Structure Models
ERIC Educational Resources Information Center
Raykov, Tenko
2005-01-01
A bias-corrected estimator of noncentrality parameters of covariance structure models is discussed. The approach represents an application of the bootstrap methodology for purposes of bias correction, and utilizes the relation between average of resample conventional noncentrality parameter estimates and their sample counterpart. The…
Determination of Structural Parameters of Thin-Film Photocatalytic Materials by BDS
NASA Astrophysics Data System (ADS)
Korte, Dorota; Franko, Mladen
2015-09-01
A method for determination of structural parameters of some thin-film photocatalytic materials is presented. The analysis was based on the material's thermal parameter dependences on its surface structure or porosity and was thus performed by the use of beam deflection spectroscopy (BDS) supported by theoretical analysis made in the framework of complex geometrical optics. The results obtained by BDS were than compared with those received on the basis of AFM and SEM measurements and found to be in good agreement.
NASA Technical Reports Server (NTRS)
Cooksy, A. L.; Saykally, R. J.; Brown, J. M.; Evenson, K. M.
1986-01-01
Accurate values are presented for the fine-structure intervals in the 3P ground state of neutral atomic C-12 and C-13 as obtained from laser magnetic resonance spectroscopy. The rigorous analysis of C-13 hyperfine structure, the measurement of resonant fields for C-12 transitions at several additional far-infrared laser frequencies, and the increased precision of the C-12 measurements, permit significant improvement in the evaluation of these energies relative to earlier work. These results will expedite the direct and precise measurement of these transitions in interstellar sources and should assist in the determination of the interstellar C-12/C-13 abundance ratio.
Geerligs, Linda; Cam-Can; Henson, Richard N
2016-07-15
Studies of brain-wide functional connectivity or structural covariance typically use measures like the Pearson correlation coefficient, applied to data that have been averaged across voxels within regions of interest (ROIs). However, averaging across voxels may result in biased connectivity estimates when there is inhomogeneity within those ROIs, e.g., sub-regions that exhibit different patterns of functional connectivity or structural covariance. Here, we propose a new measure based on "distance correlation"; a test of multivariate dependence of high dimensional vectors, which allows for both linear and non-linear dependencies. We used simulations to show how distance correlation out-performs Pearson correlation in the face of inhomogeneous ROIs. To evaluate this new measure on real data, we use resting-state fMRI scans and T1 structural scans from 2 sessions on each of 214 participants from the Cambridge Centre for Ageing & Neuroscience (Cam-CAN) project. Pearson correlation and distance correlation showed similar average connectivity patterns, for both functional connectivity and structural covariance. Nevertheless, distance correlation was shown to be 1) more reliable across sessions, 2) more similar across participants, and 3) more robust to different sets of ROIs. Moreover, we found that the similarity between functional connectivity and structural covariance estimates was higher for distance correlation compared to Pearson correlation. We also explored the relative effects of different preprocessing options and motion artefacts on functional connectivity. Because distance correlation is easy to implement and fast to compute, it is a promising alternative to Pearson correlations for investigating ROI-based brain-wide connectivity patterns, for functional as well as structural data. PMID:27114055
Barone, Vincenzo; Biczysko, Malgorzata; Bloino, Julien; Cimino, Paola; Penocchio, Emanuele; Puzzarini, Cristina
2015-09-01
The structures and relative stabilities as well as the rotational and vibrational spectra of the three low-energy conformers of pyruvic acid (PA) have been characterized using a state-of-the-art quantum-mechanical approach designed for flexible molecules. By making use of the available experimental rotational constants for several isotopologues of the most stable PA conformer, Tc-PA, the semiexperimental equilibrium structure has been derived. The latter provides a reference for the pure theoretical determination of the equilibrium geometries for all conformers, thus confirming for these structures an accuracy of 0.001 Å and 0.1 deg for bond lengths and angles, respectively. Highly accurate relative energies of all conformers (Tc-, Tt-, and Ct-PA) and of the transition states connecting them are provided along with the thermodynamic properties at low and high temperatures, thus leading to conformational enthalpies accurate to 1 kJ mol(-1). Concerning microwave spectroscopy, rotational constants accurate to about 20 MHz are provided for the Tt- and Ct-PA conformers, together with the computed centrifugal-distortion constants and dipole moments required to simulate their rotational spectra. For Ct-PA, vibrational frequencies in the mid-infrared region accurate to 10 cm(-1) are reported along with theoretical estimates for the transitions in the near-infrared range, and the corresponding infrared spectrum including fundamental transitions, overtones, and combination bands has been simulated. In addition to the new data described above, theoretical results for the Tc- and Tt-PA conformers are compared with all available experimental data to further confirm the accuracy of the hybrid coupled-cluster/density functional theory (CC/DFT) protocol applied in the present study. Finally, we discuss in detail the accuracy of computational models fully based on double-hybrid DFT functionals (mainly at the B2PLYP/aug-cc-pVTZ level) that avoid the use of very expensive CC
Wang, Yong; Goh, Wang Ling; Chai, Kevin T-C; Mu, Xiaojing; Hong, Yan; Kropelnicki, Piotr; Je, Minkyu
2016-04-01
The parasitic effects from electromechanical resonance, coupling, and substrate losses were collected to derive a new two-port equivalent-circuit model for Lamb wave resonators, especially for those fabricated on silicon technology. The proposed model is a hybrid π-type Butterworth-Van Dyke (PiBVD) model that accounts for the above mentioned parasitic effects which are commonly observed in Lamb-wave resonators. It is a combination of interdigital capacitor of both plate capacitance and fringe capacitance, interdigital resistance, Ohmic losses in substrate, and the acoustic motional behavior of typical Modified Butterworth-Van Dyke (MBVD) model. In the case studies presented in this paper using two-port Y-parameters, the PiBVD model fitted significantly better than the typical MBVD model, strengthening the capability on characterizing both magnitude and phase of either Y11 or Y21. The accurate modelling on two-port Y-parameters makes the PiBVD model beneficial in the characterization of Lamb-wave resonators, providing accurate simulation to Lamb-wave resonators and oscillators. PMID:27131699
NASA Astrophysics Data System (ADS)
Wang, Yong; Goh, Wang Ling; Chai, Kevin T.-C.; Mu, Xiaojing; Hong, Yan; Kropelnicki, Piotr; Je, Minkyu
2016-04-01
The parasitic effects from electromechanical resonance, coupling, and substrate losses were collected to derive a new two-port equivalent-circuit model for Lamb wave resonators, especially for those fabricated on silicon technology. The proposed model is a hybrid π-type Butterworth-Van Dyke (PiBVD) model that accounts for the above mentioned parasitic effects which are commonly observed in Lamb-wave resonators. It is a combination of interdigital capacitor of both plate capacitance and fringe capacitance, interdigital resistance, Ohmic losses in substrate, and the acoustic motional behavior of typical Modified Butterworth-Van Dyke (MBVD) model. In the case studies presented in this paper using two-port Y-parameters, the PiBVD model fitted significantly better than the typical MBVD model, strengthening the capability on characterizing both magnitude and phase of either Y11 or Y21. The accurate modelling on two-port Y-parameters makes the PiBVD model beneficial in the characterization of Lamb-wave resonators, providing accurate simulation to Lamb-wave resonators and oscillators.
NASA Astrophysics Data System (ADS)
Mita, Akira
2016-04-01
The importance of the structural health monitoring system for tall buildings is now widely recognized by at least structural engineers and managers at large real estate companies to ensure the structural safety immediately after a large earthquake and appeal the quantitative safety of buildings to potential tenants. Some leading real estate companies decided to install the system into all tall buildings. Considering this tendency, a pilot project for the west area of Shinjuku Station supported by the Japan Science and Technology Agency was started by the author team to explore a possibility of using the system to provide safe spaces for commuters and residents. The system was installed into six tall buildings. From our experience, it turned out that viewing only from technological aspects was not sufficient for the system to be accepted and to be really useful. Safe spaces require not only the structural safety but also the soundness of key functions of the building. We need help from social scientists, medical doctors, city planners etc. to further improve the integrity of the system.
Parameter estimation in a structural acoustic system with fully nonlinear coupling conditions
NASA Technical Reports Server (NTRS)
Banks, H. T.; Smith, Ralph C.
1994-01-01
A methodology for estimating physical parameters in a class of structural acoustic systems is presented. The general model under consideration consists of an interior cavity which is separated from an exterior noise source by an enclosing elastic structure. Piezoceramic patches are bonded to or embedded in the structure; these can be used both as actuators and sensors in applications ranging from the control of interior noise levels to the determination of structural flaws through nondestructive evaluation techniques. The presence and excitation of patches, however, changes the geometry and material properties of the structure as well as involves unknown patch parameters, thus necessitating the development of parameter estimation techniques which are applicable in this coupled setting. In developing a framework for approximation, parameter estimation and implementation, strong consideration is given to the fact that the input operator is unbonded due to the discrete nature of the patches. Moreover, the model is weakly nonlinear. As a result of the coupling mechanism between the structural vibrations and the interior acoustic dynamics. Within this context, an illustrating model is given, well-posedness and approximations results are discussed and an applicable parameter estimation methodology is presented. The scheme is then illustrated through several numerical examples with simulations modeling a variety of commonly used structural acoustic techniques for systems excitations and data collection.
Correlation and assessment of structural airplane crash data with flight parameters at impact
NASA Technical Reports Server (NTRS)
Carden, H. D.
1982-01-01
Crash deceleration pulse data from a crash dynamics program on general aviation airplanes and from transport crash data were analyzed. Structural airplane crash data and flight parameters at impact were correlated. Uncoupled equations for the normal and longitudinal floor impulses in the cabin area of the airplane were derived, and analytical expressions for structural crushing during impact and horizontal slide out were also determined. Agreement was found between experimental and analytical data for general aviation and transport airplanes over a relatively wide range of impact parameter. Two possible applications of the impulse data are presented: a postcrash evaluation of crash test parameters and an assumed crash scenario.
Early Type Galaxies and Structural Parameters from ESO Public Survey KiDS
NASA Astrophysics Data System (ADS)
Roy, N.; Napolitano, N. R.; La Barbera, F.; Tortora, C.; Getman, F.; Radovich, M.; Capaccioli, M.
The Kilo Degree survey (KiDS) is a large-scale optical imaging survey carried out with the VLT Survey Telescope (VST), which is the ideal tool for galaxy evolution studies. We expect to observe millions of galaxies for which we extract the structural parameters in four wavebands (u, g, r and i). This sample will represent the largest dataset with measured structural parameters up to a redshift z = 0. 5. In this paper we will introduce the sample, and describe the 2D fitting procedure using the 2DPHOT environment and the validation of the parameters with an external catalog.
NASA Astrophysics Data System (ADS)
Bera, Mahua; Banerjee, Jayeta; Ray, Mina
2014-02-01
Metallic film thickness optimization in mono- and bimetallic plasmonic structures has been carried out in order to determine the correct device parameters. Different resonance parameters, such as reflectivity, phase, field enhancement, and the complex amplitude reflectance Argand diagram (CARAD), have been investigated for the proposed optimization procedure. Comparison of mono- and bimetallic plasmonic structures has been carried out in the context of these resonance parameters with simultaneous angular and spectral interrogation. Differential phase analysis has also been performed and its application to sensing has been discussed along with a proposed interferometric set-up.
Woon, D.E.; Dunning, T.H. Jr.; Peterson, K.A.
1996-04-01
Augmented correlation consistent basis sets of double (aug-cc-pVDZ), triple (aug-cc-pVTZ), and modified quadruple zeta (aug-cc-pVQZ{prime}) quality have been employed to describe the N{sub 2}{endash}HF potential energy surface at the Hartree{endash}Fock level and with single reference correlated wave functions including Mo/ller{endash}Plesset perturbation theory (MP2, MP3, MP4) and coupled cluster methods [CCSD, CCSD(T)]. The most accurate computed equilibrium binding energies {ital D}{sub {ital e}} are (with counterpoise correction) 810 cm{sup {minus}1} (MP4/aug-cc-pVQZ{prime}) and 788 cm{sup {minus}1} [CCSD(T)/aug-cc-pVQZ{prime}]. Estimated complete basis set limits of 814 cm{sup {minus}1} (MP4) and 793 cm{sup {minus}1} [CCSD(T)] indicate that the large basis set results are essentially converged. Harmonic frequencies and zero-point energies were determined through the aug-cc-pVTZ level. Combining the zero point energies computed at the aug-cc-pVTZ level with the equilibrium binding energies computed at the aug-cc-pVQZ{prime} level, we predict {ital D}{sub 0} values of 322 and 296 cm{sup {minus}1}, respectively, at the MP4 and CCSD(T) levels of theory. Using experimental anharmonic frequencies, on the other hand, the CCSD(T) value of {ital D}{sub 0} is increased to 415 cm{sup {minus}1}, in good agreement with the experimental value recently reported by Miller and co-workers, 398{plus_minus}2 cm{sup {minus}1}. {copyright} {ital 1996 American Institute of Physics.}
Li, Xiaowei; Liu, Taigang; Tao, Peiying; Wang, Chunhua; Chen, Lanming
2015-12-01
Structural class characterizes the overall folding type of a protein or its domain. Many methods have been proposed to improve the prediction accuracy of protein structural class in recent years, but it is still a challenge for the low-similarity sequences. In this study, we introduce a feature extraction technique based on auto cross covariance (ACC) transformation of position-specific score matrix (PSSM) to represent a protein sequence. Then support vector machine-recursive feature elimination (SVM-RFE) is adopted to select top K features according to their importance and these features are input to a support vector machine (SVM) to conduct the prediction. Performance evaluation of the proposed method is performed using the jackknife test on three low-similarity datasets, i.e., D640, 1189 and 25PDB. By means of this method, the overall accuracies of 97.2%, 96.2%, and 93.3% are achieved on these three datasets, which are higher than those of most existing methods. This suggests that the proposed method could serve as a very cost-effective tool for predicting protein structural class especially for low-similarity datasets. PMID:26460680
van Deelen, Rutger; Hasekamp, Otto P; Landgraf, Jochen
2007-01-10
We present what we believe to be a novel approach to simulating the spectral fine structure (<1 nm) in measurements of spectrometers such as the Global Ozone Monitoring Experiment (GOME). GOME measures the Earth's radiance spectra and daily solar irradiance spectra from which a reflectivity spectrum is commonly extracted. The high-frequency structures contained in such a spectrum are, apart from atmospheric absorption, caused by Raman scattering and by a shift between the solar irradiance and the Earth's radiance spectrum. Normally, an a priori high-resolution solar spectrum is used to simulate these structures. We present an alternative method in which all the required information on the solar spectrum is retrieved from the GOME measurements. We investigate two approaches for the spectral range of 390-400 nm. First, a solar spectrum is reconstructed on a fine spectral grid from the GOME solar measurement. This approach leads to undersampling errors of up to 0.5% in the modeling of the Earth's radiance spectra. Second, a combination of the solar measurement and one of the Earth's radiance measurement is used to retrieve a solar spectrum. This approach effectively removes the undersampling error and results in residuals close to the GOME measurement noise of 0.1%. PMID:17268571
NASA Astrophysics Data System (ADS)
Chen, Hanghui; Millis, Andrew J.
2016-05-01
We systematically compare predictions of various exchange correlation functionals for the structural and magnetic properties of perovskite Sr1 -xBaxMnO3 (0 ≤x ≤1 )—a representative class of multiferroic oxides. The local spin density approximation (LSDA) and spin-dependent generalized gradient approximation with Perdew-Burke-Ernzerhof parametrization (sPBE) make substantial different predictions for ferroelectric atomic distortions, tetragonality, and ground state magnetic ordering. Neither approximation quantitatively reproduces all the measured structural and magnetic properties of perovskite Sr0.5Ba0.5MnO3 . The spin-dependent generalized gradient approximation with Perdew-Burke-Ernzerhof revised for solids parametrization (sPBEsol) and the charge-only Perdew-Burke-Ernzerhof parametrized generalized gradient approximation with Hubbard U and Hund's J extensions both provide overall better agreement with measured structural and magnetic properties of Sr0.5Ba0.5MnO3 , compared to LSDA and sPBE. Using these two methods, we find that different from previous predictions, perovskite BaMnO3 has large Mn off-center displacements and is close to a ferromagnetic-to-antiferromagnetic phase boundary, making it a promising candidate to induce effective giant magnetoelectric effects and to achieve cross-field control of polarization and magnetism.