Showing 1–48 of 48 results for author: Arnold, N

Performance of Higher-Order Networks in Reconstructing Sequential Paths: from Micro to Macro Scale

Authors: Kevin Teo, Naomi Arnold, Andrew Hone, István Zoltán Kiss

Abstract: Activities such as the movement of passengers and goods, the transfer of physical or digital assets, web navigation and even successive passes in football, result in timestamped paths through a physical or virtual network. The need to analyse such paths has produced a new modelling paradigm in the form of higher-order networks which are able to capture temporal and topological characteristics of s… ▽ More Activities such as the movement of passengers and goods, the transfer of physical or digital assets, web navigation and even successive passes in football, result in timestamped paths through a physical or virtual network. The need to analyse such paths has produced a new modelling paradigm in the form of higher-order networks which are able to capture temporal and topological characteristics of sequential data. This has been complemented by sequence mining approaches, a key example being sequential motifs measuring the prevalence of recurrent subsequences. Previous work on higher-order networks has focused on how to identify the optimal order for a path dataset, where the order can be thought of as the number of steps of memory encoded in the model. In this paper, we build on these approaches to consider which orders are necessary to reproduce different path characteristics, from path lengths to counts of sequential motifs, viewing paths generated from different higher-order models as null models which capture features of the data up to a certain order, and randomised otherwise. Furthermore, we provide an important extension to motif counting, whereby cases with self-loops, starting nodes, and ending nodes of paths are taken into consideration. Conducting a thorough analysis using path lengths and sequential motifs on a diverse range of path datasets, we show that our approach can shed light on precisely where models of different order overperform or underperform, and what this may imply about the original path data. △ Less

Submitted 29 July, 2024; originally announced July 2024.

Comments: 23 pages, 8 figures, submitted to the Journal of Complex Networks

Insights and caveats from mining local and global temporal motifs in cryptocurrency transaction networks

Authors: Naomi A. Arnold, Peijie Zhong, Cheick Tidiane Ba, Ben Steer, Raul Mondragon, Felix Cuadrado, Renaud Lambiotte, Richard G. Clegg

Abstract: Distributed ledger technologies have opened up a wealth of fine-grained transaction data from cryptocurrencies like Bitcoin and Ethereum. This allows research into problems like anomaly detection, anti-money laundering, pattern mining and activity clustering (where data from traditional currencies is rarely available). The formalism of temporal networks offers a natural way of representing this da… ▽ More Distributed ledger technologies have opened up a wealth of fine-grained transaction data from cryptocurrencies like Bitcoin and Ethereum. This allows research into problems like anomaly detection, anti-money laundering, pattern mining and activity clustering (where data from traditional currencies is rarely available). The formalism of temporal networks offers a natural way of representing this data and offers access to a wealth of metrics and models. However, the large scale of the data presents a challenge using standard graph analysis techniques. We use temporal motifs to analyse two Bitcoin datasets and one NFT dataset, using sequences of three transactions and up to three users. We show that the commonly used technique of simply counting temporal motifs over all users and all time can give misleading conclusions. Here we also study the motifs contributed by each user and discover that the motif distribution is heavy-tailed and that the key players have diverse motif signatures. We study the motifs that occur in different time periods and find events and anomalous activity that cannot be seen just by a count on the whole dataset. Studying motif completion time reveals dynamics driven by human behaviour as well as algorithmic behaviour. △ Less

Submitted 14 February, 2024; originally announced February 2024.

Raphtory: The temporal graph engine for Rust and Python

Authors: Ben Steer, Naomi Arnold, Cheick Tidiane Ba, Renaud Lambiotte, Haaroon Yousaf, Lucas Jeub, Fabian Murariu, Shivam Kapoor, Pedro Rico, Rachel Chan, Louis Chan, James Alford, Richard G. Clegg, Felix Cuadrado, Matthew Russell Barnes, Peijie Zhong, John N. Pougué Biyong, Alhamza Alnaimi

Abstract: Raphtory is a platform for building and analysing temporal networks. The library includes methods for creating networks from a variety of data sources; algorithms to explore their structure and evolution; and an extensible GraphQL server for deployment of applications built on top. Raphtory's core engine is built in Rust, for efficiency, with Python interfaces, for ease of use. Raphtory is develop… ▽ More Raphtory is a platform for building and analysing temporal networks. The library includes methods for creating networks from a variety of data sources; algorithms to explore their structure and evolution; and an extensible GraphQL server for deployment of applications built on top. Raphtory's core engine is built in Rust, for efficiency, with Python interfaces, for ease of use. Raphtory is developed by network scientists, with a background in Physics, Applied Mathematics, Engineering and Computer Science, for use across academia and industry. △ Less

Submitted 3 January, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Non-Markovian paths and cycles in NFT trades

Authors: Haaroon Yousaf, Naomi A. Arnold, Renaud Lambiotte, Timothy LaRock, Richard G. Clegg, Peijie Zhong, Alhamza Alnaimi, Ben Steer

Abstract: Recent years have witnessed the availability of richer and richer datasets in a variety of domains, where signals often have a multi-modal nature, blending temporal, relational and semantic information. Within this context, several works have shown that standard network models are sometimes not sufficient to properly capture the complexity of real-world interacting systems. For this reason, differ… ▽ More Recent years have witnessed the availability of richer and richer datasets in a variety of domains, where signals often have a multi-modal nature, blending temporal, relational and semantic information. Within this context, several works have shown that standard network models are sometimes not sufficient to properly capture the complexity of real-world interacting systems. For this reason, different attempts have been made to enrich the network language, leading to the emerging field of higher-order networks. In this work, we investigate the possibility of applying methods from higher-order networks to extract information from the online trade of Non-fungible tokens (NFTs), leveraging on their intrinsic temporal and non-Markovian nature. While NFTs as a technology open up the realms for many exciting applications, its future is marred by challenges of proof of ownership, scams, wash trading and possible money laundering. We demonstrate that by investigating time-respecting non-Markovian paths exhibited by NFT trades, we provide a practical path-based approach to fraud detection. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: Abstract accepted as a workshop presentation at Complex Data Blocks, a satellite of the 2022 Conference of Complex Systems

Using a Bayesian approach to reconstruct graph statistics after edge sampling

Authors: Naomi A. Arnold, Raul J. Mondragon, Richard G. Clegg

Abstract: Often, due to prohibitively large size or to limits to data collecting APIs, it is not possible to work with a complete network dataset and sampling is required. A type of sampling which is consistent with Twitter API restrictions is uniform edge sampling. In this paper, we propose a methodology for the recovery of two fundamental network properties from an edge-sampled network: the degree distrib… ▽ More Often, due to prohibitively large size or to limits to data collecting APIs, it is not possible to work with a complete network dataset and sampling is required. A type of sampling which is consistent with Twitter API restrictions is uniform edge sampling. In this paper, we propose a methodology for the recovery of two fundamental network properties from an edge-sampled network: the degree distribution and the triangle count (we estimate the totals for the network and the counts associated with each edge). We use a Bayesian approach and show a range of methods for constructing a prior which does not require assumptions about the original network. Our approach is tested on two synthetic and three real datasets with diverse sizes, degree distributions, degree-degree correlations and triangle count distributions. △ Less

Submitted 26 June, 2023; v1 submitted 24 July, 2022; originally announced July 2022.

Comments: Extended version of the paper accepted in Complex Networks 2022

Local $L^2$-bounded commuting projections in FEEC

Authors: Douglas N. Arnold, Johnny Gúzman

Abstract: We construct local projections into canonical finite element spaces that appear in the finite element exterior calculus. These projections are bounded in $L^2$ and commute with the exterior derivative. We construct local projections into canonical finite element spaces that appear in the finite element exterior calculus. These projections are bounded in $L^2$ and commute with the exterior derivative. △ Less

Submitted 31 March, 2021; originally announced April 2021.

MSC Class: 65N30

The landscape law for tight binding Hamiltonians

Authors: Douglas N. Arnold, Marcel Filoche, Svitlana Mayboroda, Wei Wang, Shiwen Zhang

Abstract: The present paper extends the landscape theory pioneered in [FM, ADFJM2, DFM] to the tight-binding Schrödinger operator on $\Z^d$. In particular, we establish upper and lower bounds for the integrated density of states in terms of the counting function based upon the localization landscape. The present paper extends the landscape theory pioneered in [FM, ADFJM2, DFM] to the tight-binding Schrödinger operator on $\Z^d$. In particular, we establish upper and lower bounds for the integrated density of states in terms of the counting function based upon the localization landscape. △ Less

Submitted 8 June, 2021; v1 submitted 8 January, 2021; originally announced January 2021.

Sharp estimates for the integrated density of states in Anderson tight-binding models

Authors: Perceval Desforges, Svitlana Mayboroda, Shiwen Zhang, Guy David, Douglas N. Arnold, Wei Wang, Marcel Filoche

Abstract: Recent work [G. David, M. Filoche, and S. Mayboroda, arXiv:1909.10558[Adv. Math. (to be published)]] has proved the existence of bounds from above and below for the Integrated Density of States (IDOS) of the Schrödinger operator throughout the spectrum, called the landscape law. These bounds involve dimensional constants whose optimal values are yet to be determined. Here, we investigate the accur… ▽ More Recent work [G. David, M. Filoche, and S. Mayboroda, arXiv:1909.10558[Adv. Math. (to be published)]] has proved the existence of bounds from above and below for the Integrated Density of States (IDOS) of the Schrödinger operator throughout the spectrum, called the landscape law. These bounds involve dimensional constants whose optimal values are yet to be determined. Here, we investigate the accuracy of the landscape law in 1D and 2D tight-binding Anderson models, with binary or uniform random distributions. We show, in particular, that in 1D, the IDOS can be approximated with high accuracy through a single formula involving a remarkably simple multiplicative energy shift. In 2D, the same idea applies but the prefactor has to be changed between the bottom and top parts of the spectrum. △ Less

Submitted 18 July, 2021; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: 8 pages, 8 figues

MSC Class: 35J10; 65N30

Journal ref: Phys. Rev. A 104, 012207 (2021)

Moving with the Times: Investigating the Alt-Right Network Gab with Temporal Interaction Graphs

Authors: Naomi A. Arnold, Benjamin A. Steer, Imane Hafnaoui, Hugo A. Parada G., Raul J. Mondragon, Felix Cuadrado, Richard G. Clegg

Abstract: Gab is an online social network often associated with the alt-right political movement and users barred from other networks. It presents an interesting opportunity for research because near-complete data is available from day one of the network's creation. In this paper, we investigate the evolution of the user interaction graph, that is the graph where a link represents a user interacting with an… ▽ More Gab is an online social network often associated with the alt-right political movement and users barred from other networks. It presents an interesting opportunity for research because near-complete data is available from day one of the network's creation. In this paper, we investigate the evolution of the user interaction graph, that is the graph where a link represents a user interacting with another user at a given time. We view this graph both at different times and at different timescales. The latter is achieved by using sliding windows on the graph which gives a novel perspective on social network data. The Gab network is relatively slowly growing over the period of months but subject to large bursts of arrivals over hours and days. We identify plausible events that are of interest to the Gab community associated with the most obvious such bursts. The network is characterised by interactions between `strangers' rather than by reinforcing links between `friends'. Gab usage follows the diurnal cycle of the predominantly US and Europe based users. At off-peak hours the Gab interaction network fragments into sub-networks with absolutely no interaction between them. A small group of users are highly influential across larger timescales, but a substantial number of users gain influence for short periods of time. Temporal analysis at different timescales gives new insights above and beyond what could be found on static graphs. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 10 pages, 9 figures

Complexes from complexes

Authors: Douglas N. Arnold, Kaibo Hu

Abstract: This paper is concerned with the derivation and properties of differential complexes arising from a variety of problems in differential equations, with applications in continuum mechanics, relativity, and other fields. We present a systematic procedure which, starting from well-understood differential complexes such as the de Rham complex, derives new complexes and deduces the properties of the ne… ▽ More This paper is concerned with the derivation and properties of differential complexes arising from a variety of problems in differential equations, with applications in continuum mechanics, relativity, and other fields. We present a systematic procedure which, starting from well-understood differential complexes such as the de Rham complex, derives new complexes and deduces the properties of the new complexes from the old. We relate the cohomology of the output complex to that of the input complexes and show that the new complex has closed ranges, and, consequently, satisfies a Hodge decomposition, Poincaré type inequalities, well-posed Hodge-Laplacian boundary value problems, regular decomposition, and compactness properties on general Lipschitz domains. △ Less

Submitted 25 January, 2021; v1 submitted 25 May, 2020; originally announced May 2020.

Comments: 31 pages. This preprint corresponds to the version accepted by Foundations of Computational Mathematics

MSC Class: 58J10 (Primary) 58A12; 58A14; 35J58 (Secondary)

Search for the double-beta decay of 82Se to the excited states of 82Kr with NEMO-3

Authors: The NEMO-3 collaboration R. Arnold, C. Augier, A. S. Barabash, A. Basharina-Freshville, S. Blondel, S. Blot, M. Bongrand, D. Boursette, R. Breier, V. Brudanin, J. Busto, A. J. Caffrey, S. Calvez, M. Cascella, C. Cerna, J. P. Cesar, A. Chapon, E. Chauveau, A. Chopra, L. Dawson, D. Duchesneau, D. Durand, V. Egorov, G. Eurin, J. J. Evans , et al. (82 additional authors not shown)

Abstract: The double-beta decay of 82Se to the 0+1 excited state of 82Kr has been studied with the NEMO-3 detector using 0.93 kg of enriched 82Se measured for 4.75 y, corresponding to an exposure of 4.42 kg y. A dedicated analysis to reconstruct the gamma-rays has been performed to search for events in the 2e2g channel. No evidence of a 2nbb decay to the 0+1 state has been observed and a limit of T2n 1/2(82… ▽ More The double-beta decay of 82Se to the 0+1 excited state of 82Kr has been studied with the NEMO-3 detector using 0.93 kg of enriched 82Se measured for 4.75 y, corresponding to an exposure of 4.42 kg y. A dedicated analysis to reconstruct the gamma-rays has been performed to search for events in the 2e2g channel. No evidence of a 2nbb decay to the 0+1 state has been observed and a limit of T2n 1/2(82Se; 0+gs -> 0+1) > 1.3 1021 y at 90% CL has been set. Concerning the 0nbb decay to the 0+1 state, a limit for this decay has been obtained with T0n 1/2(82Se; 0+g s -> 0+1) > 2.3 1022 y at 90% CL, independently from the 2nbb decay process. These results are obtained for the first time with a tracko-calo detector, reconstructing every particle in the final state. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Journal ref: Nuclear Physics A Volume 996, April 2020, 121701

The optical tweezer of ferroelectric skyrmions

Authors: X. -G. Wang, L. Chotorlishvili, V. K. Dugaev, A. Ernst, I. Maznichenko, N. Arnold, Chenglong Jia, J. Berakdar, I. Mertig, J. Barnaś

Abstract: Strong magneto-electric coupling in two-dimensional helical materials leads to a peculiar type of topologically protected solutions -- skyrmions. Coupling between the net ferroelectric polarization and magnetization allows control of the magnetic texture with an external electric field. In this work we propose the model of optical tweezer -- a particular configuration of an external electric field… ▽ More Strong magneto-electric coupling in two-dimensional helical materials leads to a peculiar type of topologically protected solutions -- skyrmions. Coupling between the net ferroelectric polarization and magnetization allows control of the magnetic texture with an external electric field. In this work we propose the model of optical tweezer -- a particular configuration of an external electric field and Gaussian laser beam that can trap or release the skyrmions in a highly controlled manner. Functionality of such a tweezer is visualized by micromagnetic simulations and model analysis. △ Less

Submitted 6 January, 2020; originally announced January 2020.

Software Development for APS-Upgrade LLRF

Authors: T. Madden, N. Arnold, T. Berenc, E. Breeding, T. Fors, G. Shen, S. Veseli, Y. Yang, K. Vodopivec, W. Yoder

Abstract: An extensive service-based software system has been developed for real-time data acquisition (DAQ) on numerous subsystems of the Advanced Photon Source (APS) accelerator. The software, called DAQ, is based upon EPICS 3 Area-Detector with added EPICS 4 services, with EPICS 3 + 4 being called EPICS 7. DAQ is comprised of several network services including file storage and processing as well as numer… ▽ More An extensive service-based software system has been developed for real-time data acquisition (DAQ) on numerous subsystems of the Advanced Photon Source (APS) accelerator. The software, called DAQ, is based upon EPICS 3 Area-Detector with added EPICS 4 services, with EPICS 3 + 4 being called EPICS 7. DAQ is comprised of several network services including file storage and processing as well as numerous EPICS 7 Input-Output Controller (IOC) applications that collect data and publish to EPICS 7 from accelerator subsystems, including power supplies, diagnostics, timing, and RF systems. For the APS upgrade extensions to the DAQ software have been developed for digital Low Level RF (LLRF) systems supporting multiple hardware platforms including Lawrence Berkeley Laboratory's LLRF4 hardware and Micro Telecommunications Computing Architecture. EPICS 7 IOCs for LLRF run on multiple threads for real-time data collection, visualization , system control, and data file storage that can be triggered on accelerator beam dumps for system debug. A Python-based oscilloscope application has been developed for real-time data visualization. △ Less

Submitted 17 October, 2019; v1 submitted 10 October, 2019; originally announced October 2019.

Comments: Talk/Poster presented at LLRF Workshop 2019 (LLRF2019, arXiv:1909.06754)

Report number: LLRF2019/213

Likelihood-based approach to discriminate mixtures of network models that vary in time

Authors: Naomi A. Arnold, Raul J. Mondragon, Richard G. Clegg

Abstract: Discriminating between competing explanatory models as to which is more likely responsible for the growth of a network is a problem of fundamental importance for network science. The rules governing this growth are attributed to mechanisms such as preferential attachment and triangle closure, with a wealth of explanatory models based on these. These models are deliberately simple, commonly with th… ▽ More Discriminating between competing explanatory models as to which is more likely responsible for the growth of a network is a problem of fundamental importance for network science. The rules governing this growth are attributed to mechanisms such as preferential attachment and triangle closure, with a wealth of explanatory models based on these. These models are deliberately simple, commonly with the network growing according to a constant mechanism for its lifetime, to allow for analytical results. We use a likelihood-based framework on artificial data where the network model changes at a known point in time and demonstrate that we can recover the change point from analysis of the network. We then use real datasets and demonstrate how our framework can show the changing importance of network growth mechanisms over time. △ Less

Submitted 8 February, 2021; v1 submitted 29 September, 2019; originally announced September 2019.

Journal ref: Sci Rep 11, 5205 (2021)

The Hellan-Herrmann-Johnson method with curved elements

Authors: Douglas N. Arnold, Shawn W. Walker

Abstract: We study the finite element approximation of the Kirchhoff plate equation on domains with curved boundaries using the Hellan-Herrmann-Johnson (HHJ) method. We prove optimal convergence on domains with piecewise $C^{k+1}$ boundary for $k \geq 1$ when using a parametric (curved) HHJ space. Computational results are given that demonstrate optimal convergence and how convergence degrades when curved t… ▽ More We study the finite element approximation of the Kirchhoff plate equation on domains with curved boundaries using the Hellan-Herrmann-Johnson (HHJ) method. We prove optimal convergence on domains with piecewise $C^{k+1}$ boundary for $k \geq 1$ when using a parametric (curved) HHJ space. Computational results are given that demonstrate optimal convergence and how convergence degrades when curved triangles of insufficient polynomial degree are used. Moreover, we show that the lowest order HHJ method on a polygonal approximation of the disk does not succumb to the classic Babuška paradox, highlighting the geometrically non-conforming aspect of the HHJ method. △ Less

Submitted 25 July, 2020; v1 submitted 20 September, 2019; originally announced September 2019.

Comments: 27 pages, 4 figures, 4 tables. To appear in SIAM Journal on Numerical Analysis

MSC Class: 65N30; 35J40; 35Q72

Localization of eigenfunctions via an effective potential

Authors: Douglas N. Arnold, Guy David, Marcel Filoche, David Jerison, Svitlana Mayboroda

Abstract: We consider the localization of eigenfunctions for the operator $L=-\mbox{div} A \nabla + V$ on a Lipschitz domain $Ωおめが$ and, more generally, on manifolds with and without boundary. In earlier work, two authors of the present paper demonstrated the remarkable ability of the landscape, defined as the solution to $Lu=1$, to predict the location of the localized eigenfunctions. Here, we explain and jus… ▽ More We consider the localization of eigenfunctions for the operator $L=-\mbox{div} A \nabla + V$ on a Lipschitz domain $Ωおめが$ and, more generally, on manifolds with and without boundary. In earlier work, two authors of the present paper demonstrated the remarkable ability of the landscape, defined as the solution to $Lu=1$, to predict the location of the localized eigenfunctions. Here, we explain and justify a new framework that reveals a richly detailed portrait of the eigenfunctions and eigenvalues. We show that the reciprocal of the landscape function, $1/u$, acts as an effective potential. Hence from the single measurement of $u$, we obtain, via $1/u$, explicit bounds on the exponential decay of the eigenfunctions of the system and estimates on the distribution of eigenvalues near the bottom of the spectrum. (This version strengthens and simplifies the results of the first one by replacing a global bi-Lipschitz hypothesis on the domain with a local bi-Lipschitz hypothesis. It improves on the second version by adding pictures and numerical examples. This version is identical to the third version; all that is changed is to correct some tex mistakes in symbols in this abstract. There are no changes to the paper itself.) △ Less

Submitted 3 October, 2018; v1 submitted 6 December, 2017; originally announced December 2017.

Comments: 33 pages

Report number: Journal Communications in Partial Differential Equations 44 (2019) 1186-1216 MSC Class: 35P; 35Q40

Correlated In-Situ Low-Frequency Noise and Impedance Spectroscopy Reveal Recombination Dynamics in Organic Solar Cells using Fullerene and Non-Fullerene Acceptors

Authors: Kyle A. Luck, Vinod K. Sangwan, Patrick E. Hartnett, Heather N. Arnold, Michael R. Wasielewski, Tobin J. Marks, Mark C. Hersam

Abstract: Non-fullerene acceptors based on perylenediimides (PDIs) have garnered significant interest as an alternative to fullerene acceptors in organic photovoltaics (OPVs), but their charge transport phenomena are not well understood, especially in bulk heterojunctions (BHJs). Here, we investigate charge transport and current fluctuations by performing correlated low-frequency noise and impedance spectro… ▽ More Non-fullerene acceptors based on perylenediimides (PDIs) have garnered significant interest as an alternative to fullerene acceptors in organic photovoltaics (OPVs), but their charge transport phenomena are not well understood, especially in bulk heterojunctions (BHJs). Here, we investigate charge transport and current fluctuations by performing correlated low-frequency noise and impedance spectroscopy measurements on two BHJ OPV systems, one employing a fullerene acceptor and the other employing a dimeric PDI acceptor. In the dark, these measurements reveal that PDI-based OPVs have a greater degree of recombination in comparison to fullerene-based OPVs. Furthermore, for the first time in organic solar cells, 1/f noise data are fit to the Kleinpenning model to reveal underlying current fluctuations in different transport regimes. Under illumination, 1/f noise increases by approximately four orders of magnitude for the fullerene-based OPVs and three orders of magnitude for the PDI-based OPVs. An inverse correlation is also observed between noise spectral density and power conversion efficiency. Overall, these results show that low-frequency noise spectroscopy is an effective in-situ diagnostic tool to assess charge transport in emerging photovoltaic materials, thereby providing quantitative guidance for the design of next-generation solar cell materials and technologies. △ Less

Submitted 20 September, 2017; originally announced September 2017.

Comments: 37 pages, 7 figures

Finite element exterior calculus with lower-order terms

Authors: Douglas N. Arnold, Lizao Li

Abstract: The scalar and vector Laplacians are basic operators in physics and engineering. In applications, they show up frequently perturbed by lower-order terms. The effect of such perturbations on mixed finite element methods in the scalar case is well-understood, but that in the vector case is not. In this paper, we first show that surprisingly for certain elements there is degradation of the convergenc… ▽ More The scalar and vector Laplacians are basic operators in physics and engineering. In applications, they show up frequently perturbed by lower-order terms. The effect of such perturbations on mixed finite element methods in the scalar case is well-understood, but that in the vector case is not. In this paper, we first show that surprisingly for certain elements there is degradation of the convergence rates with certain lower-order terms even when both the solution and the data are smooth. We then give a systematic analysis of lower-order terms in mixed methods by extending the Finite Element Exterior Calculus (FEEC) framework, which contains the scalar, vector Laplacian, and many other elliptic operators as special cases. We prove that stable mixed discretization remains stable with lower-order terms for sufficiently fine discretization. Moreover, we derive sharp improved error estimates for each individual variable. In particular, this yields new results for the vector Laplacian problem which are useful in applications such as electromagnetism and acoustics modeling. Further our results imply many previous results for the scalar problem and thus unifies them all under the FEEC framework. △ Less

Submitted 25 February, 2016; v1 submitted 22 September, 2015; originally announced September 2015.

MSC Class: 65N30 (Primary)

Journal ref: Math. Comp. 86 (2017) 2193-2212

Effective confining potential of quantum states in disordered media

Authors: Douglas N. Arnold, Guy David, David Jerison, Svitlana Mayboroda, Marcel Filoche

Abstract: The amplitude of localized quantum states in random or disordered media may exhibit long range exponential decay. We present here a theory that unveils the existence of an effective potential which finely governs the confinement of these states. In this picture, the boundaries of the localization subregions for low energy eigenfunctions correspond to the barriers of this effective potential, and t… ▽ More The amplitude of localized quantum states in random or disordered media may exhibit long range exponential decay. We present here a theory that unveils the existence of an effective potential which finely governs the confinement of these states. In this picture, the boundaries of the localization subregions for low energy eigenfunctions correspond to the barriers of this effective potential, and the long range exponential decay characteristic of Anderson localization is explained as the consequence of multiple tunneling in the dense network of barriers created by this effective potential. Finally, we show that the Weyl's formula based on this potential turns out to be a remarkable approximation of the density of states for a large variety of one-dimensional systems, periodic or random. △ Less

Submitted 2 February, 2016; v1 submitted 11 May, 2015; originally announced May 2015.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. Lett. 116, 056602 (2016)

Mixed methods for elastodynamics with weak symmetry

Authors: Douglas N. Arnold, Jeonghun J. Lee

Abstract: We analyze the application to elastodynamic problems of mixed finite element methods for elasticity with weak symmetry. Our approach leads to a semidiscrete method which consists of a system of ordinary differential equations without algebraic constraints. Our error analysis, which is based on a new elliptic projection operator, applies to several mixed finite element spaces developed for elastost… ▽ More We analyze the application to elastodynamic problems of mixed finite element methods for elasticity with weak symmetry. Our approach leads to a semidiscrete method which consists of a system of ordinary differential equations without algebraic constraints. Our error analysis, which is based on a new elliptic projection operator, applies to several mixed finite element spaces developed for elastostatics. The error estimates we obtain are robust for nearly incompressible materials. △ Less

Submitted 21 December, 2013; originally announced December 2013.

Comments: 20 pages

MSC Class: 65N30; 74H15; 74S05

Journal ref: SIAM J. Numer. Anal. 52 (2014) 2743-2769

Low Frequency Electronic Noise in Single-Layer MoS2 Transistors

Authors: Vinod K. Sangwan, Heather N. Arnold, Deep Jariwala, Tobin J. Marks, Lincoln J. Lauhon, Mark C. Hersam

Abstract: Ubiquitous low frequency 1/f noise can be a limiting factor in the performance and application of nanoscale devices. Here, we quantitatively investigate low frequency electronic noise in single-layer transition metal dichalcogenide MoS2 field-effect transistors. The measured 1/f noise can be explained by an empirical formulation of mobility fluctuations with the Hooge parameter ranging between 0.0… ▽ More Ubiquitous low frequency 1/f noise can be a limiting factor in the performance and application of nanoscale devices. Here, we quantitatively investigate low frequency electronic noise in single-layer transition metal dichalcogenide MoS2 field-effect transistors. The measured 1/f noise can be explained by an empirical formulation of mobility fluctuations with the Hooge parameter ranging between 0.005 and 2.0 in vacuum (< 10-5 Torr). The field-effect mobility decreased and the noise amplitude increased by an order of magnitude in ambient conditions, revealing the significant influence of atmospheric adsorbates on charge transport. In addition, single Lorentzian generation-recombination noise was observed to increase by an order of magnitude as the devices were cooled from 300 K to 6.5 K. △ Less

Submitted 15 August, 2013; originally announced August 2013.

Comments: Nano Letters (2013)

Mixed finite elements for elasticity on quadrilateral meshes

Authors: Douglas N. Arnold, Gerard Awanou, Weifeng Qiu

Abstract: We present stable mixed finite elements for planar linear elasticity on general quadrilateral meshes. The symmetry of the stress tensor is imposed weakly and so there are three primary variables, the stress tensor, the displacement vector field, and the scalar rotation. We develop and analyze a stable family of methods, indexed by an integer $r \geq 2$ and with rate of convergence in the $L^2$ nor… ▽ More We present stable mixed finite elements for planar linear elasticity on general quadrilateral meshes. The symmetry of the stress tensor is imposed weakly and so there are three primary variables, the stress tensor, the displacement vector field, and the scalar rotation. We develop and analyze a stable family of methods, indexed by an integer $r \geq 2$ and with rate of convergence in the $L^2$ norm of order $r$ for all the variables. The methods use Raviart-Thomas elements for the stress, piecewise tensor product polynomials for the displacement, and piecewise polynomials for the rotation. We also present a simple first order element, not belonging to this family. It uses the lowest order BDM elements for the stress, and piecewise constants for the displacement and rotation, and achieves first order convergence for all three variables. △ Less

Submitted 19 April, 2014; v1 submitted 28 June, 2013; originally announced June 2013.

Comments: v2 is the revision made in response to the referees

MSC Class: 65N30

Extrinsic and Intrinsic Photoresponse in Monodisperse Carbon Nanotube Thin Film Transistors

Authors: Erik Sczygelski, Vinod K. Sangwan, Chung-Chiang Wu, Heather N. Arnold, Ken Everaerts, Tobin J. Marks, Mark C. Hersam, Lincoln J. Lauhon

Abstract: Spectroscopic, time-resolved scanning photocurrent microscopy is shown to distinguish the intrinsic photoresponse of monodisperse semiconducting (99%) single-walled carbon nanotubes (SWCNTs) from the extrinsic photoresponse of the substrate. A persistent positive photocurrent induced by near-IR excitation is attributed to the generation of free carriers by inter-band excitation in SWCNTs. For shor… ▽ More Spectroscopic, time-resolved scanning photocurrent microscopy is shown to distinguish the intrinsic photoresponse of monodisperse semiconducting (99%) single-walled carbon nanotubes (SWCNTs) from the extrinsic photoresponse of the substrate. A persistent positive photocurrent induced by near-IR excitation is attributed to the generation of free carriers by inter-band excitation in SWCNTs. For shorter excitation wavelengths, absorption by the Si substrate generates two types of photocurrent: a transient positive photoresponse, identified as a displacement current, and a persistent negative photocurrent that arises from photogating of the SWCNT thin film. △ Less

Submitted 9 February, 2013; originally announced February 2013.

Finite element differential forms on curvilinear cubic meshes and their approximation properties

Authors: Douglas N. Arnold, Daniele Boffi, Francesca Bonizzoni

Abstract: We study the approximation properties of a wide class of finite element differential forms on curvilinear cubic meshes in n dimensions. Specifically, we consider meshes in which each element is the image of a cubical reference element under a diffeomorphism, and finite element spaces in which the shape functions and degrees of freedom are obtained from the reference element by pullback of differen… ▽ More We study the approximation properties of a wide class of finite element differential forms on curvilinear cubic meshes in n dimensions. Specifically, we consider meshes in which each element is the image of a cubical reference element under a diffeomorphism, and finite element spaces in which the shape functions and degrees of freedom are obtained from the reference element by pullback of differential forms. In the case where the diffeomorphisms from the reference element are all affine, i.e., mesh consists of parallelotopes, it is standard that the rate of convergence in L2 exceeds by one the degree of the largest full polynomial space contained in the reference space of shape functions. When the diffeomorphism is multilinear, the rate of convergence for the same space of reference shape function may degrade severely, the more so when the form degree is larger. The main result of the paper gives a sufficient condition on the reference shape functions to obtain a given rate of convergence. △ Less

Submitted 30 April, 2014; v1 submitted 28 December, 2012; originally announced December 2012.

Comments: 17 pages, 1 figure; v2: changes in response to referee reports; v3: minor additional changes, this version accepted for Numerische Mathematik; v3: very minor updates, this version corresponds to the final published version

MSC Class: 65N30; 41A25; 41A10; 41A15; 41A63

Journal ref: Numer. Math. 129 (2015) 1-20

On the consistency of the combinatorial codifferential

Authors: Douglas N. Arnold, Richard S. Falk, Johnny Guzmán, Gantumur Tsogtgerel

Abstract: In 1976, Dodziuk and Patodi employed Whitney forms to define a combinatorial codifferential operator on cochains, and they raised the question whether it is consistent in the sense that for a smooth enough differential form the combinatorial codifferential of the associated cochain converges to the exterior codifferential of the form as the triangulation is refined. In 1991, Smits proved this to b… ▽ More In 1976, Dodziuk and Patodi employed Whitney forms to define a combinatorial codifferential operator on cochains, and they raised the question whether it is consistent in the sense that for a smooth enough differential form the combinatorial codifferential of the associated cochain converges to the exterior codifferential of the form as the triangulation is refined. In 1991, Smits proved this to be the case for the combinatorial codifferential applied to 1-forms in two dimensions under the additional assumption that the initial triangulation is refined in a completely regular fashion, by dividing each triangle into four similar triangles. In this paper we extend Smits's result to arbitrary dimensions, showing that the combinatorial codifferential on 1-forms is consistent if the triangulations are uniform or piecewise uniform in a certain precise sense. We also show that this restriction on the triangulations is needed, giving a counterexample in which a different regular refinement procedure, namely Whitney's standard subdivision, is used. Further, we show by numerical example that for 2-forms in three dimensions, the combinatorial codifferential is not consistent even for the most regular subdivision process. △ Less

Submitted 4 February, 2014; v1 submitted 18 December, 2012; originally announced December 2012.

Comments: 16 pages, 7 figures, 5 tables; version 2 has minor stylistic corrections, corresponds to version to appear in Transactions of the AMS 2014

MSC Class: Primary 58A10; 65N30; Secondary 39A12; 57Q55; 58A14

Journal ref: Trans. Amer. Math. Soc. 366 (2014) 5487-5502

Nonconforming tetrahedral mixed finite elements for elasticity

Authors: Douglas N. Arnold, Gerard Awanou, Ragnar Winther

Abstract: This paper presents a nonconforming finite element approximation of the space of symmetric tensors with square integrable divergence, on tetrahedral meshes. Used for stress approximation together with the full space of piecewise linear vector fields for displacement, this gives a stable mixed finite element method which is shown to be linearly convergent for both the stress and displacement, and w… ▽ More This paper presents a nonconforming finite element approximation of the space of symmetric tensors with square integrable divergence, on tetrahedral meshes. Used for stress approximation together with the full space of piecewise linear vector fields for displacement, this gives a stable mixed finite element method which is shown to be linearly convergent for both the stress and displacement, and which is significantly simpler than any stable conforming mixed finite element method. The method may be viewed as the three-dimensional analogue of a previously developed element in two dimensions. As in that case, a variant of the method is proposed as well, in which the displacement approximation is reduced to piecewise rigid motions and the stress space is reduced accordingly, but the linear convergence is retained. △ Less

Submitted 21 May, 2013; v1 submitted 23 October, 2012; originally announced October 2012.

Comments: 13 pages, 2 figures

MSC Class: 65N30; secondary: 74S05

Journal ref: Math. Models Methods Appl. Sci., 24:1-14, 2014

Finite element exterior calculus for parabolic problems

Authors: Douglas N. Arnold, Hongtao Chen

Abstract: In this paper, we consider the extension of the finite element exterior calculus from elliptic problems, in which the Hodge Laplacian is an appropriate model problem, to parabolic problems, for which we take the Hodge heat equation as our model problem. The numerical method we study is a Galerkin method based on a mixed variational formulation and using as subspaces the same spaces of finite eleme… ▽ More In this paper, we consider the extension of the finite element exterior calculus from elliptic problems, in which the Hodge Laplacian is an appropriate model problem, to parabolic problems, for which we take the Hodge heat equation as our model problem. The numerical method we study is a Galerkin method based on a mixed variational formulation and using as subspaces the same spaces of finite element differential forms which are used for elliptic problems. We analyze both the semidiscrete and a fully-discrete numerical scheme. △ Less

Submitted 5 September, 2012; originally announced September 2012.

Comments: 17 pages

MSC Class: 65N30

Journal ref: ESAIM Math. Model. Numer. Anal. 51 (2017) 17-34

Spaces of finite element differential forms

Authors: Douglas N. Arnold

Abstract: We discuss the construction of finite element spaces of differential forms which satisfy the crucial assumptions of the finite element exterior calculus, namely that they can be assembled into subcomplexes of the de Rham complex which admit commuting projections. We present two families of spaces in the case of simplicial meshes, and two other families in the case of cubical meshes. We make use of… ▽ More We discuss the construction of finite element spaces of differential forms which satisfy the crucial assumptions of the finite element exterior calculus, namely that they can be assembled into subcomplexes of the de Rham complex which admit commuting projections. We present two families of spaces in the case of simplicial meshes, and two other families in the case of cubical meshes. We make use of the exterior calculus and the Koszul complex to define and understand the spaces. These tools allow us to treat a wide variety of situations, which are often treated separately, in a unified fashion. △ Less

Submitted 27 November, 2012; v1 submitted 9 August, 2012; originally announced August 2012.

Comments: To appear in: Analysis and Numerics of Partial Differential Equations, U. Gianazza, F. Brezzi, P. Colli Franzone, and G. Gilardi, eds., Springer 2013. v2: a few minor typos corrected. v3: a few more typo corrections

MSC Class: 65N30

Journal ref: U. Gianazza, F. Brezzi, P. Colli Franzone, and G. Gilardi, editors, Analysis and Numerics of Partial Differential Equations, pages 117-140. Springer, 2013

Finite element differential forms on cubical meshes

Authors: Douglas N. Arnold, Gerard Awanou

Abstract: We develop a family of finite element spaces of differential forms defined on cubical meshes in any number of dimensions. The family contains elements of all polynomial degrees and all form degrees. In two dimensions, these include the serendipity finite elements and the rectangular BDM elements. In three dimensions they include a recent generalization of the serendipity spaces, and new H(curl) an… ▽ More We develop a family of finite element spaces of differential forms defined on cubical meshes in any number of dimensions. The family contains elements of all polynomial degrees and all form degrees. In two dimensions, these include the serendipity finite elements and the rectangular BDM elements. In three dimensions they include a recent generalization of the serendipity spaces, and new H(curl) and H(div) finite element spaces. Spaces in the family can be combined to give finite element subcomplexes of the de Rham complex which satisfy the basic hypotheses of the finite element exterior calculus, and hence can be used for stable discretization of a variety of problems. The construction and properties of the spaces are established in a uniform manner using finite element exterior calculus. △ Less

Submitted 29 July, 2013; v1 submitted 11 April, 2012; originally announced April 2012.

Comments: v2: as accepted by Mathematics of Computation after minor revisions; v3: this version corresponds to the final version for Math. Comp., after copyediting and galley proofs

MSC Class: 65N30

Journal ref: Math. Comp. 83 (2014) 1551-1570

Mathematicians take a stand

Authors: Douglas N. Arnold, Henry Cohn

Abstract: We survey the reasons for the ongoing boycott of the publisher Elsevier. We examine Elsevier's pricing and bundling policies, restrictions on dissemination by authors, and lapses in ethics and peer review, and we conclude with thoughts about the future of mathematical publishing. We survey the reasons for the ongoing boycott of the publisher Elsevier. We examine Elsevier's pricing and bundling policies, restrictions on dissemination by authors, and lapses in ethics and peer review, and we conclude with thoughts about the future of mathematical publishing. △ Less

Submitted 20 June, 2012; v1 submitted 5 April, 2012; originally announced April 2012.

Comments: 5 pages

Journal ref: Notices of the American Mathematical Society 59 (2012), 828-833

Mixed finite element approximation of the vector Laplacian with Dirichlet boundary conditions

Authors: Douglas N. Arnold, Richard S. Falk, Jay Gopalakrishnan

Abstract: We consider the finite element solution of the vector Laplace equation on a domain in two dimensions. For various choices of boundary conditions, it is known that a mixed finite element method, in which the rotation of the solution is introduced as a second unknown, is advantageous, and appropriate choices of mixed finite element spaces lead to a stable, optimally convergent discretization. Howeve… ▽ More We consider the finite element solution of the vector Laplace equation on a domain in two dimensions. For various choices of boundary conditions, it is known that a mixed finite element method, in which the rotation of the solution is introduced as a second unknown, is advantageous, and appropriate choices of mixed finite element spaces lead to a stable, optimally convergent discretization. However, the theory that leads to these conclusions does not apply to the case of Dirichlet boundary conditions, in which both components of the solution vanish on the boundary. We show, by computational example, that indeed such mixed finite elements do not perform optimally in this case, and we analyze the suboptimal convergence that does occur. As we indicate, these results have implications for the solution of the biharmonic equation and of the Stokes equations using a mixed formulation involving the vorticity. △ Less

Submitted 16 September, 2011; originally announced September 2011.

MSC Class: 65N30

Journal ref: Math. Models Methods Appl. Sci., 22(9):26 pages, 2012

The serendipity family of finite elements

Authors: Douglas N. Arnold, Gerard Awanou

Abstract: We give a new, simple, dimension-independent definition of the serendipity finite element family. The shape functions are the span of all monomials which are linear in at least s-r of the variables where s is the degree of the monomial or, equivalently, whose superlinear degree (total degree with respect to variables entering at least quadratically) is at most r. The degrees of freedom are given b… ▽ More We give a new, simple, dimension-independent definition of the serendipity finite element family. The shape functions are the span of all monomials which are linear in at least s-r of the variables where s is the degree of the monomial or, equivalently, whose superlinear degree (total degree with respect to variables entering at least quadratically) is at most r. The degrees of freedom are given by moments of degree at most r-2d on each face of dimension d. We establish unisolvence and a geometric decomposition of the space. △ Less

Submitted 3 January, 2011; originally announced January 2011.

Comments: 7 pages, 1 figure, 1 table. To appear in Foundations in Computational Mathematics, 2011

MSC Class: 65N30

Journal ref: Found. Comput. Math., 11:337-344, 2011

Nefarious Numbers

Authors: Douglas N. Arnold, Kristine K. Fowler

Abstract: We investigate the journal impact factor, focusing on the applied mathematics category. We discuss impact factor manipulation and demonstrate that the impact factor gives an inaccurate view of journal quality, which is poorly correlated with expert opinion. We investigate the journal impact factor, focusing on the applied mathematics category. We discuss impact factor manipulation and demonstrate that the impact factor gives an inaccurate view of journal quality, which is poorly correlated with expert opinion. △ Less

Submitted 17 November, 2010; v1 submitted 1 October, 2010; originally announced October 2010.

Comments: 5 pages, 1 figure

Journal ref: Notices Amer. Math. Soc., 58(3):434-437, 2011

Stability of Lagrange elements for the mixed Laplacian

Authors: Douglas N. Arnold, Marie E. Rognes

Abstract: The stability properties of simple element choices for the mixed formulation of the Laplacian are investigated numerically. The element choices studied use vector Lagrange elements, i.e., the space of continuous piecewise polynomial vector fields of degree at most r, for the vector variable, and the divergence of this space, which consists of discontinuous piecewise polynomials of one degree low… ▽ More The stability properties of simple element choices for the mixed formulation of the Laplacian are investigated numerically. The element choices studied use vector Lagrange elements, i.e., the space of continuous piecewise polynomial vector fields of degree at most r, for the vector variable, and the divergence of this space, which consists of discontinuous piecewise polynomials of one degree lower, for the scalar variable. For polynomial degrees r equal 2 or 3, this pair of spaces was found to be stable for all mesh families tested. In particular, it is stable on diagonal mesh families, in contrast to its behaviour for the Stokes equations. For degree r equal 1, stability holds for some meshes, but not for others. Additionally, convergence was observed precisely for the methods that were observed to be stable. However, it seems that optimal order L2 estimates for the vector variable, known to hold for r>3, do not hold for lower degrees. △ Less

Submitted 20 July, 2009; originally announced July 2009.

Comments: to appear in Calcolo

MSC Class: 65N30

Journal ref: Calcolo 46 (2009) 245-260

Finite element exterior calculus: from Hodge theory to numerical stability

Authors: Douglas N. Arnold, Richard S. Falk, Ragnar Winther

Abstract: This article reports on the confluence of two streams of research, one emanating from the fields of numerical analysis and scientific computation, the other from topology and geometry. In it we consider the numerical discretization of partial differential equations that are related to differential complexes so that de Rham cohomology and Hodge theory are key tools for the continuous problem. Aft… ▽ More This article reports on the confluence of two streams of research, one emanating from the fields of numerical analysis and scientific computation, the other from topology and geometry. In it we consider the numerical discretization of partial differential equations that are related to differential complexes so that de Rham cohomology and Hodge theory are key tools for the continuous problem. After a brief introduction to finite element methods, the discretization methods we consider, we develop an abstract Hilbert space framework for analyzing stability and convergence. In this framework, the differential complex is represented by a complex of Hilbert spaces and stability is obtained by transferring Hodge theoretic structures from the continuous level to the discrete. We show stable discretization is achieved if the finite element spaces satisfy two hypotheses: they form a subcomplex and there exists a bounded cochain projection from the full complex to the subcomplex. Next, we consider the most canonical example of the abstract theory, in which the Hilbert complex is the de Rham complex of a domain in Euclidean space. We use the Koszul complex to construct two families of finite element differential forms, show that these can be arranged in subcomplexes of the de Rham complex in numerous ways, and for each construct a bounded cochain projection. The abstract theory therefore applies to give the stability and convergence of finite element approximations of the Hodge Laplacian. Other applications are considered as well, especially to the equations of elasticity. Background material is included to make the presentation self-contained for a variety of readers. △ Less

Submitted 26 August, 2009; v1 submitted 23 June, 2009; originally announced June 2009.

Comments: 74 pages, 8 figures; reorganized introductory material, added additional example and references; final version accepted by Bulletin of the AMS, added references to codes and adjusted some diagrams. Bulletin of the American Mathematical Society, to appear 2010

MSC Class: 65N30; 58A14

Journal ref: Bull. Amer. Math. Soc. (N.S.), 47:281-354, 2010

Geometric decompositions and local bases for spaces of finite element differential forms

Authors: Douglas N. Arnold, Richard S. Falk, Ragnar Winther

Abstract: We study the two primary families of spaces of finite element differential forms with respect to a simplicial mesh in any number of space dimensions. These spaces are generalizations of the classical finite element spaces for vector fields, frequently referred to as Raviart-Thomas, Brezzi-Douglas-Marini, and Nedelec spaces. In the present paper, we derive geometric decompositions of these spaces… ▽ More We study the two primary families of spaces of finite element differential forms with respect to a simplicial mesh in any number of space dimensions. These spaces are generalizations of the classical finite element spaces for vector fields, frequently referred to as Raviart-Thomas, Brezzi-Douglas-Marini, and Nedelec spaces. In the present paper, we derive geometric decompositions of these spaces which lead directly to explicit local bases for them, generalizing the Bernstein basis for ordinary Lagrange finite elements. The approach applies to both families of finite element spaces, for arbitrary polynomial degree, arbitrary order of the differential forms, and an arbitrary simplicial triangulation in any number of space dimensions. A prominent role in the construction is played by the notion of a consistent family of extension operators, which expresses in an abstract framework a sufficient condition for deriving a geometric decomposition of a finite element space leading to a local basis. △ Less

Submitted 6 June, 2008; originally announced June 2008.

Comments: 25 pages, 2 figures, submitted to Computer Methods in Applied Mechanics and Engineering

MSC Class: 65N30

Journal ref: Comput. Methods Appl. Mech. Engrg., 198:1660-1672, 2009

Habitat of early life: Solar X-ray and UV radiation at Earth's surface 4-3.5 billion years ago

Authors: I. Cnossen, J. Sanz-Forcada, F. Favata, O. Witasse, T. Zegers, N. F. Arnold

Abstract: Solar X-ray and UV radiation (0.1-320 nm) received at Earth's surface is an important aspect of the circumstances under which life formed on Earth. The quantity that is received depends on two main variables: the emission of radiation by the young Sun and its extinction through absorption and scattering by the Earth's early atmosphere. The spectrum emitted by the Sun when life formed, between 4… ▽ More Solar X-ray and UV radiation (0.1-320 nm) received at Earth's surface is an important aspect of the circumstances under which life formed on Earth. The quantity that is received depends on two main variables: the emission of radiation by the young Sun and its extinction through absorption and scattering by the Earth's early atmosphere. The spectrum emitted by the Sun when life formed, between 4 and 3.5 Ga, was modeled here, including the effects of flares and activity cycles, using a solar-like star that has the same age now as the Sun had 4-3.5 Ga. Atmospheric extinction was calculated using the Beer-Lambert law, assuming several density profiles for the atmosphere of the Archean Earth. We found that almost all radiation with a wavelength shorter than 200 nm is attenuated effectively, even by very tenuous atmospheres. Longer-wavelength radiation is progressively less well attenuated, and its extinction is more sensitive to atmospheric composition. Minor atmospheric components, such as methane, ozone, water vapor, etc., have only negligible effects, but changes in CO2 concentration can cause large differences in surface flux. Differences due to variability in solar emission are small compared to this. In all cases surface radiation levels on the Archean Earth were several orders of magnitude higher in the 200-300 nm wavelength range than current levels in this range. That means that any form of life that might have been present at Earth's surface 4-3.5 Ga must have been exposed to much higher quantities of damaging radiation than at present. △ Less

Submitted 20 February, 2007; originally announced February 2007.

Comments: 19 pages, published by JGR-Planets

Journal ref: J.Geophys.Res.112:E02008,2007

Mixed finite element methods for linear elasticity with weakly imposed symmetry

Authors: Douglas N. Arnold, Richard S. Falk, Ragnar Winther

Abstract: In this paper, we construct new finite element methods for the approximation of the equations of linear elasticity in three space dimensions that produce direct approximations to both stresses and displacements. The methods are based on a modified form of the Hellinger--Reissner variational principle that only weakly imposes the symmetry condition on the stresses. Although this approach has been… ▽ More In this paper, we construct new finite element methods for the approximation of the equations of linear elasticity in three space dimensions that produce direct approximations to both stresses and displacements. The methods are based on a modified form of the Hellinger--Reissner variational principle that only weakly imposes the symmetry condition on the stresses. Although this approach has been previously used by a number of authors, a key new ingredient here is a constructive derivation of the elasticity complex starting from the de Rham complex. By mimicking this construction in the discrete case, we derive new mixed finite elements for elasticity in a systematic manner from known discretizations of the de Rham complex. These elements appear to be simpler than the ones previously derived. For example, we construct stable discretizations which use only piecewise linear elements to approximate the stress field and piecewise constant functions to approximate the displacement field. △ Less

Submitted 18 January, 2007; originally announced January 2007.

Comments: to appear in Mathematics of Computation

MSC Class: 65N30; 74S05

Journal ref: Math. Comput., 76:1699-1723, 2007

Finite elements for symmetric tensors in three dimension

Authors: Douglas N. Arnold, Gerard Awanou, Ragnar Winther

Abstract: We construct finite element subspaces of the space of symmetric tensors with square-integrable divergence on a three-dimensional domain. These spaces can be used to approximate the stress field in the classical Hellinger--Reissner mixed formulation of the elasticty equations, when standard discontinous finite element spaces are used to approximate the displacement field. These finite element spa… ▽ More We construct finite element subspaces of the space of symmetric tensors with square-integrable divergence on a three-dimensional domain. These spaces can be used to approximate the stress field in the classical Hellinger--Reissner mixed formulation of the elasticty equations, when standard discontinous finite element spaces are used to approximate the displacement field. These finite element spaces are defined with respect to an arbitrary simplicial triangulation of the domain, and there is one for each positive value of the polynomial degree used for the displacements. For each degree, these provide a stable finite element discretization. The construction of the spaces is closely tied to discretizations of the elasticity complex, and can be viewed as the three-dimensional analogue of the triangular element family for plane elasticity previously proposed by Arnold and Winther. △ Less

Submitted 15 May, 2007; v1 submitted 18 January, 2007; originally announced January 2007.

MSC Class: 65N30; 74S05

Journal ref: Math. Comput., 77:1229-1251, 2008

Boundary conditions for the Einstein-Christoffel formulation of Einstein's equations

Authors: Douglas N. Arnold, Nicolae Tarfulea

Abstract: Specifying boundary conditions continues to be a challenge in numerical relativity in order to obtain a long time convergent numerical simulation of Einstein's equations in domains with artificial boundaries. In this paper, we address this problem for the Einstein--Christoffel (EC) symmetric hyperbolic formulation of Einstein's equations linearized around flat spacetime. First, we prescribe simp… ▽ More Specifying boundary conditions continues to be a challenge in numerical relativity in order to obtain a long time convergent numerical simulation of Einstein's equations in domains with artificial boundaries. In this paper, we address this problem for the Einstein--Christoffel (EC) symmetric hyperbolic formulation of Einstein's equations linearized around flat spacetime. First, we prescribe simple boundary conditions that make the problem well posed and preserve the constraints. Next, we indicate boundary conditions for a system that extends the linearized EC system by including the momentum constraints and whose solution solves Einstein's equations in a bounded domain. △ Less

Submitted 1 November, 2006; originally announced November 2006.

Axially symmetric focusing as a cuspoid diffraction catastrophe: Scalar and vector cases and comparison with the theory of Mie

Authors: Johannes Kofler, Nikita Arnold

Abstract: An analytical description of arbitrary strongly aberrated axially symmetric focusing is developed. This is done by matching the solution of geometrical optics with a wave pattern which is universal for the underlying ray structure. The corresponding canonical integral is the Bessoid integral, which is a three dimensional generalization of the Pearcey integral that approximates the field near an… ▽ More An analytical description of arbitrary strongly aberrated axially symmetric focusing is developed. This is done by matching the solution of geometrical optics with a wave pattern which is universal for the underlying ray structure. The corresponding canonical integral is the Bessoid integral, which is a three dimensional generalization of the Pearcey integral that approximates the field near an arbitrary two-dimensional cusp. We first develop the description for scalar fields and then generalize it to the vector case. As a practical example the formalism is applied to the focusing of light by transparent dielectric spheres with a few wavelengths in diameter. The results demonstrate good agreement with the Mie theory down to Mie parameters of about 30. Compact analytical expressions are derived for the intensity on the axis and the position of the diffraction focus both for the general case and for the focusing by microspheres. The high intensity region is narrower than for an ideal lens of the same aperture at the expense of longitudinal localization and has a polarization dependent fine structure, which can be explained quantitatively. The results are relevant for aerosol and colloid science where natural light focusing occurs and can be used in laser micro- and nano-processing of materials. △ Less

Submitted 2 June, 2006; originally announced June 2006.

Comments: 17 pages, 11 figures

Journal ref: Phys. Rev. B 73, 235401 (2006)

Differential complexes and numerical stability

Authors: Douglas N. Arnold

Abstract: Differential complexes such as the de Rham complex have recently come to play an important role in the design and analysis of numerical methods for partial differential equations. The design of stable discretizations of systems of partial differential equations often hinges on capturing subtle aspects of the structure of the system in the discretization. In many cases the differential geometric… ▽ More Differential complexes such as the de Rham complex have recently come to play an important role in the design and analysis of numerical methods for partial differential equations. The design of stable discretizations of systems of partial differential equations often hinges on capturing subtle aspects of the structure of the system in the discretization. In many cases the differential geometric structure captured by a differential complex has proven to be a key element, and a discrete differential complex which is appropriately related to the original complex is essential. This new geometric viewpoint has provided a unifying understanding of a variety of innovative numerical methods developed over recent decades and pointed the way to stable discretizations of problems for which none were previously known, and it appears likely to play an important role in attacking some currently intractable problems in numerical PDE. △ Less

Submitted 30 November, 2002; originally announced December 2002.

Report number: ICM-2002 MSC Class: 65N12

Journal ref: Proceedings of the ICM, Beijing 2002, vol. 1, 137--157

New first-order formulation for the Einstein equations

Authors: A. M. Alekseenko, D. N. Arnold

Abstract: We derive a new first-order formulation for Einstein's equations which involves fewer unknowns than other first-order formulations that have been proposed. The new formulation is based on the 3+1 decomposition with arbitrary lapse and shift. In the reduction to first order form only 8 particular combinations of the 18 first derivatives of the spatial metric are introduced. In the case of lineari… ▽ More We derive a new first-order formulation for Einstein's equations which involves fewer unknowns than other first-order formulations that have been proposed. The new formulation is based on the 3+1 decomposition with arbitrary lapse and shift. In the reduction to first order form only 8 particular combinations of the 18 first derivatives of the spatial metric are introduced. In the case of linearization about Minkowski space, the new formulation consists of symmetric hyperbolic system in 14 unknowns, namely the components of the extrinsic curvature perturbation and the 8 new variables, from whose solution the metric perturbation can be computed by integration. △ Less

Submitted 16 May, 2004; v1 submitted 21 October, 2002; originally announced October 2002.

Comments: 6 pages

Journal ref: Phys.Rev. D68 (2003) 064013

The Advanced Photon Source Injector Test Stand Control System

Authors: J. F. Maclean, N. D. Arnold

Abstract: The Advanced Photon Source (APS) primary and backup injectors consist of two thermionic-cathode rf guns. These guns are being upgraded to provide improved performance, to improve ease of maintenance, and to reduce downtime required for repair or replacement of a failed injector. As part of the process, an injector test stand is being prepared. This stand is effectively independent of the APS lin… ▽ More The Advanced Photon Source (APS) primary and backup injectors consist of two thermionic-cathode rf guns. These guns are being upgraded to provide improved performance, to improve ease of maintenance, and to reduce downtime required for repair or replacement of a failed injector. As part of the process, an injector test stand is being prepared. This stand is effectively independent of the APS linac and will allow for complete characterization and validation of an injector prior to its installation into the APS linac. A modular control system for the test stand has been developed using standard APS control solutions with EPICS to deliver a flexible and comprehensive control system. The modularity of the system will allow both future expansion of test stand functionality and evaluation of new control techniques and solutions. △ Less

Submitted 1 February, 2002; v1 submitted 8 November, 2001; originally announced November 2001.

Comments: Poster paper (TUAP015) at ICALEPCS 2001, 3 pages, 2 figures, pdf

Journal ref: eConf C011127 (2001) TUAP015

Testing and Implementation Progress on the Advanced Photon Source (APS) Linear Accelerator (Linac) High-Power S-band Switching System

Authors: A. E. Grelick, N. Arnold, S. Berg, D. Dohan, G. Goeppner, Y. W. Kang, A. Nassiri, S. Pasky, G. Pile, T. Smith, S. J. Stein

Abstract: An S-band linear accelerator is the source of particles and the front end of the Advanced Photon Source injector. In addition, it supports a low-energy undulator test line (LEUTL) and drives a free-electron laser (FEL). A waveguide-switching and distribution system is now under construction. The system configuration was revised to be consistent with the recent change to electron-only operation.… ▽ More An S-band linear accelerator is the source of particles and the front end of the Advanced Photon Source injector. In addition, it supports a low-energy undulator test line (LEUTL) and drives a free-electron laser (FEL). A waveguide-switching and distribution system is now under construction. The system configuration was revised to be consistent with the recent change to electron-only operation. There are now six modulator-klystron subsystems, two of which are being configured to act as hot spares for two S-band transmitters each, so that no single failure will prevent injector operation. The two subsystems are also used to support additional LEUTL capabilities and off-line testing. Design considerations for the waveguide-switching subsystem, topology selection, control and protection provisions, high-power test results, and current status are described △ Less

Submitted 18 August, 2000; originally announced August 2000.

Comments: Linac 2000 paper No. THE07 3 pages with 3 figures

Journal ref: eConf C000821 (2000) THE07

RF System Upgrades to the Advanced Photon Source Linear Accelerator in Support of the Fel Operation

Authors: T. L. Smith, M. H. Cho, A. E. Grelick, G. Pile, A. Nassiri, N. Arnold

Abstract: The S-band linear accelerator, which was built to be the source of particles and the front end of the Advanced Photon Source injector, is now also being used to support a low-energy undulator test line (LEUTL) and to drive a free-electron laser (FEL). The more severe rf stability requirements of the FEL have resulted in an effort to identify sources of phase and amplitude instability and impleme… ▽ More The S-band linear accelerator, which was built to be the source of particles and the front end of the Advanced Photon Source injector, is now also being used to support a low-energy undulator test line (LEUTL) and to drive a free-electron laser (FEL). The more severe rf stability requirements of the FEL have resulted in an effort to identify sources of phase and amplitude instability and implement corresponding upgrades to the rf generation chain and the measurement system. Test data and improvements implemented and planned are described △ Less

Submitted 17 August, 2000; originally announced August 2000.

Comments: LC 2000 (3 pages, 6 figures)

Journal ref: eConf C000821 (2000) THc10

Approximation by quadrilateral finite elements

Authors: Douglas N. Arnold, Daniele Boffi, Richard S. Falk

Abstract: We consider the approximation properties of finite element spaces on quadrilateral meshes. The finite element spaces are constructed starting with a given finite dimensional space of functions on a square reference element, which is then transformed to a space of functions on each convex quadrilateral element via a bilinear isomorphism of the square onto the element. It is known that for affine… ▽ More We consider the approximation properties of finite element spaces on quadrilateral meshes. The finite element spaces are constructed starting with a given finite dimensional space of functions on a square reference element, which is then transformed to a space of functions on each convex quadrilateral element via a bilinear isomorphism of the square onto the element. It is known that for affine isomorphisms, a necessary and sufficient condition for approximation of order r+1 in L2 and order r in H1 is that the given space of functions on the reference element contain all polynomial functions of total degree at most r. In the case of bilinear isomorphisms, it is known that the same estimates hold if the function space contains all polynomial functions of separate degree r. We show, by means of a counterexample, that this latter condition is also necessary. As applications we demonstrate degradation of the convergence order on quadrilateral meshes as compared to rectangular meshes for serendipity finite elements and for various mixed and nonconforming finite elements. △ Less

Submitted 3 May, 2000; originally announced May 2000.

Comments: 13 pages, 2 figures, submitted to Math. Comp

MSC Class: 65N30 (Primary) 41A10; 41A25; 41A27; 41A63 (Secondary)

Journal ref: Math. Comp., 71(239):909-922, 2002

Adaptive Finite Elements and Colliding Black Holes

Authors: Douglas N. Arnold, Arup Mukherjee, Luc Pouly

Abstract: According to the theory of general relativity, the relative acceleration of masses generates gravitational radiation. Although gravitational radiation has not yet been detected, it is believed that extremely violent cosmic events, such as the collision of black holes, should generate gravity waves of sufficient amplitude to detect on earth. The massive Laser Interferometer Gravitational-wave Obs… ▽ More According to the theory of general relativity, the relative acceleration of masses generates gravitational radiation. Although gravitational radiation has not yet been detected, it is believed that extremely violent cosmic events, such as the collision of black holes, should generate gravity waves of sufficient amplitude to detect on earth. The massive Laser Interferometer Gravitational-wave Observatory, or LIGO, is now being constructed to detect gravity waves. Consequently there is great interest in the computer simulation of black hole collisions and similar events, based on the numerical solution of the Einstein field equations. In this note we introduce the scientific, mathematical, and computational problems and discuss the development of a computer code to solve the initial data problem for colliding black holes, a nonlinear elliptic boundary value problem posed in an unbounded three dimensional domain which is a key step in solving the full field equations. The code is based on finite elements, adaptive meshes, and a multigrid solution process. Here we will particularly emphasize the mathematical and algorithmic issues arising in the generation of adaptive tetrahedral meshes. △ Less

Submitted 15 September, 1997; originally announced September 1997.

Comments: 15 pages, 6 figures, to be published in "Numerical Analysis 1997: Proceedings of the 17th Biennial Conference on Numerical Analysis"