149 results on '"Global fitting"'
Search Results
2. The influence of the way of regression on the results obtained by the receptorial responsiveness method (RRM), a procedure to estimate a change in the concentration of a pharmacological agonist near the receptor.
- Author
-
Ovari, Ignac, Viczjan, Gabor, Erdei, Tamas, Takacs, Barbara, Tarjanyi, Vera, Zsuga, Judit, Szucs, Miklos, Szilvassy, Zoltan, Juhasz, Bela, and Gesztelyi, Rudolf
- Subjects
CURVE fitting ,ADENOSINES - Abstract
The receptorial responsiveness method (RRM) enables the estimation of a change in concentration of an (even degradable) agonist, near its receptor, via curve fitting to (at least) two concentration-effect (E/c) curves of a stable agonist. One curve should be generated before this change, and the other afterwards, in the same system. It follows that RRM yields a surrogate parameter ("cx") as the concentration of the stable agonist being equieffective with the change in concentration of the other agonist. However, regression can be conducted several ways, which can affect the accuracy, precision and ease-of-use. This study utilized data of previous ex vivo investigations. Known concentrations of stable agonists were estimated with RRM by performing individual (local) or global fitting, this latter with one or two model(s), using a logarithmic (logcx) or a nonlogarithmic (cx) parameter (the latter in a complex or in a simplified equation), with ordinary least-squares or robust regression, and with an "allat-once" or "pairwise" fitting manner. We found that the simplified model containing logcx was superior to all alternative models. The most complicated individual regression was the most accurate, followed closely by the moderately complicated two-model global regression and then by the easy-to-perform onemodel global regression. The two-model global fitting was the most precise, followed by the individual fitting (closely) and by the one-model global fitting (from afar). Pairwise fitting (two E/c curves at once) improved the estimation. Thus, the two-model global fitting, performed pairwise, and the individual fitting are recommended for RRM, using the simplified model containing logcx. [ABSTRACT FROM AUTHOR]
- Published
- 2024
- Full Text
- View/download PDF
3. Kinetics of the cross‐reaction of CH3O2 + HO2 radicals measured in the Highly Instrumented Reactor for Atmospheric Chemistry.
- Author
-
Østerstrøm, Freja F., Onel, Lavinia, Brennan, Alexander, Parr, Joseph M., Whalley, Lisa K., Seakins, Paul W., and Heard, Dwayne E.
- Subjects
- *
RADICALS (Chemistry) , *ATMOSPHERIC chemistry , *CROSS reactions (Immunology) , *CHEMICAL kinetics , *GAS phase reactions - Abstract
The sensitive Fluorescence Assay by Gas Expansion (FAGE) method has been used to detect methyl peroxy (CH3O2) and hydroperoxyl (HO2) radicals after their conversion by titration with excess NO to methoxy (CH3O) and hydroxyl (OH) radicals, respectively, to study the kinetics of the reaction of CH3O2 + HO2 radicals. The rate coefficient of the reaction was measured in the Highly Instrumented Reactor for Atmospheric Chemistry (HIRAC) at 1000 mbar of synthetic air at T = 268–344 K, selectively detecting both radicals. Using a numerical model to fit both CH3O2 and HO2 radical temporal decays globally at each temperature investigated, rate coefficients for the reaction have been obtained. The room temperature rate coefficient was found to be kCH3O2 +HO2(295 K) = (4.6 ± 0.7) × 10−12 molecule−1 cm3 s−1 (2σ errors) and the temperature dependence of the rate coefficient can be characterized in Arrhenius form by kCH3O2 + HO2(268 K < T < 344 K) = (5.1 ± 2.1) × 10−13 × exp((637 ± 121)/T) cm3 molecule−1 s−1. The rate coefficients obtained here are 14%–16% lower than the literature recommended values with an uncertainty which is reduced significantly compared to previous reports. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
4. The kinetic mechanism of DNA strand separation by high-fidelity DNA methyltransferase, CcrM
- Author
-
Konttinen, Olivia Rae
- Subjects
Biochemistry ,Bioinformatics ,Biophysics ,CcrM ,DNA methyltransferase ,DNA recognition ,Enzyme kinetics ,Global fitting ,Strand separation - Abstract
DNA methyltransferases are responsible for transcriptional regulation, cell cycle progression, DNA repair, DNA protection, tumor suppression, and several other important biological processes. Aberrant bacterial DNA methylation can lead to cell death and loss of protection against viral infection; in humans this leads to cancer, autoimmune diseases, metabolic disorders, and neurological disorders. Thus, DNA methyltransferases are common drug targets for cancer therapeutics and novel antibiotics.The conformational transitions in DNA and protein that govern recognition, substrate accessibility, and catalysis are fundamental to understanding the mechanisms that regulate biological processes. The bacterial N6-adenine cell-cycle regulated DNA methyltransferase, CcrM, is the first DNA methyltransferase shown to rely on a unique DNA recognition mechanism in which the DNA strands are separated and most recognition interactions appear to involve the target strand. Strand-separation is emerging as a novel DNA recognition mechanism but the underlying mechanisms and quantitative contribution of strand-separation to fidelity remain obscure for any enzyme (CRISPR-Cas9 and RNA polymerase sigma factor). This work uncovers the fundamental steps governing CcrM's DNA strand separation and high-fidelity DNA recognition mechanism. We relied on mutational analysis of highly conserved residues in the C-terminal domain, Loop-2B, Loop-45, and the active site to probe the function of structurally implicated protein moieties. We collected stopped-flow kinetic fluorescence to monitor transitions in DNA and protein and relied on rigorous global data fitting to understand the states that regulate catalysis in CcrM.We incorporated Pyrrolo-dC into cognate and noncognate DNA to monitor the kinetics of strand-separation and used tryptophan fluorescence to follow protein conformational changes. Both signals are biphasic and global fitting showed that the faster phase of DNA strand-separation was coincident with the protein conformational transition. Non-cognate sequences did not display strand-separation and methylation was reduced >300-fold, providing evidence that strand-separation is a major determinant of selectivity. Analysis of an R350A mutant (C-term domain) showed that the enzyme conformational step can occur without strand-separation, so the two events are uncoupled. A stabilizing role for the methyl-donor (SAM) is proposed; the cofactor interacts with a critical loop which is inserted between the DNA strands, thereby stabilizing the strand-separated conformation.Loops 2B and 45 are inserted between the strand-separated DNA interface. During strand-separation, residues within Loops 2B, 45, and 6E contact the target DNA strand that undergoes methylation. R44 and R129 (Loop-2B and Loop-45, respectively) when mutated to Alanine, disrupt strand-separation and are catalytically inactive. The highly conserved Loop-45 residue F125, which is positioned between the separated DNA strands, is also essential for maintaining the strand-separated intermediate; replacement of F125 with Alanine, Leucine, and Tryptophan results in various perturbations of strand-separation that are correlated to the bulkiness of the substituted residue. Global fitting for each mutant shows that generation and stabilization of DNA strand-separation are perturbed, providing a functional role for these loops in generating and stabilizing the strand-separated intermediate, which is essential for discrimination and catalysis.Employing a fluorescent adenine analog (6MAP) at the target position to monitor base flipping, we resolved that target adenine base flipping follows DNA strand separation and is followed by fast methylation and fast product DNA release. A W57F mutant (active site) displayed an unaltered rate of base flipping as monitored by 6MAP fluorescence but greatly reduced rate of methylation, showing that base-flipping and methylation can be uncoupled. In addition, single-stranded DNA bypasses the DNA strand separation step, while rates of base flipping measured by 6MAP fluorescence and DNA methylation are similar to dsDNA. Global data fitting for each model resolves that base flipping of the target adenine is the rate-limiting step in catalysis.The results presented here are broadly applicable to the study of other N6-adenine methyltransferases that contain the structural moieties implicated in strand-separation (Loop-2B, Loop-45, and the C-terminal domain), which are found widely dispersed across many bacterial phyla, including human and animal pathogens. Insights into CcrM’s mechanism of DNA strand-separation are likely to clarify strand-separation mechanisms for other enzymes such as CRISPR-Cas9 and RNA polymerase sigma factor. Additionally, the elevated understanding of CcrM’s strand-separation mechanism could be useful for the development of selective CcrM inhibitors as novel antibiotics.
- Published
- 2024
5. Quantitative NMR analysis of the mechanism and kinetics of chaperone Hsp104 action on amyloid-β42 aggregation and fibril formation.
- Author
-
Ghosh, Shreya, Tugarinov, Vitali, and Clore, G. Marius
- Subjects
- *
ALZHEIMER'S disease , *QUANTITATIVE research , *NUCLEAR forces (Physics) , *PEPTIDES , *MONOMERS - Abstract
The chaperone Hsp104, a member of the Hsp100/Clp family of translocases, prevents fibril formation of a variety of amyloidogenic peptides in a paradoxically substoichiometric manner. To understand the mechanism whereby Hsp104 inhibits fibril formation, we probed the interaction of Hsp104 with the Alzheimer's amyloid-ß42 (Aβ42) peptide using a variety of biophysical techniques. Hsp104 is highly effective at suppressing the formation of Thioflavin T (ThT) reactive mature fibrils that are readily observed by atomic force (AFM) and electron (EM) microscopies. Quantitative kinetic analysis and global fitting was performed on serially recorded 1H-15N correlation spectra to monitor the disappearance of Aβ42 monomers during the course of aggregation over a wide range of Hsp104 concentrations. Under the conditions employed (50 µM Aβ42 at 20 ℃), Aβ42 aggregation occurs by a branching mechanism: an irreversible on-pathway leading to mature fibrils that entails primary and secondary nucleation and saturating elongation; and a reversible off-pathway to form nonfibrillar oligomers, unreactive to ThT and too large to be observed directly by NMR, but too small to be visualized by AFM or EM. Hsp104 binds reversibly with nanomolar affinity to sparsely populated Aβ42 nuclei present in nanomolar concentrations, generated by primary and secondary nucleation, thereby completely inhibiting on-pathway fibril formation at substoichiometric ratios of Hsp104 to Aβ42 monomers. Tight binding to sparsely populated nuclei likely constitutes a general mechanism for substoichiometric inhibition of fibrillization by a variety of chaperones. Hsp104 also impacts off-pathway oligomerization but to a much smaller degree initially reducing and then increasing the rate of off-pathway oligomerization. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
6. Simultant: simultaneous curve fitting of functions and differential equations using analytical gradient calculations
- Author
-
Julius B. Kirkegaard
- Subjects
Data analysis ,Simultaneous fitting ,Global fitting ,Parameter sharing ,Differential equations ,Computer applications to medicine. Medical informatics ,R858-859.7 ,Biology (General) ,QH301-705.5 - Abstract
Abstract Background The initial step in comparing mathematical models to experimental data is to do a fit. This process can be complicated when either the mathematical models are not analytically solvable (e.g. because of nonlinear differential equations) or when the relation between data and models is complex (e.g. when some fitting parameters must be shared between many data sets). Results We introduce Simultant, a software package that allows complex fitting setups to be easily defined using a simple graphical user interface. Fitting functions can be defined directly as mathematical expressions or indirectly as the solution to specified ordinary differential equations. Analytical gradients of these functions, including the solution of differential equations, are automatically calculated to provide fast fitting even for functions with many parameters. The software enables easy definition of complex fitting setups in which parameters can be shared across both data sets and models to allow simultaneous fits to be performed. Conclusions Simultant exploits differentiable programming and simplifies modern fitting approaches in a unified graphical interface.
- Published
- 2022
- Full Text
- View/download PDF
7. Simultant: simultaneous curve fitting of functions and differential equations using analytical gradient calculations.
- Author
-
Kirkegaard, Julius B.
- Subjects
DIFFERENTIAL equations ,CURVE fitting ,NONLINEAR differential equations ,GRAPHICAL user interfaces ,MATHEMATICAL models - Abstract
Background: The initial step in comparing mathematical models to experimental data is to do a fit. This process can be complicated when either the mathematical models are not analytically solvable (e.g. because of nonlinear differential equations) or when the relation between data and models is complex (e.g. when some fitting parameters must be shared between many data sets). Results: We introduce Simultant, a software package that allows complex fitting setups to be easily defined using a simple graphical user interface. Fitting functions can be defined directly as mathematical expressions or indirectly as the solution to specified ordinary differential equations. Analytical gradients of these functions, including the solution of differential equations, are automatically calculated to provide fast fitting even for functions with many parameters. The software enables easy definition of complex fitting setups in which parameters can be shared across both data sets and models to allow simultaneous fits to be performed. Conclusions: Simultant exploits differentiable programming and simplifies modern fitting approaches in a unified graphical interface. [ABSTRACT FROM AUTHOR]
- Published
- 2022
- Full Text
- View/download PDF
8. Structured Stochastic Curve Fitting without Gradient Calculation.
- Author
-
Chen J
- Abstract
Optimization of parameters and hyperparameters is a general process for any data analysis. Because not all models are mathematically well-behaved, stochastic optimization can be useful in many analyses by randomly choosing parameters in each optimization iteration. Many such algorithms have been reported and applied in chemistry data analysis, but the one reported here is interesting to check out, where a naïve algorithm searches each parameter sequentially and randomly in its bounds. Then it picks the best for the next iteration. Thus, one can ignore irrational solution of the model itself or its gradient in parameter space., Competing Interests: Disclosure statement The author declares no competing financial interest.
- Published
- 2024
- Full Text
- View/download PDF
9. Intelligently optimized global analysis of time resolved spectra with particle swarm optimization.
- Author
-
Ma, Lin and Jiang, Lianlian
- Subjects
- *
PARTICLE swarm optimization , *TIME-resolved spectroscopy , *GLOBAL analysis (Mathematics) , *MANUAL labor , *GLOBAL optimization , *EXCITED states - Abstract
[Display omitted] • Particle swarm optimization automates global fitting of time resolved spectral data. • Particle swarm optimization eliminates manually setting initial parameters in conventional fitting. • Automation makes global analysis efficient and accessible for non-experts. Time-resolved spectroscopy, especially transient absorption spectroscopy (TAS), provides valuable insights to excited state dynamics. Analyzing TAS data involves fitting complex kinetic traces at various probe wavelengths using different rate equations. Conventional TAS global fitting methods require domain experts to establish physically valid models and provide good initial guesses to generate converged solutions. This poses challenges for non-experts who seek to utilize TAS, thus limiting its broader application and impact. To address this problem, we propose an intelligent optimization framework based on the particle swarm optimization (PSO) algorithm. In the proposed method, the PSO algorithm acts as the global fitting method to find the optimal values of the target variables or unknown parameters in the kinetics models. The target solution is optimized by iteratively updating candidate solutions with respect to an objective feedback signal. We demonstrated the effectiveness of the proposed PSO-based global fitting method with both synthetic and experimental datasets. The results show that our proposed method can successfully find the optimal target values in the global fitting process automatically, thus eliminating the iterative manual labor traditionally required. The proposed intelligent optimization framework provides a novel approach for automatic global fitting of TAS data, which significantly enhances the accessibility and utilization of the TAS methodology. [ABSTRACT FROM AUTHOR]
- Published
- 2024
- Full Text
- View/download PDF
10. Exploring charge transfer mechanisms and optical properties in vdW heterostructures of MoS2 and Bi2Se3 at nanoscale regime.
- Author
-
Chaudhary, Amit Kumar, Sharma, Prince, Rana, Archana, Jain, Sanyam, Saini, Saurabh K., Kumar, Kapil, Kuldeep, Singh, Rajiv K., Goswami, Lalit, Tanwar, Praveen K., Basheed, G.A., Gupta, Govind, and Kumar, Mahesh
- Subjects
- *
CHARGE transfer , *HETEROSTRUCTURES , *OPTICAL properties , *OPTOELECTRONIC devices , *CHARGE carriers , *ZETA potential - Abstract
We report here the two-dimensional layered material of MoS 2 /Bi 2 Se 3 vdW heterostructures formed by ultrasonication of different concentrations of exfoliated MoS 2 and Bi 2 Se 3 nanosheets. The vdW heterostructures show unique optical characteristics compared to individual exfoliated nanosheets. The enhanced lifetime of heterostructures explained the charge transfer with increasing efficacy for producing optoelectronic devices. Raman spectroscopy confirms the presence of intrinsic Raman active modes of Bi 2 Se 3 , MoS 2 , and respective modes of vdW heterostructures. TEM images show characteristic lattice planes in both materials and clear overlapping regimes, confirming successful heterojunction formation. Further, Zeta potential confirms the enhanced stability of the heterostructures formed while the significant PL quenching, strongly suggests the presence of charge transfer processes within the heterostructures. Ultrafast pump-probe spectroscopy further confirms the non-linear charge transfer and relaxation mechanism of the heterostructures. We employed Surface Xplorer and Glotaran software for model fitting and Glotaran's global fitting revealed unresolved charge carrier dynamics. [ABSTRACT FROM AUTHOR]
- Published
- 2024
- Full Text
- View/download PDF
11. 钙钛矿CsPbBr3中光生载流子的动力学行为.
- Author
-
宋宗鹏, 章琪, 蒋凌峰, 裴继红, and 阮双琢
- Abstract
Copyright of Journal of Shenzhen University Science & Engineering is the property of Editorial Department of Journal of Shenzhen University Science & Engineering and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)
- Published
- 2019
- Full Text
- View/download PDF
12. Global analysis of complex PELDOR time traces.
- Author
-
Rein, Stephan, Lewe, Philipp, Andrade, Susana L., Kacprzak, Sylwia, and Weber, Stefan
- Subjects
- *
ELECTRONS , *RESONANCE , *ELECTRON spin , *DATA analysis , *NOISE - Abstract
Graphical abstract Highlights • Sophisticated PELDOR/DEER data analysis software for Windows, macOS and Linux. • Global analysis of complex distance distributions. • Optimized data acquisition procedure saves measurement time. Abstract Pulsed electron–electron double resonance (PELDOR, alternatively called DEER for double electron–electron resonance) pulse sequences allow for the detection of echo decay curves that are modulated by dipole–dipole-coupling frequencies of interacting electron spins. With increasing distance between them, the echo decay needs to be monitored over a progressively extended time period. However, since the echo intensity typically falls off exponentially with increasing time, this might be problematic with respect to the minimum signal-to-noise ratio required for a sound data analysis. In this contribution we present the new PELDOR analysis tool GloPel (Glo bal analysis of PEL DOR data), an open-source Python-based application, that allows to extract improved-quality distance distributions from PELDOR data for which no ideal signal-to-noise ratio can be achieved for a very long observation window. By using Tikhonov regularization, GloPel allows for the simultaneous analysis of two time traces acquired for a sample in two different observation time windows, thus taking advantage of both, the typically high signal-to-noise ratio of the time trace acquired at early times of the echo decay, and the best possible background function fitted for the decay at later times, which is in most cases superimposed with considerable noise. In this way, short distances are not overseen in the higher noise of the longer time traces while long distances are not artificially shortened by limiting the observation time window of the experiment. Following our suggested data acquisition procedure, a significant reduction of the measurement time may also be achieved. [ABSTRACT FROM AUTHOR]
- Published
- 2018
- Full Text
- View/download PDF
13. Translation initiation in bacterial polysomes through ribosome loading on a standby site on a highly translated mRNA.
- Author
-
Andreeva, Irena, Belardinelli, Riccardo, and Rodnina, Marina V.
- Subjects
- *
TRANSLATION initiation factors (Biochemistry) , *RIBOSOME-inactivating proteins , *BACTERIAL RNA , *MRNA guanylyltransferase , *EXPONENTIAL functions - Abstract
During translation, consecutive ribosomes load on an mRNA and form a polysome. The first ribosome binds to a single-stranded mRNA region and moves toward the start codon, unwinding potential mRNA structures on theway. In contrast, the following ribosomes can dock at the start codon only when the first ribosome has vacated the initiation site. Here we show that loading of the second ribosome on a natural 38-nt-long 5' untranslated region of lpp mRNA, which codes for the outer membrane lipoprotein from Escherichia coli, takes place before the leading ribosome has moved away fromthe start codon. The rapid formation of this standby complex depends on the presence of ribosomal proteins S1/S2 in the leading ribosome. The early recruitment of the second ribosome to the standby site before translation by the leading ribosome and the tight coupling between translation elongation by the first ribosome and the accommodation of the second ribosome can contribute to high translational efficiency of the lpp mRNA. [ABSTRACT FROM AUTHOR]
- Published
- 2018
- Full Text
- View/download PDF
14. Highly Emissive Cerium(III) Thiocyanate Complexes Reevaluated Using a Silane‐Treated Quartz Cell
- Author
-
Yusuke Kuramochi, Shunsuke Sayama, Akiharu Satake, and Kotaro Ohminato
- Subjects
Inorganic Chemistry ,chemistry.chemical_compound ,Cerium ,Photoluminescence ,Thiocyanate ,Chemistry ,Inorganic chemistry ,Global fitting ,chemistry.chemical_element ,Silane ,Quartz - Published
- 2021
- Full Text
- View/download PDF
15. Second-sphere tuning of analogues for the ferric-hydroperoxoheme form of Mycobacterium tuberculosis MhuD.
- Author
-
Johnson, Kayla L., Graves, Amanda B., Eckhert, Kaitlyn, and Liptak, Matthew D.
- Subjects
- *
MYCOBACTERIUM tuberculosis , *HYDROGEN bonding , *HYDROXYL group , *CIRCULAR dichroism , *HOMOLYSIS - Abstract
Mycobacterium tuberculosis MhuD catalyzes the oxygenation of heme to mycobilin; experimental data presented here elucidates the novel hydroxylation reaction catalyzed by this enzyme. Analogues for the critical ferric–hydroperoxoheme (MhuD–heme–OOH) intermediate of this enzyme were characterized using UV/Vis absorption (Abs), circular dichroism (CD), and magnetic CD (MCD) spectroscopies. In order to extract electronic transition energies from these spectroscopic data, a novel global fitting model was developed for analysis of UV/Vis Abs, CD, and MCD data. A variant of MhuD was prepared, N7S, which weakens the affinity of heme-bound enzyme for a hydroperoxo analogue, azide, without significantly altering the protein secondary structure. Global fitting of spectroscopic data acquired in this study revealed that the second-sphere N7S substitution perturbs the electronic structure of two analogues for MhuD–heme–OOH: azide-inhibited MhuD (MhuD–heme–N 3) and cyanide-inhibited MhuD (MhuD–heme–CN). The ground state electronic structures of MhuD–heme–N 3 and MhuD–heme–CN were assessed using variable-temperature, variable-field MCD. Altogether, these data strongly suggest that there is a hydrogen bond between the Asn7 side-chain and the terminal oxygen of the hydroperoxo ligand in MhuD–heme–OOH. As discussed herein, this finding supports a novel hydroxylation reaction mechanism where the Asn7 side-chain guides a transient hydroxyl radical derived from homolysis of the O O bond in MhuD–heme–OOH to the β- or δ- meso carbon of the porphyrin ligand yielding β- or δ-meso-hydroxyheme, respectively. Synopsis: In this work, it is shown that Mycobacterium tuberculosis MhuD utilizes a hydrogen bond between Asn7 and the terminal oxygen of a hydroperoxo ligand to orient hydroperoxo along the β/δ-axis of heme. This orientation promotes regiospecific hydroxylation of heme to β- and δ- meso -hydroxyheme. [Display omitted] • Asn7 forms a hydrogen bond with the terminal atom of the distal ligand to heme in MhuD. • The hydrogen bond between Asn7 and the distal ligand orients the ligand within the active site. • A hydrogen bond with Asn7 perturbs the electronic structures of analogues for the ferric–hydroperoxoheme intermediate. • Asn7 guides a transient hydroxyl radical to the β- and δ-meso carbons of heme. • A novel global fitting method was developed for the analysis of UV/Vis Abs, CD, and MCD data. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
16. Global Analysis of 24 Rovibrational Bands of the Octad of the 76GeH4 Molecule
- Author
-
Elena Sergeevna Bekhtereva, N.I. Raspopova, Olga Vasilievna Gromova, Peter Sennikov, and M.A. Koshelev
- Subjects
Physics ,Global fitting ,General Physics and Astronomy ,Molecule ,Inverse ,Rotational–vibrational spectroscopy ,Atomic physics ,Standard deviation ,Variable (mathematics) - Abstract
A global fitting method has been developed for spherical top molecules. The octad of the GeH4 molecule comprising 24 interacting states has been investigated. By solving the inverse spectroscopic problem, positions of the band centers and rotational and resonant parameters are calculated. A total of 70 variable parameters allow 1184 initial experimental rovibrational energies to be reproduced with a standard deviation of 1.2·10–3 cm–1.
- Published
- 2021
- Full Text
- View/download PDF
17. Where are we with light sterile neutrinos?
- Author
-
Michael H. Shaevitz, Alejandro Diaz, Janet Conrad, Carlos Arguelles, and G. H. Collin
- Subjects
Physics ,Sterile neutrino ,Particle physics ,010308 nuclear & particles physics ,Bayesian probability ,Global fitting ,FOS: Physical sciences ,General Physics and Astronomy ,Parameter space ,01 natural sciences ,High Energy Physics - Experiment ,3. Good health ,High Energy Physics - Experiment (hep-ex) ,High Energy Physics - Phenomenology ,High Energy Physics - Phenomenology (hep-ph) ,Frequentist inference ,0103 physical sciences ,Energy spectrum ,High Energy Physics::Experiment ,Neutrino ,010306 general physics - Abstract
We review the status of searches for sterile neutrinos in the $\sim 1$ eV range, with an emphasis on the latest results from short baseline oscillation experiments and how they fit within sterile neutrino oscillation models. We present global fit results to a three-active-flavor plus one-sterile-flavor model (3+1), where we find an improvement of $\Delta \chi^2=35$ for 3 additional parameters compared to a model with no sterile neutrino. This is a 5$\sigma$ improvement, indicating that an effect that is like that of a sterile neutrino is highly preferred by the data. However we note that separate fits to the appearance and disappearance oscillation data sets within a 3+1 model do not show the expected overlapping allowed regions in parameter space. This "tension" leads us to explore two options: 3+2, where a second additional mass state is introduced, and a 3+1+decay model, where the $\nu_4$ state can decay to invisible particles. The 3+1+decay model, which is also motivated by improving compatibility with cosmological observations, yields the larger improvement, with a $\Delta \chi^2=8$ for 1 additional parameter beyond the 3+1 model, which is a $2.6\sigma$ improvement. Moreover the tension between appearance and disappearance experiments is reduced compared to 3+1, although disagreement remains. In these studies, we use a frequentist approach and also a Bayesean method of finding credible regions. With respect to this tension, we review possible problems with the global fitting method. We note multiple issues, including problems with reproducing the experimental results, especially in the case of experiments that do not provide adequate data releases. We discuss an unexpected 5 MeV excess, observed in the reactor flux energy spectrum, that may be affecting the oscillation interpretation of the short baseline reactor data. We emphasize the care that must be taken in mapping to the true neutrino energy in the case of oscillation experiments that are subject to multiple interaction modes and nuclear effects. We point to problems with the "Parameter-Goodness-of-Fit test" that is used to quantify the tension. Lastly, we point out that analyses presenting limits often receive less scrutiny that signals. While we provide a snapshot of the status of sterile neutrino searches today and global fits to their interpretation, we emphasize that this is a fast-moving field. We briefly review experiments that are expected to report new data in the immediate future. Lastly, we consider the 5-year horizon, where we propose that decay-at-rest neutrino sources are the best method of finally resolving the confusing situation., Comment: 52 pages, 36 figures. Corrected typos in Fig 12 and Table III, plus other minor corrections
- Published
- 2020
- Full Text
- View/download PDF
18. Excited-state dynamics of the 1Bu +, 3Ag −, and 1Bu − states in all-trans-spirilloxanthin as revealed by sub-5-fs time-resolved absorption spectroscopy
- Author
-
Kobayashi, Takayoshi, Nishimura, Kumiko, Rondonuwu, Ferdy S., Koyama, Yasushi, Castleman, A. W., Jr., editor, Toennies, J.P., editor, Zinth, W., editor, Kobayashi, Takayoshi, editor, Okada, Tadashi, editor, Kobayashi, Tetsuro, editor, Nelson, Keith A., editor, and De Silvestri, Sandro, editor
- Published
- 2005
- Full Text
- View/download PDF
19. Global fitting for high-accuracy multi-channel single-molecule localization
- Author
-
Ulf Matti, Yiming Li, Decheng Wu, Sheng Liu, Wei Shi, and Jonas Ries
- Subjects
Single molecule localization ,Computer science ,Microscopy ,Content (measure theory) ,Global fitting ,Graphics processing unit ,Biplane ,Algorithm ,Multi channel ,3d localization - Abstract
Multi-channel detection in single-molecule localization microscopy (SMLM) greatly increases information content for various biological applications. Here, we present globLoc, a graphics processing unit (GPU) based global fitting algorithm with flexible PSF modeling and parameter sharing, to extract maximum information from multi-channel single molecule data. We show, both in simulations and experiments, that global fitting can substantially improve the 3D localization precision for biplane and 4Pi SMLM and color assignment for ratiometric multicolor imaging.
- Published
- 2021
- Full Text
- View/download PDF
20. Alternative predictors in chaotic time series.
- Author
-
Alves, P.R.L., Duarte, L.G.S., and da Mota, L.A.C.P.
- Subjects
- *
TURBULENCE , *PROGRAMMING languages , *GNU General Public License (Free software license) , *OPEN source software , *SCALAR field theory - Abstract
In the scheme of reconstruction, non-polynomial predictors improve the forecast from chaotic time series. The algebraic manipulation in the Maple environment is the basis for obtaining of accurate predictors. Beyond the different times of prediction, the optional arguments of the computational routines optimize the running and the analysis of global mappings. New version program summary Program Title: LinMapTS Program Files doi: http://dx.doi.org/10.17632/pnhy9zymrp.1 Licensing provisions: GNU General Public License version 3 Programming language: Maple17 Journal reference of previous version: Comput. Phys. Comm. 207 (2016) 325 Does the new version supersede the previous version?: Yes Nature of problem: Time series analysis and improving forecast capability. Solution method: The method of solution is published in [1]. Restrictions: The routines employ the global variables { a i , b , X i } ; If more than 2000 vectors are employed in the global mapping the normality test is not applicable. Unusual features: The algebraic manipulation of the predictors improves the global forecast. Reasons for the new version: In the reconstruction’s scheme [2], the predictor in the global approach has a standard form [3]. From a time series { X ( 0 Δ t ) , X ( 1 Δ t ) ⋯ , X ( ( S − 1 ) Δ t ) } with S scalar quantities X ( t ) and a time interval Δ t , a state vector in the N -dimensional reconstructed phase space is (1) | x ( t ) 〉 ≐ [ X ( ( N − 1 ) T Δ t ) ⋮ X ( T Δ t ) X ( 0 ) ] . The time delay T is a parameter for the choice of the observables X ( ( N − k ) T Δ t ) ( k is an integer) that are available in the time series [4]. The predictors P ( | x r 〉 ) in the routine LinGfiTS of the package LinMapTS are linear combinations of the adjustment parameters. (2) P ( | x 〉 ) = ∑ i = 1 m a i φ i ( | x 〉 ) . Because this restricted form, the m -dimensional vector of parameters | a 〉 is a computational solution of a matrix equation. So the computational procedure requires small runtimes for the least squares minimization [1]. The routine generates only polynomial global maps, i.e. the functions φ ( | x 〉 ) can assume forms–in a reconstructed phase space with variables ( X 1 , X 2 , X 3 ) –like φ 1 = X 1 2 , φ 2 = X 1 3 X 2 2 X 3 and so on. The predictors do not admit terms such as sin ( X 1 X 2 X 3 ) or ln ( 1 + 1 X 1 X 2 X 3 ) . A prediction–denoted by X 1 P –is the result of the application of the global map (3) X 1 P = P ( | x P − 1 〉 ) , where P is the order of the last known observable in the time series. The principal focus in this new version is to extend the permissible functional forms for the global mappings. If non-polynomial terms take part in the predictors, the accuracy of the forecast can be improved. Here, the purpose is to offer the researcher best features when he intends to increase the predictions’ power from a chaotic time series . Another desired extension refers to the time of prediction. With the integer parameter τ , the future instant is given by t + τ Δ t . We have been presented this idea in the new version of the package TimeS [5]. But the runtime for generating polynomial maps is greater than the computational procedure LinGfiTS [1]. In this work, we apply 1 → τ in Eq. (3) . Then it is rewritten as (4) X 1 P = P ( | x P − τ 〉 ) . Summary of revisions: New optional arguments enable the selection of functional forms with different prediction times in the method of forecasting; Instructions to use the LinMapTS package (README.pdf), computational routines (LinMapTS.txt) and test file (LinMapTS.mw). In order to optimize the running of the programs LinGfiTS and ConfiTS , the arguments have now been rearranged. The current routines own optional arguments. Some of them are indispensable in the previous version. Nowadays, the command LinGfiTS requires only two arguments. The first is a list of reconstructed vectors–assigned as V –and the second is the order of the last vector present in the global mapping–assigned as f i n a l . The input necessary for the running of the procedure ConfiTS must have, in the following order: the global map–assigned as m a p –, the list V and the integer f i n a l . Below, we describe all arguments that take part in the new version of the package LinMapTS . • • Required arguments: – List of reconstructed vectors—assigned as V in this paper. – The vector that has, as its first component, the last known value of the time series—assigned as f i n a l in this paper. • Optional arguments: – Degree =
. This argument specifies the degree of the polynomial predictor. The default is 2. – Func = . This argument specifies the predictor for the global mapping. – Level . This argument selects the interval of the time series for the global mapping. The default is 5. – PT = . This argument specifies the value of the parameter τ . The default is 1. – Analysis = 1 . This argument is the necessary input for a graphical analysis of the global fitting. • • Required arguments: – The global map—assigned as m a p in this paper. – List of reconstructed vectors—assigned as V in this paper. – The vector that has, as its first component, the last known value of the time series—assigned as f i n a l in this paper. • Optional arguments: – Level . This argument selects the interval of the time series for the global mapping. The default is 5. – PT = . This argument specifies the value of the parameter τ . The default is 1. – Analysis = 1 . This argument is the necessary input for a graphical analysis of the residuals’ distribution and the applying of the normality test. (5) σ τ = ∑ j = 1 M ( X 1 j − ∑ i = 1 m a i φ i ( | x j − τ 〉 ) ) 2 M − 1 . The outputs remain unchanged from the original programs [1]. The routine LinGfiTS makes available a global map, whereas the procedure ConfiTS returns the expected deviation σ τ in the forecast. However, this statistical quantity now incorporates the parameter τ . Its formula (5) employs the M reconstructed vectors which take part in global fitting. As a first example of using the commands, the time parameter selected is τ = 2 . So it is necessary to include the optional argument PT=2 in the Maple prompt. Here, the file ts37.txt stores a chaotic time series for a Lorenz System . The procedure VecTS –of the TimeS package–reconstructs the state vectors. Our code still attaches all routines of this package in the present version [1,5]. This paper does not show the outputs of the computational procedures. We invite the reader to run the Maple worksheet from the additional file LinMapTS.mw . If the selection Analysis=1 is included in the commands, the routines perform analytical tasks. A graphic with the calculated and actual values of the observables provides a monitoring of the global fitting performed by the program LinGfiTS . In the command ConfiTS , the previous argument displays an analytical histogram beyond the results of Shapiro–Wilk test for the residuals ϵ j − τ in the global mapping. When the normality test gives a positive result, then the confidence level for a prediction time L τ is well-established [1]. For a prediction error ϵ ́ j + τ , this magnitude now is given by (6) L τ = 1 − 2 2 π ∫ ϵ ́ j + τ / σ τ ∞ exp { − ϵ j − τ 2 2 σ τ 2 } d ( ϵ j − τ / σ τ ) . (7) P pol ( | x 〉 ) = a 1 X 1 + a 2 X 2 + a 3 X 3 + a 4 X 1 2 + a 5 X 1 X 2 + a 6 X 1 X 3 + a 7 X 2 2 + a 8 X 2 X 3 + a 9 X 3 2 (8) P pot ( | x 〉 ) = a 1 X 1 | X 1 | | X 1 | 0.9 + a 2 X 1 | X 1 | | X 1 | 1.1 + a 3 X 2 | X 2 | | X 2 | 0.9 + a 4 X 2 | X 2 | | X 2 | 1.1 + a 5 X 3 | X 3 | | X 3 | 0.9 + a 6 X 3 | X 3 | | X 3 | 1.1 + a 7 X 1 X 2 X 3 | X 1 X 2 X 3 | | X 1 X 2 X 3 | 0.9 + a 8 X 1 X 2 X 3 | X 1 X 2 X 3 | | X 1 X 2 X 3 | 1.1 (9) P log ( | x r 〉 ) = a 1 X 1 + a 2 X 2 + a 3 X 3 + a 4 ln ( 1 + 1 10 cos ( X 1 ) ) + a 5 ln ( 1 + 1 10 cos ( X 2 ) ) + a 6 ln ( 1 + 1 10 cos ( X 3 ) ) + a 7 ln ( 1 + 1 10 sin ( X 1 ) ) + a 8 ln ( 1 + 1 10 sin ( X 2 ) ) + a 9 ln ( 1 + 1 10 sin ( X 3 ) ) . A global mapping for the entire time series presented in previous paragraphs has been explored in the first paper about the package LinMapTS [1]. We revisited this application with the inclusion of the optional argument Analysis=1 . The predictor (10) P ( | x r 〉 ) = a 1 X 1 + a 2 X 2 + a 3 X 3 + ⋯ + a 55 X 3 5 is automatically generated by the LinMapTS with the optional argument Degree=5 . The choice Level=39 assures that the interval for global mapping covers the whole series. The commands in this case are Fig. 1 shows the graphical analysis which the routines of the LinMapTS package make available from the commands above. The plotting in Fig. 1 (a) assures that the global fitting has high-quality. It is according to the accurate forecast presented in the previous version of the package [1]. The histogram–conjugated with the standard normal curve in Fig. 1 (b)–is compatible with the positive result of the Shapiro–Wilk test for the residuals’ distribution. With the non-polynomial predictors P pot (8) –assigned as F [ 2 ] –and P log (9) –assigned as F [ 3 ] –we got forecasts more accurate than the polynomial form P pol (7) –assigned as F [ 1 ] . The commands employed the optional argument Func . They have the structure shown below. Table 1 presents results of the two routines for the same time series studied in the previous example. The observable of interest in the prediction is V[723][1] . Although the alternative global mappings have required greater runtime than the polynomial function, the time measures show that the computational cost remains small. From the errors | ϵ ́ j + τ | the most accurate predictor is P log , followed by P pot . In this case study, the polynomial function P pol gives the larger error. So, the alternative functional forms had implied increasing the accuracy in forecasting. Another point in favor of the method is the comparison between the actual errors | ϵ ́ j + τ | and 3 σ τ . For the alternative predictors, the inequality | ϵ ́ j + τ | < 3 σ τ corresponds to the positive result in the normality test. However, the polynomial mapping did not permit to establish a credible confidence level L τ (6) in the prediction. From this application and many others to list, one concludes that this update represents a substantial improvement in the predictive capacity of the LinMapTS package. Acknowledgments L.G.S. Duarte and L.A.C.P. da Mota wish to thank Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ) for the Research Grant. [1] P. Alves, L. Duarte, L. da Mota, Computer Physics Communications 207 (2016) 325–340. [2] F. Takens, Detecting strange attractors in turbulence, in: D. Rand, L.-S. Young (Eds.), Dynamical Systems and Turbulence, Warwick 1980, Vol. 898 of Lecture Notes in Mathematics, Springer Berlin Heidelberg, Berlin, 1981, pp. 366–381. [3] M. Casdagli, Physica D: Nonlinear Phenomena 35 (3) (1989) 335–356. [4] D. Ruelle, Chaotic evolution and strange attractors: the statistical analysis of time series for deterministic nonlinear systems, Cambridge University Press, Cambridge New York, 1989. [5] P. Alves, L. Duarte, L. da Mota, Computer Physics Communications 207 (2016) 539–541. [ABSTRACT FROM AUTHOR] - Published
- 2017
- Full Text
- View/download PDF
21. Global Fitting Functions for Kinetics of Fe-Selective Chlorination in Ilmenite and Successive Chlorination of Beneficiated TiO2
- Author
-
Yong Sun Won, Duk-Yong Song, Yong Ha Kim, Eun-Jin Jung, and Dong-Kyu Chung
- Subjects
Materials science ,Chemical engineering ,Kinetics ,engineering ,Global fitting ,General Materials Science ,engineering.material ,Ilmenite - Published
- 2019
- Full Text
- View/download PDF
22. Accurate Measurement of Raman Depolarization Ratio in Gaseous CO2.
- Author
-
Yu-juan Jin, Yuan-qin Yu, Yu-xi Wang, Ke Lin, Xiao-guo Zhou, and Shi-lin Liu
- Abstract
The Raman depolarization ratios of gaseous CO
2 in the spectral range of 1240-1430 cm-1 are determined with a sensitive photoacoustic Raman spectroscopy, and more accurate data compared to the literature results are presented. The precision of the obtained depolarization ratio is achieved by measuring and fitting the dependence of the PARS signal intensity on the cross angle between the polarizations of two incident laser beams. [ABSTRACT FROM AUTHOR]- Published
- 2015
- Full Text
- View/download PDF
23. Improvements to facial contour detection by hierarchical fitting and regression.
- Author
-
Irie, Atsushi, Takagiwa, Mutsuki, Moriyama, Kozo, and Yamashita, Takayoshi
- Abstract
There are many methods based on shape and texture models for detecting eye and mouth contour points from facial images. They reduce the false positive rate by utilizing a global model and adapting it for a given face. Changes to facial expressions are coupled with changes to the shapes of eyes and mouth, and a global facial model in itself cannot be adapted to all human facial expressions. Therefore, a hierarchical model fitting approach has been developed, whereby the global fitting captures the facial shape using the global model and the local fitting captures the each facial parts using these local models. This can detect facial contours with high accuracy for expressions to which the global model cannot be adapted. [ABSTRACT FROM PUBLISHER]
- Published
- 2011
- Full Text
- View/download PDF
24. Improvement in global forecast for chaotic time series.
- Author
-
Alves, P.R.L., Duarte, L.G.S., and da Mota, L.A.C.P.
- Subjects
- *
CHAOS theory , *COMPUTATIONAL complexity , *TIME series analysis , *POLYNOMIALS , *ERROR analysis in mathematics , *ALGORITHMS - Abstract
In the Polynomial Global Approach to Time Series Analysis, the most costly (computationally speaking) step is the finding of the fitting polynomial. Here we present two routines that improve the forecasting. In the first, an algorithm that greatly improves this situation is introduced and implemented. The heart of this procedure is implemented on the specific routine which performs a mapping with great efficiency. In comparison with the similar procedure of the TimeS package developed by Carli et al. (2014), an enormous gain in efficiency and an increasing in accuracy are obtained. Another development in this work is the establishment of a level of confidence in global prediction with a statistical test for evaluating if the minimization performed is suitable or not. The other program presented in this article applies the Shapiro–Wilk test for checking the normality of the distribution of errors and calculates the expected deviation. The development is employed in observed and simulated time series to illustrate the performance obtained. Program summary Program title: LinMapTS Catalogue identifier: AFAJ_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AFAJ_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.: 9850 No. of bytes in distributed program, including test data, etc.: 223921 Distribution format: tar.gz Programming language: Maple 16. Computer: Any capable of running Maple. Operating system: Any capable of running Maple. Tested on Windows ME, Windows XP, Windows 7. RAM: 128 MB bytes Classification: 4.3, 4.9, 5. Nature of problem: Time series analysis and improving forecast capability. Solution method: The basis of the solution method is the result published in [1]. Restrictions: Global variables X [i] are used in the generated map; If more than 2000 vectors are employed in the global mapping the normality test is not applicable. Unusual features: When many polynomial coefficients are calculated (e.g., 55) their values can be different in distinct computers. These discrepancies do not affect significantly the accuracy in forecasting. Running time: A few seconds are required for the usual applications. References: [1] H. Carli, L. Duarte, L. da Mota, A maple package for improved global mapping forecast, Comp. Phts. Comm. 185(2014)1115 [ABSTRACT FROM AUTHOR]
- Published
- 2016
- Full Text
- View/download PDF
25. A COMPREHENSIVE STUDY OF THE PROTON STRUCTURE: FROM PDFS TO WIGNER FUNCTIONS
- Subjects
BC sum rules ,Global fitting ,Proton structure ,High Energy Physics::Lattice ,Physics ,High Energy Physics::Phenomenology ,Generalized TMDs ,Matching ,High Energy Physics::Experiment ,Quasi-distributions ,Nuclear Experiment - Abstract
It has been known since the 1930’s that protons and neutrons, collectively called as nucleons, are not “point-like” elementary particles, but rather have a substructure. Today, we know from Quantum Chromodynamics (QCD) that nucleons are made from quarks and gluons, with gluons being the elementary force carriers for strong interactions. Quarks and gluons are collectively called as partons. The substructure of the nucleons can be described in terms of parton correlation functions such as Form Factors, (1D) Parton Distribution Functions (PDFs) and their 3D generalizations in terms of Transverse Momentum-dependent parton Distributions (TMDs) and Generalized Parton Distributions (GPDs). All these functions can be derived from the even more general Generalized Transverse Momentum-dependent Distributions (GTMDs). This dissertation promises to provide an insight into all these functions from the point of view of their accessibility in experiments, from model calculations, and from their direct calculation within lattice formulations of QCD. In the first part of this dissertation, we identify physical processes to access GTMDs. By considering the exclusive double Drell-Yan process, we demonstrate, for the very first time, that quark GTMDs can be measured. We also show that exclusive double-quarkonium production in nucleon-nucleon collisions is a direct probe of gluon GTMDs. In the second part of this dissertation, we shift our focus to the “parton quasi-distributions”. Over the last few decades, lattice QCD extraction of the full x-dependence of the parton distributions has always been prohibited by the explicit time-dependence of the correlation functions. In 2013, there was a path-breaking proposal by X. Ji to calculate instead parton quasi-distributions (quasi-PDFs). The procedure of “matching” is a crucial ingredient in the lattice QCD extraction of parton distributions from the quasi-PDF approach. We address the matching for the twist-3 PDFs gT (x), e(x), and hL(x) for the very first time. We pay special attention to the challenges involved in the calculations due to the presence of singular zero-mode contributions. We also present the first-ever lattice QCD results for gT (x) and hL(x) and we discuss the impact of these results on the phenomenology. Next, we explore the general features of quasi-GPDs and quasi-PDFs in diquark spectator models. Furthermore, we address the Burkhardt-Cottingham-type sum rules for the relevant light-cone PDFs and quasi-PDFs in a model-independent manner and also check them explicitly in perturbative model calculations. The last part of this dissertation focuses on the extraction g1T (x,~k2⊥) TMD for the very first time from experimental data using Monte Carlo techniques. This dissertation therefore unravels different aspects of the distribution functions from varied perspectives.
- Published
- 2021
- Full Text
- View/download PDF
26. Multispectrum Rotational States Distribution Thermometry
- Author
-
Thomas Puppe, Wilhelm Kaenders, Rafal Wilk, Riccardo Gotti, Yuriy Mayzlin, Marco Lamperti, Paolo Laporta, Davide Gatti, Julian Robinson-Tait, Marco Marangoni, Szymon Wójtewicz, Patrick Leisching, Felix Rohde, and Bidoor Alsaif
- Subjects
Physics ,Condensed Matter::Quantum Gases ,Materials science ,Distribution (number theory) ,Spectrometer ,Absorption spectroscopy ,Physics::Instrumentation and Detectors ,Optical testing ,Global fitting ,Physics::Optics ,Cavity ring-down spectroscopy ,Computational physics ,Thermometer ,Atomic physics ,Nuclear Experiment ,Absorption (electromagnetic radiation) ,Phase modulation ,Molecular absorption - Abstract
We exploit a widely tunable comb-locked frequency-swept synthesizer to test a new optical approach to primary gas thermometry based on a global fitting of multiple molecular absorption lines of the same band at different pressures.
- Published
- 2021
27. Defining a kinetic mechanism for l-DOPA 2,3 dioxygenase, a single-domain type I extradiol dioxygenase from Streptomyces lincolnensis.
- Author
-
Colabroy, Keri L., Smith, Ian R., Vlahos, Alexander H.S., Markham, Androo J., and Jakubik, Matthew E.
- Subjects
- *
ANALYTICAL mechanics , *DOPA , *DIOXYGENASES , *STREPTOMYCES , *LINCOMYCIN , *ANTIBIOTICS - Abstract
Abstract: l-DOPA-2,3-dioxygenase from Streptomyces lincolnensis is a single domain type I extradiol dioxygenase of the vicinal oxygen chelate superfamily and catalyzes the second step in the metabolism of the propylhygric acid moiety of the antibiotic, lincomycin. In this report, the kinetic mechanism of l-DOPA dioxygenase is interrogated using stopped-flow in order to determine microscopic rate constants. Pre-steady state, progress curve and steady-state data were combined in a global kinetic analysis using KinTek Explorer in order to define and constrain a kinetic model for the type I l-DOPA dioxygenase. The data are best described by a four step mechanism, in which the cyclization of the enzymatic product is not enzyme catalyzed. [Copyright &y& Elsevier]
- Published
- 2014
- Full Text
- View/download PDF
28. Global fitting and parameter identifiability for amyloid-β aggregation with competing pathways
- Author
-
Preetam Ghosh, Vijay Rangachari, Ashwin Vaidya, Pratip Rana, and Priyankar Bose
- Subjects
0301 basic medicine ,Chemistry ,Systems biology ,Kinetics ,Global fitting ,Protein aggregation ,03 medical and health sciences ,030104 developmental biology ,0302 clinical medicine ,Ordinary differential equation ,Identifiability ,Uncertainty quantification ,Biological system ,030217 neurology & neurosurgery ,Free parameter - Abstract
Aggregation of the amyloid-$\beta$(A$\beta$) protein has been implicated in Alzheimer’s disease (AD). Since, low molecular weight A$\beta$ aggregates are hypothesized to serve as the primary toxic species in AD pathogenesis, significant research has been conducted to understand the mechanistic details of the aggregation process. We previously demonstrated that heterotypic interactions between A$\beta$ and fatty acids (FAs) can lead to competing pathways of A$\beta$ aggregation, termed as the off-pathway; this off-pathway kinetics can also be modulated by FA concentrations as captured by mass action models. We employed ensemble kinetics simulations which uses a system of Ordinary Differential Equations to model the competing on-and off-pathways of $A\beta$ aggregation that were trained and validated by biophysical experiments. However, these models had several rate constants, treated as free parameters to be estimated, which resulted in over-fitting of the model. Hence, in this paper, we present a global fitting based method to accurately identify the rate constants involved in the complex competing pathway model of $A\beta$ aggregation. We additionally employ detailed parameter identifiability tests for uncertainty quantification using the profile likelihood method. Since, the emergence of off-or on-pathway aggregates are typically controlled by a narrow set of rate constants, it is imperative to rigorously identify the proper rate constants involved in these pathways. These rate constants serve as a basis for future experiments on modulating the aggregation pathways to populate a particular possibly less toxic oligomeric species. The obtained rate constants also motivate new biophysical experiments to better understand the mechanisms of amyloid aggregation in other neurodegenerative diseases.
- Published
- 2020
- Full Text
- View/download PDF
29. Analysis of nonideality: insights from high concentration simulations of sedimentation velocity data
- Author
-
Walter F. Stafford, R. T. Wright, Peter J. Sherwood, and John J. Correia
- Subjects
0301 basic medicine ,030103 biophysics ,Concentration dependence ,Diffusion ,Biophysics ,Thermodynamics ,03 medical and health sciences ,Thermodynamic nonideality ,Physics ,High concentration ,Sedimentation velocity ,Global fitting ,Antibodies, Monoclonal ,General Medicine ,Sedimentation ,Models, Theoretical ,Software package ,2nd Virial coefficient ,Kinetics ,030104 developmental biology ,Spectrometry, Fluorescence ,Weak association ,Hydrodynamic nonideality ,Original Article ,Ultracentrifugation - Abstract
The Aviv fluorescence detection system (Aviv-FDS) has allowed the performance of sedimentation velocity experiments on therapeutic antibodies in highly concentrated environments like formulation buffers and serum. Methods were implemented in the software package SEDANAL for the analysis of nonideal, weakly associating AUC data acquired on therapeutic antibodies and proteins (Wright et al. Eur Biophys J 47:709–722, 2018, Anal Biochem 550:72–83, 2018). This involved fitting both hydrodynamic, ks, and thermodynamic, BM1, nonideality where concentration dependence is expressed as s = so/(1 + ksc) and D = Do(1 + 2BM1c)/(1 + ksc) and so and Do are values extrapolated to c = 0 (mg/ml). To gain insight into the consequences of these phenomenological parameters, we performed simulations with SEDANAL of a monoclonal antibody as a function of ks (0–100 ml/g) and BM1 (0–100 ml/g). This provides a visual understanding of the separate and joint impact of ks and BM1 on the shape of high-concentration sedimentation velocity boundaries and the challenge of their unique determination by finite element methods. In addition, mAbs undergo weak self- and hetero-association (Yang et al. Prot Sci 27:1334–1348, 2018) and thus we have simulated examples of nonideal weak association over a wide range of concentrations (1–120 mg/ml). Here we demonstrate these data are best analyzed by direct boundary global fitting to models that account for ks, BM1 and weak association. Because a typical clinical dose of mAb is 50–200 mg/ml, these results have relevance for biophysical understanding of concentrated therapeutic proteins.
- Published
- 2020
30. Global Small-Angle X-ray Scattering Data Analysis of Triacylglycerols in the Molten State (Part I)
- Author
-
Josélio Batista Vieira, Michael Rappolt, Megan Povey, Amin Sadeghpour, and Marjorie Ladd Parada
- Subjects
Materials science ,010405 organic chemistry ,Scattering ,Small-angle X-ray scattering ,Isotropy ,Global fitting ,010402 general chemistry ,01 natural sciences ,Uncorrelated ,0104 chemical sciences ,Surfaces, Coatings and Films ,Computational physics ,law.invention ,Molten state ,law ,Materials Chemistry ,Lamellar structure ,Physical and Theoretical Chemistry ,Crystallization - Abstract
The study of triacylglycerols (TAGs) in their molten state is of fundamental importance for a deeper understanding of the TAG crystallization processes, being highly relevant for both manufacturing and medical applications. Although different models have been proposed to explain the nanostructured nature of the fluid state of TAGs, none of them are fully satisfactory. In this paper, we propose a new model consisting of positionally uncorrelated lamellar TAG assemblies embedded in an isotropic medium that assist as prenucleating structures. This model was validated by applying a novel global fitting method, resulting in an excellent agreement with the small-angle X-ray scattering data. A deeper analysis of the scattering patterns at different temperatures, both in cooling and heating directions, allowed us further to detect the crystalline traces of TAGs even after heating to 40 °C and record, on cooling, the onset of crystallization at 30-25 °C. The application of the presented novel model not only explains the outstandingly structured fluid of molten TAGs, but also lays the basis for analyzing first the crystallization steps in greater detail, which is outlined in our follow-up paper "Global Small-Angle X-ray Scattering Data Analysis of Triacylglycerols in the α-Phase (Part II)".
- Published
- 2018
- Full Text
- View/download PDF
31. River extraction of SAR images via active contours driven by adaptive global fitting energies
- Author
-
Bin Han and Yiquan Wu
- Subjects
Synthetic aperture radar ,Active contour model ,ComputingMethodologies_IMAGEPROCESSINGANDCOMPUTERVISION ,Global fitting ,020206 networking & telecommunications ,02 engineering and technology ,Computer Science::Graphics ,Remote sensing (archaeology) ,0202 electrical engineering, electronic engineering, information engineering ,General Earth and Planetary Sciences ,020201 artificial intelligence & image processing ,Extraction (military) ,Satellite imagery ,Physics::Atmospheric and Oceanic Physics ,Geology ,ComputingMethodologies_COMPUTERGRAPHICS ,Remote sensing - Abstract
This article proposed a novel active contour model based on adaptive global fitting energies to achieve accurate river extraction of synthetic aperture radar (SAR) images. First, inspired by the me...
- Published
- 2018
- Full Text
- View/download PDF
32. On the kinetic mechanism of dimethylarginine dimethylaminohydrolase.
- Author
-
Johnson, Corey M. and Fast, Walter
- Subjects
- *
NITRIC-oxide synthases , *HYDROLYSIS , *STABILITY constants , *ASYMMETRIC dimethylarginine , *PSEUDOMONAS aeruginosa , *NITRIC oxide - Abstract
[Display omitted] Dimethylarginine dimethylaminohydrolase (DDAH, EC 3.5.3.18) catalyzes the hydrolysis of asymmetric N ω, N ω-dimethyl- l -arginine (ADMA), an endogenous inhibitor of human nitric oxide synthases. The active-site cysteine residue has been proposed to serve as the catalytic nucleophile, forming an S -alkylthiourea reaction intermediate, and serving as a target for covalent inhibitors. Inhibition can lead to ADMA accumulation and downstream inhibition of nitric oxide production. Prior studies have provided experimental evidence for formation of this covalent adduct but have not characterized it kinetically. Here, rapid quench-flow is used with ADMA and the DDAH from Pseudomonas aeruginosa to determine the rate constants for formation (k 2 = 17 ± 2 s−1) and decay (k 3 = 1.5 ± 0.1 s−1) of the covalent S -alkylthiourea adduct. A minimal kinetic mechanism for DDAH is proposed that supports the kinetic competence of this species as a covalent reaction intermediate and assigns the rate-limiting step in substrate turnover as hydrolysis of this intermediate. This work helps elucidate the different reactivities of S -alkylthiourea intermediates found among the mechanistically diverse pentein superfamily of guanidine-modifying enzymes and provides information useful for inhibitor development. [ABSTRACT FROM AUTHOR]
- Published
- 2022
- Full Text
- View/download PDF
33. Fluorescent-labeled antibodies: Balancing functionality and degree of labeling
- Author
-
Vira, Shaleen, Mekhedov, Elena, Humphrey, Glen, and Blank, Paul S.
- Subjects
- *
FLUORESCENT antibody technique , *AFFINITY labeling , *BIOCONJUGATES , *MONOCLONAL antibodies , *HEMAGGLUTININ , *ENZYME-linked immunosorbent assay , *ANALYTICAL biochemistry - Abstract
Abstract: A critical assumption in using labeled antibodies is that the conjugation reaction has no deleterious effects on antibody avidity. This study demonstrates that this assumption need not hold true and presents a methodology to quantitatively determine the degree of inactivation and/or changes in antibody–antigen binding that can occur with conjugation. Fluorescein isothiocyanate (FITC) was conjugated to a mouse monoclonal antibody, Fc125, against hemagglutinin (HA) using varying fluorophore/protein (F:P) labeling ratios. Antibody binding, as a function of the F:P labeling ratio, was evaluated using a kinetic enzyme-linked immunosorbent assay (ELISA) and analyzed using global fitting. A two-parameter adjustment of the antibody concentration and the maximum rate was sufficient to describe the rate changes. The concentration parameter dominated the rate changes, consistent with the hypothesis that the coupling reaction inactivated an increasing fraction of the antibody population with a smaller change (∼15% at the highest F:P ratio) in antibody–antigen binding. An optimal F:P ratio that minimized both inactivation and unlabeled antibody was calculated. This procedure can be used to prepare functional, labeled antibody reagents with defined activity and can aid in quantitative applications where the stoichiometry and functionality of the labeled antibody are critical. [Copyright &y& Elsevier]
- Published
- 2010
- Full Text
- View/download PDF
34. Global Kinetic Explorer: A new computer program for dynamic simulation and fitting of kinetic data
- Author
-
Johnson, Kenneth A., Simpson, Zachary B., and Blom, Thomas
- Subjects
- *
COMPUTER software , *ENZYME kinetics , *COMPUTER simulation , *BIOCHEMISTRY , *NUMERICAL integration , *REGRESSION analysis software - Abstract
Abstract: We describe a new dynamic kinetic simulation program that allows multiple data sets to be fit simultaneously to a single model based on numerical integration of the rate equations describing the reaction mechanism. Unlike other programs that allow fitting based on numerical integration of rate equations, in the dynamic simulation rate constants, output factors, and starting concentrations of reactants can be scrolled while observing the change in the shape of the simulated reaction curves. Fast dynamic simulation facilitates the exploration of initial parameters that serve as the starting point for nonlinear regression in fitting data and facilitates exploration of the relationships between individual constants and observable reactions. The exploration of parameter space by dynamic simulation provides a powerful tool for learning kinetics and for evaluating the extent to which parameters are constrained by the data. This feature is critical to avoid overly complex models that are not supported by the data. [Copyright &y& Elsevier]
- Published
- 2009
- Full Text
- View/download PDF
35. Ultrafast isomerization dynamics of retinal in bacteriorhodopsin as revealed by femtosecond absorption spectroscopy.
- Author
-
Wu YiShi, Zhong Sheng, Al XiCheng, Hu KunSheng, and Zhang JianPing
- Subjects
- *
PHOTOISOMERIZATION , *BACTERIORHODOPSIN , *ENERGY metabolism , *ISOMERISM , *ISOMERIZATION - Abstract
A femtosecond (fs) broad-band absorption apparatus was used to measure the early photoisomerization process of bacteriorhodopsin's (BR) photocycle to reveal the character of the important intermediate of J625 and to obtain a deeper understanding of the role of photoisomerization in BR photocycle. Two time constants of 0.5 PS (95%) and 2.0 ps (5%) were brought out by global fitting on thirty curves in the near-infrared region. We suggest that the first time constant results from the decay of I460 intermediate, and the longer component might be associated with BR isomer. The global analysis over 450, 540, 630,710 and 870 nm traces identified two time constants, ∼0.5 and ∼3 ps. The slower component can be extracted from the processes of both J625→BR568 (540 nm) and J625→K590 (630 nm), suggesting J-intermediate takes a partial cis configuration. The obvious negative feature in early delay time of 700-780 nm regions was attributed to the radiative transition (stimulated emission) from the Franck-Condon active configuration along the isomerization potential surface of all-trans-retinal. [ABSTRACT FROM AUTHOR]
- Published
- 2008
- Full Text
- View/download PDF
36. Neutron Diffraction Cryoporometry—A measurement technique for studying mesoporous materials and the phases of contained liquids and their crystalline forms
- Author
-
Webber, J. Beau W. and Dore, John C.
- Subjects
- *
NEUTRON diffraction , *NUCLEAR magnetic resonance spectroscopy , *CRYSTALLIZATION , *NANOTECHNOLOGY - Abstract
Abstract: Neutron diffraction is a standard method for determining the structure of matter on an atomic scale; NMR cryoporometry is a recent widely applicable technique for characterising structure on a 2nm to scale. An extension of these techniques is described, Neutron Diffraction Cryoporometry (NDC). The information from a set of neutron diffraction measurements of liquids and their crystalline forms in meso-pores, as a function of temperature, is displayed as a cryoporometry graph. The data may then be conveniently interpreted using the Gibbs–Thomson relationship by analogy with the existing technique, NMR cryoporometry. Clear information is thus obtained on the relationship between phase and nano-structure, in a form well suited to further analysis. This method is applied to an equilibrium study of water/ice in SBA-15 templated silicas, as model nano- to meso-structured systems. The method described here uses global pattern matching (a one-dimensional morphing algorithm inside a linear least-squares fitting algorithm) applied to the full range of the diffraction data. This is a rapid method by comparison with the conventional method of fitting individual (overlapping) peaks, and has already led to NMR observations indicating plastic (rotator phase) ice in the same system. [Copyright &y& Elsevier]
- Published
- 2008
- Full Text
- View/download PDF
37. Torsion–rotation global analysis and database for the CH318OH isotopomer of methanol
- Author
-
Fisher, J., Paciga, G., Xu, Li-Hong, Zhao, S.B., Moruzzi, G., and Lees, R.M.
- Subjects
- *
INFRARED spectra , *SPECTRUM analysis , *DEFORMATIONS (Mechanics) , *PROPERTIES of matter - Abstract
Abstract: Fourier-transform far-infrared spectra of CH3 18OH in the 15–470cm−1 region have been analyzed by means of the Ritz assignment program. The far-infrared data have been combined with the literature microwave and millimeter-wave measurements in a full global fitting of the first three torsional states (ν t =0, 1, and 2) of the CH3 18OH ground vibrational state. The fitted dataset includes 550 microwave and millimeter-wave lines and more than 17000 Fourier-transform transitions covering the quantum number ranges J ⩽30, K ⩽15, and ν t ⩽2. With incorporation of 79 adjustable parameters, the global fit achieved convergence with an overall weighted standard deviation of 1.072, essentially to within the assigned measurement uncertainties of ±50kHz for almost all of the microwave and millimeter-wave lines and ±6MHz (0.0002cm−1) to ±15MHz (0.0005cm−1) for the Fourier-transform far-infrared measurements. Based on the global fit results, a database has been compiled containing transition frequencies, quantum numbers, lower state energies and transition strengths. This database will provide support for present and future astronomical studies, such as the on-going Orion surveys in preparation for the launch of the Herschel Space Observatory, in identifying isotopic methanol contributions to interstellar spectra. [Copyright &y& Elsevier]
- Published
- 2007
- Full Text
- View/download PDF
38. Gating of the HypoPP-1 mutations: II. Effects of a calcium-channel agonist BayK 8644.
- Author
-
Kuzmenkin, Alexey, Chao Hang, Kuzmenkina, Elza, and Jurkat-Rott, Karin
- Subjects
- *
HYPOKALEMIA , *CALCIUM channels , *RABBITS , *ANIMAL models in research , *MUTAGENESIS , *IONIC structure , *DIAGNOSIS - Abstract
L-type calcium-channel mutations causing hypokalemic periodic paralysis type 1 (HypoPP-1) have pronounced “loss-of-function” features and stabilize the less-selective second open state O2, as we demonstrated in the companion paper. Here, we compared the effects of the L-type calcium-channel activator (±)BayK 8644 (BayK) on the heterologously expressed wild-type (WT) calcium channel, rabbit Cav1.2 HypoPP-1 analogs, and two double mutants (R650H/R1362H, R650H/R1362G). Our goal was to elucidate (1) whether the “loss-of-function” in HypoPP-1 can be compensated by BayK application, (2) how the less-selective open state is affected by BayK in WT and HypoPP-1 mutants, as well as (3) to gain an insight into BayK mechanism of action. Ionic currents were examined by whole-cell patch-clamp and analyzed by the global-fitting procedure. Our results imply that (1) BayK promotes channel activation, but equalized the differences among the WT and mutants, thus attenuating HypoPP-related effects on activation and deactivation; (2) BayK binds to the first open state O1, and then serves as a catalyst for O2 formation; (3) binding of BayK is impaired in the HypoPP mutants, thus affecting the formation of the less-selective second open state; (4) BayK affects cooperativity between the single HypoPP-1 mutations at all stages of the channel gating; and (5) BayK favoring of O2 lowers calcium-channel selectivity. [ABSTRACT FROM AUTHOR]
- Published
- 2007
- Full Text
- View/download PDF
39. Gating of the HypoPP-1 mutations: I. Mutant-specific effects and cooperativity.
- Author
-
Kuzmenkin, Alexey, Chao Hang, Kuzmenkina, Elza, and Jurkat-Rott, Karin
- Subjects
- *
HYPOKALEMIA , *PARALYSIS , *ACTIVE biological transport , *CALCIUM channels , *GENETIC mutation , *GENETICS - Abstract
Hypokalemic periodic paralysis type 1 (HypoPP-1) is a hereditary muscular disorder caused by point mutations in the gene encoding the voltage-gated Ca2+ channel α subunit (Cav1.1). Despite extensive research, the results on HypoPP-1 mutations are minor and controversial, as it is difficult to analyse Ca2+ channel activation macroscopically due to an existence of two open states. In this study, we heterologously expressed the wild-type and HypoPP-1 mutations introduced into the rabbit cardiac Ca2+ channel (R650H, R1362H, R1362G) in HEK-293 cells. To examine the cooperative effects of the mutations on channel gating, we expressed two double mutants (R650H/R1362H, R650H/R1362G). We performed whole-cell patch-clamp and, to obtain more information, applied a global fitting procedure whereby several current traces elicited by different potentials were simultaneously fit to the kinetic model containing four closed, two open and two inactivated states. We found that all HypoPP-1 mutations have “loss-of-function” features: D4/S4 mutations shift the equilibrium to the closed states, which results in reduced open probability, shorter openings and, therefore, in smaller currents, and the D2/S4 mutant slows the activation. In addition, HypoPP-1 histidine mutants favored the second open state O2 with a possibly lower channel selectivity. Cooperativity between the D2/S4 and D4/S4 HypoPP-1 mutations manifested in dominant effects of the D4/S4 mutations on kinetics of the double mutants, suggesting different roles of D2/S4 and D4/S4 voltage sensors in the gating of voltage-gated calcium channels. [ABSTRACT FROM AUTHOR]
- Published
- 2007
- Full Text
- View/download PDF
40. Global effective Hamiltonians of 16O13C17O and 16O13C18O improved from CW-CRDS observations in the 5900–7000cm−1 region
- Author
-
Perevalov, B.V., Kassi, S., Romanini, D., Perevalov, V.I., Tashkun, S.A., and Campargue, A.
- Subjects
- *
HAMILTONIAN systems , *POLYADIC algebras , *CARBON dioxide , *SPECTRUM analysis , *MOLECULAR spectroscopy - Abstract
Abstract: The parameters of the polyad models of the effective Hamiltonian of the 16O13C17O and 16O13C18O isotopologues of carbon dioxide have been refined by the least-squares fittings to the line positions collected from the literature. Such refinement has become necessary as the observed dataset has been significantly extended by our CW-CRDS observations in the 5900–7000cm−1 region. In the case of the 16O13C17O isotopologue, 1151 line positions of 11 bands have been used to refine the effective Hamiltonian parameters published by Chédin [A. Chédin, J. Mol. Spectrosc. 76 (1979) 430–491]. With the obtained set of parameters, the collected line positions are reproduced with a RMS (root mean squares of the residuals) equal to 0.0013cm−1. In the case of the 16O13C18O isotopologue, 61 parameters of the effective Hamiltonian were fitted to more than 6410 line positions. A weighted standard deviation of χ =1.77 and a global RMS of 0.0017cm−1, close to the experimental accuracy, were achieved. However, several rotational levels of the 31113 state (P =10) could not be reproduced in the frame of this polyad model and were then excluded from the fit. We found that these levels are affected by an anharmonic resonance interaction with the 51106 vibrational state (P =11) leading to energy shifts up to 0.060cm−1 and significant intensity transfer to several extra lines which could be detected. The coupling matrix element has been estimated to 0.11cm−1 from the detailed analysis of the experimental spectrum. This is the first evidence of an interpolyad resonance interaction in the case of the carbon dioxide molecule. In order to extend the input spectroscopic information, the weak lines left unassigned in our previous analysis of the CW-CRDS spectrum of the 13C enriched carbon dioxide [Y. Ding, P. Macko, D. Romanini, V.I. Perevalov, S.A. Tashkun, J.-L. Teffo, S.-M. Hu, A. Campargue, J. Mol. Spectrosc. 226 (2004) 146–160.] have been revisited. Thirteen 13C16O2 bands, one 16O13C17O band and two 16O13C18O bands could be newly assigned together with a number of transitions corresponding to high J values of previously observed bands. The spectroscopic constants G v, B v, and D v for the unperturbed bands have been fitted to the observed line positions. [Copyright &y& Elsevier]
- Published
- 2007
- Full Text
- View/download PDF
41. Alzheimer’s-disease-associated conformation of intrinsically disordered tau protein studied by intrinsically disordered protein liquid-phase competitive enzyme-linked immunosorbent assay
- Author
-
Skrabana, Rostislav, Skrabanova-Khuebachova, Michaela, Kontsek, Peter, and Novak, Michal
- Subjects
- *
ALZHEIMER'S disease , *PROTEINS , *GLYCOPROTEINS , *BINDING energy - Abstract
Abstract: Tau protein, the major constituent of paired helical filaments in Alzheimer’s disease, belongs to the intrinsically disordered proteins (IDPs). IDPs are an emerging group in the protein kingdom characterized by the absence of a rigid three-dimensional structure. Disordered proteins usually acquire a “functional fold” upon binding to their interaction partner(s). This property of IDPs implies the need for innovative approaches to measure their binding affinity. We have mapped and measured the Alzheimer’s-disease-associated epitope on intrinsically disordered tau protein with a novel two-step sandwich competitive enzyme-linked immunosorbent assay (ELISA). This approach allowed us to determine the binding affinity of disordered tau protein in liquid phase without any disturbance to the competitive equilibrium and without any need for covalent or noncovalent modification of tau protein. Furthermore, the global fitting method, used for the reconstruction of tau binding curves, significantly improved the assay readout. The proposed novel competitive ELISA allowed us to determine the changes in the standard Gibbs energy of binding, thus enabling measurement of tau protein conformation in the core of paired helical filaments. IDP competitive ELISA results showed, for the first time, that the tau protein C terminus of the Alzheimer’s-disease-derived paired helical filaments core subunit adopts β-turn type I′ fold and is accessible from solution. [Copyright &y& Elsevier]
- Published
- 2006
- Full Text
- View/download PDF
42. Electroweak and QCD corrections to Z and W pole observables in the standard model EFT
- Author
-
Pier Paolo Giardino and Sally Dawson
- Subjects
Physics ,Quantum chromodynamics ,Particle physics ,Electroweak interaction ,Global fitting ,Effective field theory ,Order (ring theory) ,Observable ,Standard Model - Abstract
We compute the next-to-leading order QCD and electroweak corrections to $Z$ and $W$ pole observables using the dimension-6 Standard Model effective field theory (SMEFT) and present numerical results that can easily be included in global fitting programs. Limits on SMEFT coefficient functions are presented at leading order and at next-to-leading order under several assumptions.
- Published
- 2020
- Full Text
- View/download PDF
43. Visualization of transient absorption dynamics – towards a qualitative view of complex reaction kinetics
- Author
-
Satzger, H. and Zinth, W.
- Subjects
- *
PHOTOCHEMISTRY , *ABSORPTION , *NUMERICAL analysis , *SPECTRUM analysis - Abstract
Photochemical reactions may involve complex reaction schemes and their investigation requires spectroscopic techniques extending over a broad spectral range and over several orders of magnitude in time. Two commonly used numerical procedures to evaluate such data sets – global fitting and singular value decomposition – are discussed and a method for the qualitative visualization is proposed: differentiation of the transient spectra on a logarithmic scale allows to extract special kinetic components and to obtain reasonable starting information for subsequent fitting procedures. The proposed method should be well adapted to situations where e.g. due to unstable samples the data cannot be recorded at the precision required for unambiguous analysis by standard data handling procedures. The method is applied to synthetic data sets as well as to experimental data taken from femtosecond absorption experiments on the laser dye DCM in DMSO. [Copyright &y& Elsevier]
- Published
- 2003
- Full Text
- View/download PDF
44. Object-level structured contour map extraction.
- Author
-
Bergevin, R. and Bubel, A.
- Subjects
MAPS ,CONTOURS (Cartography) ,TOPOLOGY - Abstract
A new approach is proposed to extract an object-level structured contour map from the junctions in a 2D image. Local junction structures are paired in order to initialize and globally constrain the contour extraction process. An original optimization method is applied to detect and describe the contours in respect of their inherent shape and structure. Any given contour is described using a proper number of constant-curvature primitives according to the complexity of its shape. The contour map produced is represented as a graph. The characterized contours correspond to the edges of the graph and their topological structure is described by the junctions at the vertices. An elaborated experimental evaluation illustrates the capabilities of the proposed approach in demanding situations. [Copyright &y& Elsevier]
- Published
- 2003
- Full Text
- View/download PDF
45. A Maple package for improved global mapping forecast.
- Author
-
Carli, H., Duarte, L.G.S., and da Mota, L.A.C.P.
- Subjects
- *
COMPUTER software , *MATHEMATICAL mappings , *TIME series analysis , *COMPUTATIONAL complexity , *DATA analysis , *ALGORITHMS - Abstract
We present a Maple implementation of the well known global approach to time series analysis and some further developments designed to improve the computational efficiency of the forecasting capabilities of the approach. This global approach can be summarized as being a reconstruction of the phase space, based on a time ordered series of data obtained from the system. After that, using the reconstructed vectors, a portion of this space is used to produce a mapping, a polynomial fitting, through a minimization procedure, that represents the system and can be employed to forecast further entries for the series. In the present implementation, we introduce a set of commands, tools, in order to perform all these tasks. For example, the command VecTS deals mainly with the reconstruction of the vector in the phase space. The command GfiTS deals with producing the minimization and the fitting. ForecasTS uses all these and produces the prediction of the next entries. For the non-standard algorithms, we here present two commands: IforecasTS and NiforecasTS that, respectively deal with the one-step and the -step forecasting. Finally, we introduce two further tools to aid the forecasting. The commands GfiTS and AnalysTS, basically, perform an analysis of the behavior of each portion of a series regarding the settings used on the commands just mentioned above. Program summary: Program title: TimeS Catalogue identifier: AERW_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AERW_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.: 3001 No. of bytes in distributed program, including test data, etc.: 95018 Distribution format: tar.gz Programming language: Maple 14. Computer: Any capable of running Maple Operating system: Any capable of running Maple. Tested on Windows ME, Windows XP, Windows 7. RAM: 128 MB Classification: 4.3, 4.9, 5 Nature of problem: Time series analysis and improving forecast capability. Solution method: The method of solution is partially based on a result published in [1]. Restrictions: If the time series that is being analyzed presents a great amount of noise or if the dynamical system behind the time series is of high dimensionality (Dim≫3), then the method may not work well. Unusual features: Our implementation can, in the cases where the dynamics behind the time series is given by a system of low dimensionality, greatly improve the forecast. Running time: This depends strongly on the command that is being used. References: [1] Barbosa, L.M.C.R., Duarte, L.G.S., Linhares, C.A. and da Mota, L.A.C.P., Improving the global fitting method on nonlinear time series analysis, Phys. Rev. E 74, 026702 (2006). [ABSTRACT FROM AUTHOR]
- Published
- 2014
- Full Text
- View/download PDF
46. Global Fitting of the Response Surface via Estimating Multiple Contours of a Simulator
- Author
-
Feng Yang, C. Devon Lin, and Pritam Ranjan
- Subjects
Statistics and Probability ,Surface (mathematics) ,FOS: Computer and information sciences ,Physical system ,Machine Learning (stat.ML) ,Statistics - Applications ,Statistics - Computation ,01 natural sciences ,Methodology (stat.ME) ,010104 statistics & probability ,Statistics - Machine Learning ,0502 economics and business ,Applications (stat.AP) ,0101 mathematics ,Aerospace ,Statistics - Methodology ,Simulation ,Selection (genetic algorithm) ,Computation (stat.CO) ,Mathematics ,business.industry ,05 social sciences ,Process (computing) ,Global fitting ,13. Climate action ,050211 marketing ,business - Abstract
Computer simulators are nowadays widely used to understand complex physical systems in many areas such as aerospace, renewable energy, climate modeling, and manufacturing. One fundamental issue in the study of computer simulators is known as experimental design, that is, how to select the input settings where the computer simulator is run and the corresponding response is collected. Extra care should be taken in the selection process because computer simulators can be computationally expensive to run. The selection shall acknowledge and achieve the goal of the analysis. This article focuses on the goal of producing more accurate prediction which is important for risk assessment and decision making. We propose two new methods of design approaches that sequentially select input settings to achieve this goal. The approaches make novel applications of simultaneous and sequential contour estimations. Numerical examples are employed to demonstrate the effectiveness of the proposed approaches., Comment: 24 pages
- Published
- 2019
- Full Text
- View/download PDF
47. A new method for improved global mapping forecast.
- Author
-
Alves, P.R.L., Duarte, L.G.S., and da Mota, L.A.C.P.
- Subjects
- *
MATHEMATICAL mappings , *TIME series analysis , *PHASE space , *COMPUTATIONAL complexity , *PROGRAMMING languages , *DIFFERENTIAL equations - Abstract
The Maple package TimeS for time series analysis has a new feature and an improvement in forecasting by phase space reconstruction. An optional argument in the computational routines that allows the researcher to choose the different number of steps ahead to forecast. This update extends the running of the package with this new feature for the current versions of the Maple software too. New version program summary Program Title: TimeS Program Files doi: 10.17632/nhtmjc8yp8.1 Licensing provisions: GNU General Public License 3 Programming language: Maple17 Journal reference of previous version: Comput. Phys. Comm. 185 (2014) 1115 Does the new version supersede the previous version?: Yes. Nature of problem: Time series analysis and improving forecast capability. Solution method: The method of solution is published in [1]. Reasons for the new version: For a phenomenon described by an unknown low-dimensional dynamical system defined by a set of coupled differential equations x ̇ i = f i ( x ) , i = 1 , … , n , one generates a map M for which the time variable is increased by δ t : x i ( P + 1 ) = F i ( x ( P ) , δ t ) . Though another map M ¯ that approaches M , X ¯ i ( P + 1 ) is close to x i ( P + 1 ) when δ t → 0 [1]. (1) X ¯ i ( P + 1 ) = F ¯ i ( x ( P ) , δ t ) = ∑ k = 0 N ( δ t ) k k ! X k [ x i ( P ) ] . In the reconstruction’s scheme by delays method [2], x ( P ) is state vector reconstructed. Its components have regular spacings in the time series. The approach X ¯ i ( P + 1 ) corresponds to the prediction for the value x i ( P + 1 ) [3]. If one applies δ t → τ δ t ′ in the mapping (1) , then (2) X ¯ i ( P + τ ) = F ¯ i ( x ( P ) , τ δ t ′ ) = ∑ k = 0 N ( τ δ t ′ ) k k ! X k [ x i ( P ) ] . This application enables the choose for different times of prediction, e.g., in chaotic time series [4]. So, other entries beyond the nearest neighbor to the last known value can be predicted and analyzed. Here, we extend the capacity of forecasting and analysis of the computational routines by means the optional argument called PT . It corresponds to parameter τ in the mapping (2) . The researcher can include this option in the commands GfiTS , ForecasTS and IforecasTS . In the NIforecasTS command, our package was already trying to tackle a similar calculation. But it employs the global mapping corresponding to the next entry in the time series, i.e., only the parameter τ = 1 take part in the global fitting. The program predicts the N steps by the next step (corresponding to N = 1 ) and using it as the N = 2 data and so on, up to the actual value of N that we are trying to forecast. Because this different conception, the routines NIforecasTS , AnalysTS and GrafiTS have not changed with respect to their running. The use of the commands in this version of the TimeS package closely follows the previous version. But the programming logistic in the GifTS routine required an update for the correct running of the new features in the current versions of the Maple software. We considered our internal procedure gerpoly deprecated–it does not work properly in the Maple 17 release–and it has been deleted in this version. The polynomials are now generated inside the GfiTS routine. There is a change in the optional argument that selects the part of the time series for the global mapping. The choosing now refers to the last vector to be used in the global mapping. The argument IniPoint was substituted by Final for all routines in this update. The goal of this new option is to relate more easily the global mapping with the position on the time series of interest in the forecasting. Another slight modification is present for the error calculation by the optional argument Poptions . Instead of the percent error, this analysis option now prints the actual error. Let us consider the same time series of the original paper [1]. It is stored in the file ‘ts37.txt’ and corresponds to the dynamical variable X of the Lorenz System [5]. Below, we present the Maple worksheet for the reconstruction of phase space, the global mapping, the forecasting and the improvement of the forecast. The respective routines in this example are VecTS , GfiTS , ForecasTS and IforecasTS . The date to be forecast is dat[402] . − 0.03991323747 X 1 2 + 0.1157254791 X 1 X 2 − 0.05432130896 X 1 X 3 − 0.03977733053 X 2 2 + 0.03040753173 X 2 X 3 − 0.005038168581 X 3 2 + 1.487438057 X 1 − 0.3816889428 X 2 + 0.09738575830 X 3 9.2926021229.231602460 9.2926021229.278554589 . Thus, for the forecasting and analysis of the second entry in a given time series (i.e., τ = 2 ), the argument is PT=2 and so on. In order to compare the same prediction by the command NIforecasTS , we put its respective prompt too: 9.333733567 . The argument Nsteps=2 specifies the parameter τ = 2 and the map Mapag[1] corresponds to τ = 1 above. In this application, the prediction with the new feature in the IforecasTS is more accurate–9.279 compared to 9.333 for the true value 9.293–than the result of the command NIforecasTS . Besides being an alternative for improved accuracy, the new feature enlarges the possibilities in time series analysis by the TimeS package. Summary of revisions: Modification of the way the polynomials needed are generated, making the running of the program compatible with the new version (release 17 up) Maple; The introduction of the possibility of the calculating for the prediction N steps ahead that improved our previous similar calculation; Instructions to use the TimeS package (README.pdf), computational routines (TimeS.txt) and test file (TimeS.mw). Restrictions: If the time series that is being analyzed presents a great amount of noise or if the dynamical system behind the time series is of high dimensionality ( D i m ≫ 3 ), then the method may not work well. Unusual features: Our implementation can, in the cases where the dynamics behind the time series is given by a system of low dimensionality, greatly improve the forecast. Running time: It depends strongly on the command used. [1] H. Carli, L. Duarte, L. da Mota, A maple package for improved global mapping forecast, Computer Physics Communications 185 (3) (2014) 1115–1129. [2] N.H. Packard, J.P. Crutchfield, J.D. Farmer, R.S. Shaw, Geometry from a time series, Phys. Rev. Lett. 45 (1980) 712–716. [3] H. Kantz, T. Schreiber, Nonlinear Time Series Analysis, Cambridge nonlinear science series, Cambridge University Press, 2004. [4] J.D. Farmer, J.J. Sidorowich, Predicting chaotic time series, Phys. Rev. Lett. 59 (1987) 845–848. [5] E. Lorenz, Deterministic nonperiodic flow, Journal of the Atmospheric Sciences 20 (2) (1963) 130–141. [ABSTRACT FROM AUTHOR]
- Published
- 2016
- Full Text
- View/download PDF
48. Reconstruction of the Disassembly Pathway of an Icosahedral Viral Capsid and Shape Determination of Two Successive Intermediates
- Author
-
Mehdi Zeghal, Guillaume Tresset, Prabal K. Maiti, Didier Law-Hine, Sylvain Prévost, Stéphane Bressanelli, Virginie Bailleux, Anil Kumar Sahoo, Doru Constantin, Laboratoire de Physique des Solides (LPS), Centre National de la Recherche Scientifique (CNRS)-Université Paris-Sud - Paris 11 (UP11), Department of Physics, Indian Institute of Science, European Synchrotron Radiation Facility (ESRF), Institut de Biologie Intégrative de la Cellule (I2BC), Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Université Paris-Saclay-Centre National de la Recherche Scientifique (CNRS), ANR-10-LABX-0039,PALM,Physics: Atoms, Light, Matter(2010), Laboratoire de Physique des Solides (UMR CNRS 8502 - Université Paris-Sud) ( LPS ), European Synchrotron Radiation Facility ( ESRF ), Institut de Biologie Intégrative de la Cellule ( I2BC ), Université Paris-Saclay-Centre National de la Recherche Scientifique ( CNRS ) -Commissariat à l'énergie atomique et aux énergies alternatives ( CEA ) -Université Paris-Sud - Paris 11 ( UP11 ), ANR-10-LABX-0039-PALM,LabEx PALM,Laboratoire d’Excellence Physics Atoms Light Mater, and Université Paris-Sud - Paris 11 (UP11)-Centre National de la Recherche Scientifique (CNRS)
- Subjects
Icosahedral symmetry ,[SDV]Life Sciences [q-bio] ,kinetic pathway ,law.invention ,Capsid ,X-Ray Diffraction ,Ab initio quantum chemistry methods ,law ,Icosahedral viral capsid ,Scattering, Small Angle ,RNA Viruses ,General Materials Science ,Physical and Theoretical Chemistry ,[ SDV ] Life Sciences [q-bio] ,Chemistry ,Physics ,Virion ,Global fitting ,modeling ,disassembly ,Synchrotron ,Virus ,Crystallography ,Biophysics ,time-resolved small-angle X-ray scattering ,Spatiotemporal resolution - Abstract
International audience; Viral capsids derived from an icosahedral plant virus widely used in physical and nanotechnological investigations were fully dissociated into dimers by a rapid change of pH. The process was probed in vitro at high spatiotemporal resolution by time-resolved small-angle X-ray scattering using a high brilliance synchrotron source. A powerful custom-made global fitting algorithm allowed us to reconstruct the most likely pathway parametrized by a set of stoichiometric coefficients and to determine the shape of two successive intermediates by ab initio calculations. None of these two unexpected intermediates was previously identified in self-assembly experiments, which suggests that the disassembly pathway is not a mirror image of the assembly pathway. These findings shed new light on the mechanisms and the reversibility of the assembly/disassembly of natural and synthetic virus-based systems. They also demonstrate that both the structure and dynamics of an increasing number of intermediate species become accessible to experiments.
- Published
- 2015
- Full Text
- View/download PDF
49. 'Dunham-type' coefficients for the ground state of 14N16O molecule: global fitting of the experimental energy levels
- Author
-
O. N. Sulakshina and Yu. G. Borkov
- Subjects
Physics ,symbols.namesake ,symbols ,Global fitting ,Molecule ,Atomic physics ,Hamiltonian (quantum mechanics) ,Ground state ,Global model ,Diatomic molecule - Abstract
The global fitting of the experimental energy levels of unresolved Λ-doublets for the ground state of 14 N 16 O molecule is done. The dataset of 1789 experimental energy levels covering the 0-35665 cm -1 interval for 23 vibrational states was obtained using the fundamental Rydberg-Ritz combination principle. A global model of analysis with vibrational dependences of the parameters of the effective Hamiltonian was used for the theoretical treatment of the diatomic molecule in 2 Π electronic state. As a result of the fit a set of the "Dunham-type" molecular parameters was obtained. They reproduce the dataset of the experimental energy levels to the precision of the experimental ones. The found set of the parameters was compared with previous set determined by C.Amiot.
- Published
- 2017
- Full Text
- View/download PDF
50. Statistical Mechanics and Kinetics of Amyloid Fibrillation
- Author
-
Chiu Fan Lee, Ya Jing Huang, and Liu Hong
- Subjects
Fibrillation ,Amyloid ,Chemistry ,Kinetics ,medicine ,Global fitting ,Biophysics ,A protein ,Statistical mechanics ,Computational biology ,medicine.symptom ,Amyloid fibril - Abstract
Amyloid fibrillation is a protein self-assembly phenomenon that is intimately related to well-known human neurodegenerative diseases. During the past few decades, striking advances have been achieved in our understanding of the physical origin of this phenomenon and they constitute the contents of this review. Starting from a minimal model of amyloid fibrils, we explore systematically the equilibrium and kinetic aspects of amyloid fibrillation in both dilute and semi-dilute limits. We then incorporate further molecular mechanisms into the analyses. We also discuss the mathematical foundation of kinetic modeling based on chemical mass-action equations, the quantitative linkage with experimental measurements, as well as the procedure to perform global fitting.
- Published
- 2017
- Full Text
- View/download PDF
Catalog
Discovery Service for Jio Institute Digital Library
For full access to our library's resources, please sign in.