Showing 1–35 of 35 results for author: Nobile, A

Notes on Nash modification

Authors: A. Nobile

Abstract: The Nash blowing-up (or modification) of an algebraic variety $X$ is a canonical process that produces a proper, birational morphism $πぱい: X' \to X$ of varieties. It is expected that the singularities of $X'$ will be better than those of $X$. In the mid-1970's, it was proved that in characteristic zero, $πぱい$ is an isomorphism if and only if $X$ is nonsingular, which is false in positive characteristi… ▽ More The Nash blowing-up (or modification) of an algebraic variety $X$ is a canonical process that produces a proper, birational morphism $πぱい: X' \to X$ of varieties. It is expected that the singularities of $X'$ will be better than those of $X$. In the mid-1970's, it was proved that in characteristic zero, $πぱい$ is an isomorphism if and only if $X$ is nonsingular, which is false in positive characteristic. The focus of this article is on several subsequent studies on this subject. Topics covered include: (a) the extension of the mentioned theorem to the case where $X$ is normal, in any characteristic, (b) the introduction and study of Nash modifications of higher order, (c) the case where the variety $X$ is toric, where more precise results can be obtained and (d) desingularization properties of the Nash process. △ Less

Submitted 13 April, 2024; originally announced April 2024.

Comments: 19 pages

MSC Class: 14E15 (Primary); 14B05 (Secondary)

Nonvolatile Tuning of Bragg Structures Using Transparent Phase-Change Materials

Authors: Nicholas A. Nobile, Chuanyu Lian, Hongyi Sun, Yi-Siou Huang, Brian Mills, Cosmin Constantin Popescu, Dennis Callahan, Juejun Hu, Carlos A. Ríos Ocampo, Nathan Youngblood

Abstract: Bragg gratings offer high-performance filtering and routing of light on-chip through a periodic modulation of a waveguide's effective refractive index. Here, we model and experimentally demonstrate the use of Sb2Se3, a nonvolatile and transparent phase-change material, to tune the resonance conditions in two devices which leverage periodic Bragg gratings: a stopband filter and Fabry-Perot cavity.… ▽ More Bragg gratings offer high-performance filtering and routing of light on-chip through a periodic modulation of a waveguide's effective refractive index. Here, we model and experimentally demonstrate the use of Sb2Se3, a nonvolatile and transparent phase-change material, to tune the resonance conditions in two devices which leverage periodic Bragg gratings: a stopband filter and Fabry-Perot cavity. Through simulations, we show that similar refractive indices between silicon and amorphous Sb2Se3 can be used to induce broadband transparency, while the crystalline state can enhance the index contrast in these Bragg devices. Our experimental results show the promise and limitations of this design approach and highlight specific fabrication challenges which need to be addressed in future implementations. △ Less

Submitted 26 June, 2023; originally announced June 2023.

From formal smoothings to geometric smoothings

Authors: Alessandro Nobile

Abstract: Let X be a projective, equidimensional, singular scheme over an algebraically closed field. Then the existence of a geometric smoothing (i.e. a family of deformations of X over a smooth base curve whose generic fibre is smooth) implies the existence of a formal smoothing as defined by Tziolas. In this paper we address the reverse question giving sufficient conditions on X that guarantee the conver… ▽ More Let X be a projective, equidimensional, singular scheme over an algebraically closed field. Then the existence of a geometric smoothing (i.e. a family of deformations of X over a smooth base curve whose generic fibre is smooth) implies the existence of a formal smoothing as defined by Tziolas. In this paper we address the reverse question giving sufficient conditions on X that guarantee the converse, i.e. formal smoothability implies geometric smoothability. This is useful in light of Tziolas' results giving sufficient criteria for the existence of formal smoothings. △ Less

Submitted 16 January, 2023; v1 submitted 25 November, 2022; originally announced November 2022.

MSC Class: 14D15

Journal ref: Rendiconti di Matematica e delle sue Applicazioni, Volume 44 (2023), 181 -- 210

Time-resolved temperature mapping leveraging the strong thermo-optic effect in phase-change devices

Authors: Nicholas A. Nobile, John R. Erickson, Carlos Ríos, Yifei Zhang, Juejun Hu, Steven A. Vitale, Feng Xiong, Nathan Youngblood

Abstract: Optical phase-change materials are highly promising for emerging applications such as tunable metasurfaces, reconfigurable photonic circuits, and non-von Neumann computing. However, these materials typically require both high melting temperatures and fast quenching rates to reversibly switch between their crystalline and amorphous phases, a significant challenge for large-scale integration. Here,… ▽ More Optical phase-change materials are highly promising for emerging applications such as tunable metasurfaces, reconfigurable photonic circuits, and non-von Neumann computing. However, these materials typically require both high melting temperatures and fast quenching rates to reversibly switch between their crystalline and amorphous phases, a significant challenge for large-scale integration. Here, we present an experimental technique which leverages the thermo-optic effect in GST to enable both spatial and temporal thermal measurements of two common electro-thermal microheater designs currently used by the phase-change community. Our approach shows excellent agreement between experimental results and numerical simulations and provides a non-invasive method for rapid characterization of electrically programmable phase-change devices. △ Less

Submitted 14 October, 2022; originally announced October 2022.

M/M/1 queue in two alternating environments and its heavy traffic approximation

Authors: Antonio Di Crescenzo, Virginia Giorno, Balasubramanian Krishna Kumar, Amelia G. Nobile

Abstract: We investigate an M/M/1 queue operating in two switching environments, where the switch is governed by a two-state time-homogeneous Markov chain. This model allows to describe a system that is subject to regular operating phases alternating with anomalous working phases or random repairing periods. We first obtain the steady-state distribution of the process in terms of a generalized mixture of tw… ▽ More We investigate an M/M/1 queue operating in two switching environments, where the switch is governed by a two-state time-homogeneous Markov chain. This model allows to describe a system that is subject to regular operating phases alternating with anomalous working phases or random repairing periods. We first obtain the steady-state distribution of the process in terms of a generalized mixture of two geometric distributions. In the special case when only one kind of switch is allowed, we analyze the transient distribution, and investigate the busy period problem. The analysis is also performed by means of a suitable heavy-traffic approximation which leads to a continuous random process. Its distribution satisfies a partial differential equation with randomly alternating infinitesimal moments. For the approximating process we determine the steady-state distribution, the transient distribution and a first-passage-time density. △ Less

Submitted 19 March, 2021; originally announced March 2021.

Comments: 31 pages, 18 figures

MSC Class: 60K25; 60K37; 60J60; 60J70

Journal ref: Journal of Mathematical Analysis and Applications, 465 (2018), pp. 973-1001

A continuous-time Ehrenfest model with catastrophes and its jump-diffusion approximation

Authors: Selvamuthu Dharmaraja, Antonio Di Crescenzo, Virginia Giorno, Amelia G. Nobile

Abstract: We consider a continuous-time Ehrenfest model defined over the integers from -N to N, and subject to catastrophes occurring at constant rate. The effect of each catastrophe instantaneously resets the process to state 0. We investigate both the transient and steady-state probabilities of the above model. Further, the first passage time through state 0 is discussed. We perform a jump-diffusion appro… ▽ More We consider a continuous-time Ehrenfest model defined over the integers from -N to N, and subject to catastrophes occurring at constant rate. The effect of each catastrophe instantaneously resets the process to state 0. We investigate both the transient and steady-state probabilities of the above model. Further, the first passage time through state 0 is discussed. We perform a jump-diffusion approximation of the above model, which leads to the Ornstein-Uhlenbeck process with catastrophes. The underlying jump-diffusion process is finally studied, with special attention to the symmetric case arising when the Ehrenfest model has equal upward and downward transition rates. △ Less

Submitted 19 March, 2021; originally announced March 2021.

Comments: 21 pages, 10 figures

MSC Class: 60J80; 60J27; 60J60

Journal ref: Journal of Statistical Physics 161 (2015), pp. 326-345

On the Nash problem for surfaces in positive characteristic

Authors: Augusto Nobile

Abstract: This paper seeks to prove the bijectivity of the "Nash mapping" from the set of irreducible components of the scheme parametrizing analytic arcs on an algebraic surface $X$ whose origin is a singular point, into the set of irreducible components of the exceptional locus of a minimal desingularization $X'$ of $X$ when the base field has positive characteristic. The idea is to view the surface as a… ▽ More This paper seeks to prove the bijectivity of the "Nash mapping" from the set of irreducible components of the scheme parametrizing analytic arcs on an algebraic surface $X$ whose origin is a singular point, into the set of irreducible components of the exceptional locus of a minimal desingularization $X'$ of $X$ when the base field has positive characteristic. The idea is to view the surface as a specialization of another defined over a field of characteristic zero. A number of related results are proved. Among them, the construction of a scheme of arcs for a suitable one parameter family of surfaces which, by using a theorem of M. Artin on lifting of normal surface singularities to characteristic zero, seems a reasonable candidate to be the desired tool. But there are some points, necessary for a complete proof, which are not verified yet. △ Less

Submitted 1 December, 2018; originally announced December 2018.

MSC Class: 14E18; 14E15; 14B05; 14B07; 14E99; 14J99

Resolution Algorithms and Deformations

Authors: Augusto Nobile

Abstract: An algorithm for resolution of singularities in characteristic zero is described. It is expressed in terms of multi-ideals, that essentially are defined as a finite sequence of pairs, each one consiting of a sheaf of ideals and a positive integer. This approach is particularly simple and, as indicated by some results shown here, it seems suitable for applications to a good theory of simultaneous a… ▽ More An algorithm for resolution of singularities in characteristic zero is described. It is expressed in terms of multi-ideals, that essentially are defined as a finite sequence of pairs, each one consiting of a sheaf of ideals and a positive integer. This approach is particularly simple and, as indicated by some results shown here, it seems suitable for applications to a good theory of simultaneous algorithmic resolution of singularities, specially for families parametrized by the spectrum of an artinian ring. △ Less

Submitted 9 April, 2013; originally announced April 2013.

Nucleon axial charge and pion decay constant from two-flavor lattice QCD

Authors: R. Horsley, Y. Nakamura, A. Nobile, P. E. L. Rakow, G. Schierholz, J. M. Zanotti

Abstract: The axial charge of the nucleon $g_A$ and the pion decay constant $f_πぱい$ are computed in two-flavor lattice QCD. The simulations are carried out on lattices of various volumes and lattice spacings. Results are reported for pion masses as low as $m_πぱい=130\,\mbox{MeV}$. Both quantities, $g_A$ and $f_πぱい$, suffer from large finite size effects, which to leading order ChEFT and ChPT turn out to be identic… ▽ More The axial charge of the nucleon $g_A$ and the pion decay constant $f_πぱい$ are computed in two-flavor lattice QCD. The simulations are carried out on lattices of various volumes and lattice spacings. Results are reported for pion masses as low as $m_πぱい=130\,\mbox{MeV}$. Both quantities, $g_A$ and $f_πぱい$, suffer from large finite size effects, which to leading order ChEFT and ChPT turn out to be identical. By considering the naturally renormalized ratio $g_A/f_πぱい$, we observe a universal behavior as a function of decreasing quark mass. From extrapolating the ratio to the physical point, we find $g_A^R=1.29(5)(3)$, using the physical value of $f_πぱい$ as input and $r_0=0.50(1)$ to set the scale. In a subsequent calculation we attempt to extrapolate $g_A$ and $f_πぱい$ separately to the infinite volume. Both volume and quark mass dependencies of $g_A$ and $f_πぱい$ are found to be well decribed by ChEFT and ChPT. We find at the physical point $g_A^R=1.24(4)$ and $f_πぱい^R=89.6(1.1)(1.8)\,\mbox{MeV}$. Both sets of results are in good agreement with experiment. As a by-product we obtain the low-energy constant $\bar{l}_4=4.2(1)$. △ Less

Submitted 1 March, 2014; v1 submitted 9 February, 2013; originally announced February 2013.

Comments: 16 pages, 8 figures; published version (Physics Letters B)

Report number: DESY 13-017, Edinburgh 2013/01, LTH 969

<x>_{u-d} from lattice QCD at nearly physical quark masses

Authors: Gunnar S. Bali, Sara Collins, Mridupawan Deka, Benjamin Glaessle, Meinulf Gockeler, Johannes Najjar, Andrea Nobile, Dirk Pleiter, Andreas Schafer, Andre Sternbeck

Abstract: We determine the second Mellin moment of the isovector quark parton distribution function <x>_{u-d} from lattice QCD with N_f=2 sea quark flavours, employing the non-perturbatively improved Wilson-Sheikholeslami-Wohlert action at a pseudoscalar mass of 157(6) MeV. The result is converted non-perturbatively to the RI'-MOM scheme and then perturbatively to the MSbar scheme at a scale mu = 2 GeV. As… ▽ More We determine the second Mellin moment of the isovector quark parton distribution function <x>_{u-d} from lattice QCD with N_f=2 sea quark flavours, employing the non-perturbatively improved Wilson-Sheikholeslami-Wohlert action at a pseudoscalar mass of 157(6) MeV. The result is converted non-perturbatively to the RI'-MOM scheme and then perturbatively to the MSbar scheme at a scale mu = 2 GeV. As the quark mass is reduced we find the lattice prediction to approach the value extracted from experiments. △ Less

Submitted 11 September, 2012; v1 submitted 4 July, 2012; originally announced July 2012.

Comments: 4 pages, 3 figures, v2: minor updates including journal refs

Nucleon mass and sigma term from lattice QCD with two light fermion flavors

Authors: G. S. Bali, P. C. Bruns, S. Collins, M. Deka, B. Gläßle, M. Göckeler, L. Greil, T. R. Hemmert, R. Horsley, J. Najjar, Y. Nakamura, A. Nobile, D. Pleiter, P. E. L. Rakow, A. Schäfer, R. Schiel, G. Schierholz, A. Sternbeck, J. Zanotti

Abstract: We analyze Nf=2 nucleon mass data with respect to their dependence on the pion mass down to mpi = 157 MeV and compare it with predictions from covariant baryon chiral perturbation theory (BChPT). A novel feature of our approach is that we fit the nucleon mass data simultaneously with the directly obtained pion-nucleon sigma-term. Our lattice data below mpi = 435 MeV is well described by O(p^4) BCh… ▽ More We analyze Nf=2 nucleon mass data with respect to their dependence on the pion mass down to mpi = 157 MeV and compare it with predictions from covariant baryon chiral perturbation theory (BChPT). A novel feature of our approach is that we fit the nucleon mass data simultaneously with the directly obtained pion-nucleon sigma-term. Our lattice data below mpi = 435 MeV is well described by O(p^4) BChPT and we find sigma=37(8)(6) MeV for the sigma-term at the physical point. Using the nucleon mass to set the scale we obtain a Sommer parameter of r_0=0.501(10)(11) fm. △ Less

Submitted 27 August, 2012; v1 submitted 29 June, 2012; originally announced June 2012.

Comments: 26 pages, 11 figures, 5 tables. Version to appear in NPB with a few more details on the fit parameters

Report number: ADP-12-29/T796; DESY 12-105; Edinburgh 2012/11

Journal ref: Nucl. Phys. B866 (2013) 1-25

Deflation and Flexible SAP-Preconditioning of GMRES in Lattice QCD Simulation

Authors: Andreas Frommer, Andrea Nobile, Paul Zingler

Abstract: The simulation of lattice QCD on massively parallel computers stimulated the development of scalable algorithms for the solution of sparse linear systems. We tackle the problem of the Wilson-Dirac operator inversion by combining a Schwarz alternating procedure (SAP) in multiplicative form with a flexible variant of the GMRES-DR algorithm. We show that restarted GMRES is not able to converge when t… ▽ More The simulation of lattice QCD on massively parallel computers stimulated the development of scalable algorithms for the solution of sparse linear systems. We tackle the problem of the Wilson-Dirac operator inversion by combining a Schwarz alternating procedure (SAP) in multiplicative form with a flexible variant of the GMRES-DR algorithm. We show that restarted GMRES is not able to converge when the system is poorly conditioned. By adding deflation in the form of the FGMRES-DR algorithm, an important fraction of the information produced by the iterates is kept between successive restarts leading to convergence in cases in which FGMRES stagnates. △ Less

Submitted 24 April, 2012; originally announced April 2012.

Comments: 26 pages

First moments of the nucleon generalized parton distributions from lattice QCD

Authors: A. Sternbeck, M. Göckeler, Ph. Hägler, R. Horsley, Y. Nakamura, A. Nobile, D. Pleiter, P. E. L. Rakow, A. Schäfer, G. Schierholz, J. Zanotti

Abstract: We report on our lattice calculations of the nucleon's generalized parton distributions (GPDs), concentrating on their first moments for the case of N_f=2. Due to recent progress on the numerical side we are able to present results for the generalized form factors at pion masses as low as 260 MeV. We perform a fit to one-loop covariant baryon chiral perturbation theory with encouraging results. We report on our lattice calculations of the nucleon's generalized parton distributions (GPDs), concentrating on their first moments for the case of N_f=2. Due to recent progress on the numerical side we are able to present results for the generalized form factors at pion masses as low as 260 MeV. We perform a fit to one-loop covariant baryon chiral perturbation theory with encouraging results. △ Less

Submitted 29 March, 2012; originally announced March 2012.

Comments: 7 pages, 4 figures. Contribution to the XXIX International Symposium on Lattice Field Theory, July 10-16, 2011, Squaw Valley, Lake Tahoe, California

Report number: DESY 12-038, Edinburgh 2012/03, LTH 941

Journal ref: PoS(Lattice 2011)177

Strangeness Contribution to the Proton Spin from Lattice QCD

Authors: Gunnar S. Bali, Sara Collins, Meinulf Gockeler, Roger Horsley, Yoshifumi Nakamura, Andrea Nobile, Dirk Pleiter, P. E. L. Rakow, Andreas Schafer, Gerrit Schierholz, James M. Zanotti

Abstract: We compute the strangeness and light-quark contributions Delta s, Delta u and Delta d to the proton spin in nf=2 lattice QCD at a pion mass of about 285 MeV and at a lattice spacing a approx 0.073 fm, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson action. We carry out the renormalization of these matrix elements which involves mixing between contributions from different quark… ▽ More We compute the strangeness and light-quark contributions Delta s, Delta u and Delta d to the proton spin in nf=2 lattice QCD at a pion mass of about 285 MeV and at a lattice spacing a approx 0.073 fm, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson action. We carry out the renormalization of these matrix elements which involves mixing between contributions from different quark flavours. Our main result is the small negative value MSbar scheme value Delta s = -0.020(10)(4), at a scale mu = 7.4 GeV, of the strangeness contribution to the nucleon spin. The second error is an estimate of the uncertainty, due to the missing extrapolation to the physical point. △ Less

Submitted 18 May, 2012; v1 submitted 14 December, 2011; originally announced December 2011.

Comments: 5 pages, 2 figures, V2: minor changes, to enhance the readability of the article. References updated. "The" removed from the title

Report number: Adelaide ADP-11-43/T765; Edinburgh 2011/39; Liverpool LTH 934

Journal ref: Phys. Rev. Lett. 108 (2012) 222001

A lattice study of the strangeness content of the nucleon

Authors: G. S. Bali, S. Collins, M. Gockeler, R. Horsley, Y. Nakamura, A. Nobile, D. Pleiter, P. E. L. Rakow, A. Sternbeck, A. Schafer, G. Schierholz, J. M. Zanotti

Abstract: We determine the quark contributions to the nucleon spin Delta s, Delta u and Delta d as well as their contributions to the nucleon mass, the sigma-terms. This is done by computing both, the quark line connected and disconnected contributions to the respective matrix elements, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson Fermionic action. We simulate n_F=2 mass degenerate se… ▽ More We determine the quark contributions to the nucleon spin Delta s, Delta u and Delta d as well as their contributions to the nucleon mass, the sigma-terms. This is done by computing both, the quark line connected and disconnected contributions to the respective matrix elements, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson Fermionic action. We simulate n_F=2 mass degenerate sea quarks with a pion mass of about 285 MeV and a lattice spacing a = 0.073 fm. The renormalization of the matrix elements involves mixing between contributions from different quark flavours. The pion-nucleon sigma-term is extrapolated to physical quark masses exploiting the sea quark mass dependence of the nucleon mass. We obtain the renormalized value sigma_{piN}=38(12) MeV at the physical point and the strangeness fraction f_{Ts}=sigma_s/m_N=0.012(14)(+10-3) at our larger than physical sea quark mass. For the strangeness contribution to the nucleon spin we obtain in the MSbar scheme at the renormalization scale of 2.71 GeV Delta s = -0.020(10)(2). △ Less

Submitted 30 November, 2011; originally announced December 2011.

Comments: 7 pages, 3 figures, Invited Talk at the 33rd Erice School on Nuclear Physics, Erice, 16-24 September 2011, Italy

Journal ref: Prog. Part. Nucl. Phys. 67 (2012) 467

The strange and light quark contributions to the nucleon mass from Lattice QCD

Authors: Gunnar S. Bali, Sara Collins, Meinulf Gockeler, Roger Horsley, Yoshifumi Nakamura, Andrea Nobile, Dirk Pleiter, P. E. L. Rakow, Andreas Schafer, Gerrit Schierholz, Andre Sternbeck, James M. Zanotti

Abstract: We determine the strangeness and light quark fractions of the nucleon mass by computing the quark line connected and disconnected contributions to the matrix elements m_q <N|qbar q|N> in lattice QCD, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson Fermionic action. We simulate n_F=2 mass degenerate sea quarks with a pion mass of about 285 MeV and a lattice spacing a approx 0.07… ▽ More We determine the strangeness and light quark fractions of the nucleon mass by computing the quark line connected and disconnected contributions to the matrix elements m_q <N|qbar q|N> in lattice QCD, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson Fermionic action. We simulate n_F=2 mass degenerate sea quarks with a pion mass of about 285 MeV and a lattice spacing a approx 0.073 fm. The renormalization of the matrix elements involves mixing between contributions from different quark flavours. The pion-nucleon sigma-term is extrapolated to physical quark masses exploiting the sea quark mass dependence of the nucleon mass. We obtain the renormalized values σしぐま_{πぱいN} = 38(12) MeV at the physical point and f_{T_s}=σしぐま_s/m_N= 0.012(14)^{+10}_{-3} for the strangeness contribution at our larger than physical sea quark mass. △ Less

Submitted 27 February, 2012; v1 submitted 7 November, 2011; originally announced November 2011.

Comments: 11 pages, 7 figures, v2: 1 reference added, some corrections

Report number: Adelaide ADP-11-33/T755; Edinburgh 2011/33; Liverpool LTH 930

Journal ref: Phys. Rev. D 85 (2012) 054502

Solving the Dirac equation on QPACE

Authors: Andrea Nobile

Abstract: We discuss the implementation and optimization challenges for a Wilson-Dirac solver with Clover term on QPACE, a parallel machine based on Cell processors and a torus network. We choose the mixed-precision Schwarz preconditioned FGCR algorithm in order to circumvent network bandwidth and latency constraints, to make efficient use of the multicore parallelism and on-chip memory, and to achieve flex… ▽ More We discuss the implementation and optimization challenges for a Wilson-Dirac solver with Clover term on QPACE, a parallel machine based on Cell processors and a torus network. We choose the mixed-precision Schwarz preconditioned FGCR algorithm in order to circumvent network bandwidth and latency constraints, to make efficient use of the multicore parallelism and on-chip memory, and to achieve flexibility in the choice of lattice sizes. We present benchmarks on up to 256 QPACE nodes showing an aggregate sustained performance of about 10 TFlops for the complete solver and very good scaling. △ Less

Submitted 20 September, 2011; originally announced September 2011.

Comments: 7 pages, Lattice 2010 Proceedings

Dirac and Pauli form factors from lattice QCD

Authors: S. Collins, M. Göckeler, Ph. Hägler, R. Horsley, Y. Nakamura, A. Nobile, D. Pleiter, P. E. L. Rakow, A. Schäfer, G. Schierholz, W. Schroers, H. Stüben, F. Winter, J. M. Zanotti

Abstract: We present a comprehensive analysis of the electromagnetic form factors of the nucleon from a lattice simulation with two flavors of dynamical O(a)-improved Wilson fermions. A key feature of our calculation is that we make use of an extensive ensemble of lattice gauge field configurations with four different lattice spacings, multiple volumes, and pion masses down to m_πぱい~ 180 MeV. We find that by… ▽ More We present a comprehensive analysis of the electromagnetic form factors of the nucleon from a lattice simulation with two flavors of dynamical O(a)-improved Wilson fermions. A key feature of our calculation is that we make use of an extensive ensemble of lattice gauge field configurations with four different lattice spacings, multiple volumes, and pion masses down to m_πぱい~ 180 MeV. We find that by employing Kelly-inspired parametrizations for the Q^2-dependence of the form factors, we are able to obtain stable fits over our complete ensemble. Dirac and Pauli radii and the anomalous magnetic moments of the nucleon are extracted and results at light quark masses provide evidence for chiral non-analytic behavior in these fundamental observables. △ Less

Submitted 17 June, 2011; originally announced June 2011.

Comments: 31 pages, 62 figures

Report number: DESY 11-102, Edinburgh 2011/18, MKPH-T-11-13

Lattice QCD Applications on QPACE

Authors: Y. Nakamura, A. Nobile, D. Pleiter, H. Simma, T. Streuer, T. Wettig, F. Winter

Abstract: QPACE is a novel massively parallel architecture optimized for lattice QCD simulations. A single QPACE node is based on the IBM PowerXCell 8i processor. The nodes are interconnected by a custom 3-dimensional torus network implemented on an FPGA. The compute power of the processor is provided by 8 Synergistic Processing Units. Making efficient use of these accelerator cores in scientific applicatio… ▽ More QPACE is a novel massively parallel architecture optimized for lattice QCD simulations. A single QPACE node is based on the IBM PowerXCell 8i processor. The nodes are interconnected by a custom 3-dimensional torus network implemented on an FPGA. The compute power of the processor is provided by 8 Synergistic Processing Units. Making efficient use of these accelerator cores in scientific applications is challenging. In this paper we describe our strategies for porting applications to the QPACE architecture and report on performance numbers. △ Less

Submitted 7 March, 2011; originally announced March 2011.

Comments: 10 pages, 4 figures, accepted for presentation at ICCS 2011 (Tsukuba)

A double-ended queue with catastrophes and repairs, and a jump-diffusion approximation

Authors: Antonio Di Crescenzo, Virginia Giorno, Balasubramanian Krishna Kumar, Amelia G. Nobile

Abstract: Consider a system performing a continuous-time random walk on the integers, subject to catastrophes occurring at constant rate, and followed by exponentially-distributed repair times. After any repair the system starts anew from state zero. We study both the transient and steady-state probability laws of the stochastic process that describes the state of the system. We then derive a heavy-traffic… ▽ More Consider a system performing a continuous-time random walk on the integers, subject to catastrophes occurring at constant rate, and followed by exponentially-distributed repair times. After any repair the system starts anew from state zero. We study both the transient and steady-state probability laws of the stochastic process that describes the state of the system. We then derive a heavy-traffic approximation to the model that yields a jump-diffusion process. The latter is equivalent to a Wiener process subject to randomly occurring jumps, whose probability law is obtained. The goodness of the approximation is finally discussed. △ Less

Submitted 26 January, 2011; originally announced January 2011.

Comments: 18 pages, 5 figures, paper accepted by "Methodology and Computing in Applied Probability", the final publication is available at http://www.springerlink.com

MSC Class: 60J80; 60J25

Nucleon form factors and structure functions

Authors: S. Collins, M. Göckeler, Ph. Hägler, T. Hemmert, R. Horsley, Y. Nakamura, A. Nobile, H. Perlt, D. Pleiter, P. E. L. Rakow, A. Schäfer, G. Schierholz, A. Sternbeck, H. Stüben, F. Winter, J. M. Zanotti

Abstract: We give an update on our ongoing efforts to compute the nucleon's form factors and moments of structure functions using Nf=2 flavours of non-perturbatively improved Clover fermions. We focus on new results obtained on gauge configurations where the pseudo-scalar meson mass is in the range of 170-270 MeV. We will compare our results with various estimates obtained from chiral effective theories sin… ▽ More We give an update on our ongoing efforts to compute the nucleon's form factors and moments of structure functions using Nf=2 flavours of non-perturbatively improved Clover fermions. We focus on new results obtained on gauge configurations where the pseudo-scalar meson mass is in the range of 170-270 MeV. We will compare our results with various estimates obtained from chiral effective theories since we have some overlap with the quark mass region where results from such theories are believed to be applicable. △ Less

Submitted 12 January, 2011; originally announced January 2011.

Comments: Contribution to Lattice2010, Villasimius, Italy

Report number: DESY 10-204, Edinburgh 2010/30

Disconnected contributions to hadronic structure

Authors: Sara Collins, Gunnar Bali, Andrea Nobile, Andreas Schafer, Yoshifumi Nakamura, James Zanotti

Abstract: We present an update of an on-going project to determine the disconnected contributions to hadronic structure, specifically, the scalar matrix element and the quark contribution to the spin of the nucleon. We present an update of an on-going project to determine the disconnected contributions to hadronic structure, specifically, the scalar matrix element and the quark contribution to the spin of the nucleon. △ Less

Submitted 9 November, 2010; originally announced November 2010.

Comments: 7 pages, 7 figures, talk given at 28th International Symposium on Lattice Field Theory (Lattice 2010), Villasimius, Italy, 14-19 June 2010

Journal ref: PoS LATTICE 2010: 134, 2010

Algorithmic equiresolution of possibly non-reduced families of singularities

Authors: Augusto Nobile

Abstract: This paper studies the concept of algorithmic equiresolution of a family of embedded varieties or ideals, which means a simultaneous resolution of such a family compatible with a given (suitable) algorithm of resolution in characteristic zero. The paper's approach is more indirect: it primarily considers the more general case of families of basic objects (or marked ideals). A definition of algori… ▽ More This paper studies the concept of algorithmic equiresolution of a family of embedded varieties or ideals, which means a simultaneous resolution of such a family compatible with a given (suitable) algorithm of resolution in characteristic zero. The paper's approach is more indirect: it primarily considers the more general case of families of basic objects (or marked ideals). A definition of algorithmic equiresolution is proposed, which applies to families whose parameter space T may be non-reduced, e.g., the spectrum of a suitable artinian ring. Other definitions of algorithmic equiresolution are also discussed. These are geometrically very natural, but the parameter space T of the family must be assumed regular. It is proven that when T is regular all the proposed definitions are equivalent. △ Less

Submitted 5 May, 2010; originally announced May 2010.

QPACE -- a QCD parallel computer based on Cell processors

Authors: H. Baier, H. Boettiger, M. Drochner, N. Eicker, U. Fischer, Z. Fodor, A. Frommer, C. Gomez, G. Goldrian, S. Heybrock, D. Hierl, M. Hüsken, T. Huth, B. Krill, J. Lauritsen, T. Lippert, T. Maurer, B. Mendl, N. Meyer, A. Nobile, I. Ouda, M. Pivanti, D. Pleiter, M. Ries, A. Schäfer , et al. (10 additional authors not shown)

Abstract: QPACE is a novel parallel computer which has been developed to be primarily used for lattice QCD simulations. The compute power is provided by the IBM PowerXCell 8i processor, an enhanced version of the Cell processor that is used in the Playstation 3. The QPACE nodes are interconnected by a custom, application optimized 3-dimensional torus network implemented on an FPGA. To achieve the very hig… ▽ More QPACE is a novel parallel computer which has been developed to be primarily used for lattice QCD simulations. The compute power is provided by the IBM PowerXCell 8i processor, an enhanced version of the Cell processor that is used in the Playstation 3. The QPACE nodes are interconnected by a custom, application optimized 3-dimensional torus network implemented on an FPGA. To achieve the very high packaging density of 26 TFlops per rack a new water cooling concept has been developed and successfully realized. In this paper we give an overview of the architecture and highlight some important technical details of the system. Furthermore, we provide initial performance results and report on the installation of 8 QPACE racks providing an aggregate peak performance of 200 TFlops. △ Less

Submitted 23 December, 2009; v1 submitted 11 November, 2009; originally announced November 2009.

Comments: 21 pages. Poster by T. Maurer and plenary talk by D. Pleiter presented at the "XXVII International Symposium on Lattice Field Theory", July 26-31 2009, Peking University, Beijing, China. Information on recent Green500 ranking added and list of authors extended

Journal ref: PoS LAT2009:001,2009

Algorithmic equiresolution of deformations of embedded varieties

Authors: Augusto Nobile

Abstract: A theory of simultaneous resolution of singularities for families of embedded varieties (over a field of characteristic zero) parametrized by the spectrum of a suitable artinian ring, and compatible with a given algorithm of resolution, is presented. As usually, this a simple consequence of a similar theory for analogous families of basic objects, to which the main portion of this article is dev… ▽ More A theory of simultaneous resolution of singularities for families of embedded varieties (over a field of characteristic zero) parametrized by the spectrum of a suitable artinian ring, and compatible with a given algorithm of resolution, is presented. As usually, this a simple consequence of a similar theory for analogous families of basic objects, to which the main portion of this article is devoted. △ Less

Submitted 23 April, 2009; originally announced April 2009.

MSC Class: 14B05;14D99

Status of the QPACE Project

Authors: H. Baier, H. Boettiger, M. Drochner, N. Eicker, U. Fischer, Z. Fodor, G. Goldrian, S. Heybrock, D. Hierl, T. Huth, B. Krill, J. Lauritsen, T. Lippert, T. Maurer, J. McFadden, N. Meyer, A. Nobile, I. Ouda, M. Pivanti, D. Pleiter, A. Schäfer, H. Schick, F. Schifano, H. Simma, S. Solbrig , et al. (5 additional authors not shown)

Abstract: We give an overview of the QPACE project, which is pursuing the development of a massively parallel, scalable supercomputer for LQCD. The machine is a three-dimensional torus of identical processing nodes, based on the PowerXCell 8i processor. The nodes are connected by an FPGA-based, application-optimized network processor attached to the PowerXCell 8i processor. We present a performance analys… ▽ More We give an overview of the QPACE project, which is pursuing the development of a massively parallel, scalable supercomputer for LQCD. The machine is a three-dimensional torus of identical processing nodes, based on the PowerXCell 8i processor. The nodes are connected by an FPGA-based, application-optimized network processor attached to the PowerXCell 8i processor. We present a performance analysis of lattice QCD codes on QPACE and corresponding hardware benchmarks. △ Less

Submitted 9 October, 2008; originally announced October 2008.

Comments: 7 pages, poster presented at the XXVI International Symposium on Lattice Field Theory, July 14-19, 2008, Williamsburg, Virginia, USA

Journal ref: PoS LATTICE2008:039,2008

Bayesian finite mixtures: a note on prior specification and posterior computation

Authors: Agostino Nobile

Abstract: A new method for the computation of the posterior distribution of the number k of components in a finite mixture is presented. Two aspects of prior specification are also studied: an argument is made for the use of a Poisson(1) distribution as the prior for k; and methods are given for the selection of hyperparameter values in the mixture of normals model, with natural conjugate priors on the co… ▽ More A new method for the computation of the posterior distribution of the number k of components in a finite mixture is presented. Two aspects of prior specification are also studied: an argument is made for the use of a Poisson(1) distribution as the prior for k; and methods are given for the selection of hyperparameter values in the mixture of normals model, with natural conjugate priors on the components parameters. △ Less

Submitted 3 November, 2007; originally announced November 2007.

Comments: 14 pages, 5 figures

Report number: Univ of Glasgow, Dept of Statistics, Tech report 05-3

QCD on the Cell Broadband Engine

Authors: F. Belletti, G. Bilardi, M. Drochner, N. Eicker, Z. Fodor, D. Hierl, H. Kaldass, T. Lippert, T. Maurer, N. Meyer, A. Nobile, D. Pleiter, A. Schaefer, F. Schifano, H. Simma, S. Solbrig, T. Streuer, R. Tripiccione, T. Wettig

Abstract: We evaluate IBM's Enhanced Cell Broadband Engine (BE) as a possible building block of a new generation of lattice QCD machines. The Enhanced Cell BE will provide full support of double-precision floating-point arithmetics, including IEEE-compliant rounding. We have developed a performance model and applied it to relevant lattice QCD kernels. The performance estimates are supported by micro- and… ▽ More We evaluate IBM's Enhanced Cell Broadband Engine (BE) as a possible building block of a new generation of lattice QCD machines. The Enhanced Cell BE will provide full support of double-precision floating-point arithmetics, including IEEE-compliant rounding. We have developed a performance model and applied it to relevant lattice QCD kernels. The performance estimates are supported by micro- and application-benchmarks that have been obtained on currently available Cell BE-based computers, such as IBM QS20 blades and PlayStation 3. The results are encouraging and show that this processor is an interesting option for lattice QCD applications. For a massively parallel machine on the basis of the Cell BE, an application-optimized network needs to be developed. △ Less

Submitted 12 October, 2007; originally announced October 2007.

Comments: 7 pages, 3 figures, contribution to Lattice 2007 (Regensburg, Germany)

Journal ref: PoSLAT2007:039,2007

Study of a model for the folding of a small protein

Authors: Andrea Nobile, Federico Rapuano

Abstract: We describe the results obtained from an improved model for protein folding. We find that a good agreement with the native structure of a 46 residue long, five-letter protein segment is obtained by carefully tuning the parameters of the self-avoiding energy. In particular we find an improved free-energy profile. We also compare the efficiency of the multidimensional replica exchange method with… ▽ More We describe the results obtained from an improved model for protein folding. We find that a good agreement with the native structure of a 46 residue long, five-letter protein segment is obtained by carefully tuning the parameters of the self-avoiding energy. In particular we find an improved free-energy profile. We also compare the efficiency of the multidimensional replica exchange method with the widely used parallel tempering. △ Less

Submitted 30 November, 2005; v1 submitted 22 March, 2005; originally announced March 2005.

Comments: typos corrected, one figure added

On the posterior distribution of the number of components in a finite mixture

Authors: Agostino Nobile

Abstract: The posterior distribution of the number of components k in a finite mixture satisfies a set of inequality constraints. The result holds irrespective of the parametric form of the mixture components and under assumptions on the prior distribution weaker than those routinely made in the literature on Bayesian analysis of finite mixtures. The inequality constraints can be used to perform an ``inte… ▽ More The posterior distribution of the number of components k in a finite mixture satisfies a set of inequality constraints. The result holds irrespective of the parametric form of the mixture components and under assumptions on the prior distribution weaker than those routinely made in the literature on Bayesian analysis of finite mixtures. The inequality constraints can be used to perform an ``internal'' consistency check of MCMC estimates of the posterior distribution of k and to provide improved estimates which are required to satisfy the constraints. Bounds on the posterior probability of k components are derived using the constraints. Implications on prior distribution specification and on the adequacy of the posterior distribution of k as a tool for selecting an adequate number of components in the mixture are also explored. △ Less

Submitted 29 March, 2005; originally announced March 2005.

Comments: Published at http://dx.doi.org/10.1214/009053604000000788 in the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AOS-AOS241 MSC Class: 62F15. (Primary)

Journal ref: Annals of Statistics 2004, Vol. 32, No. 5, 2044-2073

Simulations of Gaussian Processes and Neuronal Modeling

Authors: Elvira Di Nardo, Amelia G. Nobile, Enrica Pirozzi, Luigi M. Ricciardi

Abstract: The research work outlined in the present note highlights the essential role played by the simulation procedures implemented by us on CINECA supercomputers to complement the mathematical investigations carried within our group over the past several years. The ultimate target of our research is the understanding of certain crucial features of the information processing and transmission by single… ▽ More The research work outlined in the present note highlights the essential role played by the simulation procedures implemented by us on CINECA supercomputers to complement the mathematical investigations carried within our group over the past several years. The ultimate target of our research is the understanding of certain crucial features of the information processing and transmission by single neurons embedded in complex networks. More specifically, here we provide a bird's eye look of some analytical, numerical and simulation results on the asymptotic behavior of first passage time densities for Gaussian processes, both of a Markov and of a non-Markov type. Several figures indicate significant similarities or diversities between computational and simulated results. △ Less

Submitted 3 December, 2004; v1 submitted 1 December, 2004; originally announced December 2004.

Comments: Extended version of the paper published in Science and Supercomputing at CINECA - Report 2003 (Garofalo F., Moretti M., Voli M., eds.), 375-381. ISBN 88-86037-13-9

MSC Class: 60G15; 60G10; 65C50

On the first-visit-time problem for birth and death processes with catastrophes

Authors: A. Di Crescenzo, V. Giorno, A. G. Nobile, L. M. Ricciardi

Abstract: For a birth-death process subject to catastrophes, defined on the state-space $S=\{r,r+1,r+2,...\}$, with $r$ a positive integer or zero, the first-visit time to a state $k\in S$ is considered and the Laplace transform of its probability density function is determined, use of which is then made to obtain mean and variance. The Laplace transform of the probability density function of the first ef… ▽ More For a birth-death process subject to catastrophes, defined on the state-space $S=\{r,r+1,r+2,...\}$, with $r$ a positive integer or zero, the first-visit time to a state $k\in S$ is considered and the Laplace transform of its probability density function is determined, use of which is then made to obtain mean and variance. The Laplace transform of the probability density function of the first effective catastrophe occurrence time and its expected value are also obtained. Some extensions to time-non-homogeneous processes are then provided. Finally, certain additional results concerning the determination of the steady-state distribution and the representation of the transition probabilities are worked out, while some applications to particular birth-death processes are shown in the Appendix. △ Less

Submitted 15 July, 2003; originally announced July 2003.

Comments: 23 pages, no figures

MSC Class: 60J27; 60J80

Towards the Modeling of Neuronal Firing by Gaussian Processes

Authors: E. Di Nardo, A. G. Nobile, E. Pirozzi, L. M. Ricciardi

Abstract: This paper focuses on the outline of some computational methods for the approximate solution of the integral equations for the neuronal firing probability density and an algorithm for the generation of sample-paths in order to construct histograms estimating the firing densities. Our results originate from the study of non-Markov stationary Gaussian neuronal models with the aim to determine the… ▽ More This paper focuses on the outline of some computational methods for the approximate solution of the integral equations for the neuronal firing probability density and an algorithm for the generation of sample-paths in order to construct histograms estimating the firing densities. Our results originate from the study of non-Markov stationary Gaussian neuronal models with the aim to determine the neuron's firing probability density function. A parallel algorithm has been implemented in order to simulate large numbers of sample paths of Gaussian processes characterized by damped oscillatory covariances in the presence of time dependent boundaries. The analysis based on the simulation procedure provides an alternative research tool when closed-form results or analytic evaluation of the neuronal firing densities are not available. △ Less

Submitted 30 May, 2003; originally announced May 2003.

Comments: 10 pages, 3 figures, to be published in Scientiae Mathematicae Japonicae

MSC Class: 60G15; 60G10; 92C20; 68U20; 65C50

On the asymptotic behavior of first passage time densities for stationary Gaussian processes

Authors: E. Di Nardo, A. G. Nobile, E. Pirozzi, L. M. Ricciardi

Abstract: Making use of a Rice-like series expansion, for a class of stationary Gaussian processes the asymptotic behavior of the first passage time probability density function through certain time-varying boundaries, including periodic boundaries, is determined. Sufficient conditions are then given such that the density asymptotically exhibits an exponential behavior when the boundary is either asymptot… ▽ More Making use of a Rice-like series expansion, for a class of stationary Gaussian processes the asymptotic behavior of the first passage time probability density function through certain time-varying boundaries, including periodic boundaries, is determined. Sufficient conditions are then given such that the density asymptotically exhibits an exponential behavior when the boundary is either asymptotically constant or asymptotically periodic. △ Less

Submitted 30 May, 2003; v1 submitted 16 May, 2003; originally announced May 2003.

Comments: 21 pages, 7 figures, to be published in Methodology and Computing in Applied Probability

MSC Class: 60G15; 60G10; 60G40

On Algorithmic Equiresolution and Stratification of Hilbert Schemes

Authors: S. Encinas, A. Nobile, O. Villamayor

Abstract: Given an algorithm of resolution of singularities satisfying certain conditions (``good algorithms''), natural notions of simultaneous algorithmic resolution, or equiresolution, for families of embedded schemes (parametrized by a reduced scheme $T$) are proposed. It is proved that these conditions are equivalent. Something similar is done for families of sheaves of ideals, here the goal is algor… ▽ More Given an algorithm of resolution of singularities satisfying certain conditions (``good algorithms''), natural notions of simultaneous algorithmic resolution, or equiresolution, for families of embedded schemes (parametrized by a reduced scheme $T$) are proposed. It is proved that these conditions are equivalent. Something similar is done for families of sheaves of ideals, here the goal is algorithmic simultaneous principalization. A consequence is that given a family of embedded schemes over a reduced $T$, this parameter scheme can be naturally expressed as a disjoint union of locally closed sets $T_{j}$, such that the induced family on each part $T_{j}$ is equisolvable. In particular, this can be applied to the Hilbert scheme of a smooth projective variety; in fact, our result shows that, in characteristic zero, the underlying topological space of any Hilbert scheme parametrizing embedded schemes can be naturally stratified in equiresolvable families. △ Less

Submitted 25 October, 2000; originally announced October 2000.

Comments: Plain TeX, 37 pages

MSC Class: 14E15; 14D99

Journal ref: Proc. London Math. Soc. (3) 86 (2003) 607-648