Energy-norm-based and goal-oriented automatic hp adaptivity for electromagnetics: Application to waveguide Discontinuities

The finite-element method (FEM) enables the use of adapted meshes. The simultaneous combination of h (local variations in element size) and p (local variations in the polynomial order of approximation) refinements, i.e., hp-adaptivity, is the most powerful and flexible type of adaptivity. In this paper, two versions of a fully automatic hp-adaptive FEM for electromagnetics are presented. The first version is based on minimizing the energy-norm of the error. The second, namely the goal oriented strategy, is based on minimizing the error of a given (user-prescribed) quantity of interest. The adaptive strategy delivers exponential convergence rates for the error, even in the presence of singularities. The hp adaptivity is presented in the context of 2-D analysis of H -plane rectangular waveguide discontinuities. Stabilized variational formulations and H(curl) FEM discretizations in terms of quadrilaterals of variable order of approximation supporting anisotropy and hanging nodes are used. Comparison of energy-norm and goal-oriented hp-adaptive strategies in the context of waveguiding problems is provided. Specifically, the scattering parameters of the discontinuity are used as goal

3D hp-Adaptive Finite Element Simulations of Bend, Step, and Magic-T Electromagnetic Waveguide Structures

Metallic rectangular waveguides are often the preferred choice on telecommunication systems and medical equipment working on the upper microwave and millimeter wave frequency bands due to the simplicity of its geometry, low losses, and the capacity to handle high powers. Waveguide translational symmetry is interrupted by the unavoidable presence of bends, transitions, and junctions, among others. This paper employs a 3D hp self-adaptive grid-refinement finite element strategy for the solution of these relevant electromagnetic waveguide problems. These structures often incorporate dielectrics, metallic screws, round corners, and so on, which may facilitate its construction or improve its design, but significantly difficult its modeling when employing semi-analytical techniques. The hp-adaptive finite element method enables accurate modeling of these structures even in the presence of complex materials and geometries. Numerical results demonstrate the suitability of the hp-adaptive method for modeling these waveguide structures, delivering errors below 0.5% with a limited number of unknowns. Solutions of waveguide problems obtained with the self-adaptive hp-FEM are of comparable accuracy to those obtained with semi-analytical techniques such as the Mode Matching method, for problems where the latest methods can be applied. At the same time, the hp-adaptive FEM enables accurate modeling of more complex waveguide structures.TEC2010-18175/TCM MTM2010-1651

Fast parallel IGA-ADS solver for time-dependent Maxwell's equations

We propose a simulator for time-dependent Maxwell's equations with linear computational cost. We employ B-spline basis functions as considered in the isogeometric analysis (IGA). We focus on non-stationary Maxwell's equations defined on a regular patch of elements. We employ the idea of alternating-directions splitting (ADS) and employ a second-order accurate time-integration scheme for the time-dependent Maxwell's equations in a weak form. After discretization, the resulting stiffness matrix exhibits a Kronecker product structure. Thus, it enables linear computational cost LU factorization. Additionally, we derive a formulation for absorbing boundary conditions (ABCs) suitable for direction splitting. We perform numerical simulations of the scattering problem (traveling pulse wave) to verify the ABC. We simulate the radiation of electromagnetic (EM) waves from the dipole antenna. We verify the order of the time integration scheme using a manufactured solution problem. We then simulate magnetotelluric measurements. Our simulator is implemented in a shared memory parallel machine, with the GALOIS library supporting the parallelization. We illustrate the parallel efficiency with strong and weak scalability tests corresponding to non-stationary Maxwell simulations.EXPERTIA (KK-2021/00048) SIGZE (KK-2021/00095) PDC2021-121093-I0

Diagnostic ability of multifocal electroretinogram in early multiple sclerosis using a new signal analysis method

Purpose To determine if a novel analysis method will increase the diagnostic value of the multifocal electroretinogram (mfERG) in diagnosing early-stage multiple sclerosis (MS). Methods We studied the mfERG signals of OD (Oculus Dexter) eyes of fifteen patients diagnosed with early-stage MS (in all cases < 12 months) and without a history of optic neuritis (ON) (F: M = 11:4), and those of six controls (F:M = 3:3). We obtained values of amplitude and latency of N1 and P1 waves, and a method to assess normalized root-mean-square error (FNRMSE) between model signals and mfERG recordings was used. Responses of each eye were analysed at a global level, and by rings, quadrants and hemispheres. AUC (area under the ROC curve) is used as discriminant factor. Results The standard method of analysis obtains further discrimination between controls and MS in ring R3 (AUC = 0.82), analysing N1 waves amplitudes. In all of the retina analysis regions, FNRMSE value shows a greater discriminating power than the standard method. The highest AUC value (AUC = 0.91) was in the superior temporal quadrant. Conclusion By analysing mfERG recordings and contrasting them with those of healthy controls it is possible to detect early-stage MS in patients without a previous history of ON

Measurement of the polarisation of W bosons produced with large transverse momentum in pp collisions at sqrt(s) = 7 TeV with the ATLAS experiment

This paper describes an analysis of the angular distribution of W->enu and W->munu decays, using data from pp collisions at sqrt(s) = 7 TeV recorded with the ATLAS detector at the LHC in 2010, corresponding to an integrated luminosity of about 35 pb^-1. Using the decay lepton transverse momentum and the missing transverse energy, the W decay angular distribution projected onto the transverse plane is obtained and analysed in terms of helicity fractions f0, fL and fR over two ranges of W transverse momentum (ptw): 35 < ptw < 50 GeV and ptw > 50 GeV. Good agreement is found with theoretical predictions. For ptw > 50 GeV, the values of f0 and fL-fR, averaged over charge and lepton flavour, are measured to be : f0 = 0.127 +/- 0.030 +/- 0.108 and fL-fR = 0.252 +/- 0.017 +/- 0.030, where the first uncertainties are statistical, and the second include all systematic effects.Comment: 19 pages plus author list (34 pages total), 9 figures, 11 tables, revised author list, matches European Journal of Physics C versio

Observation of a new chi_b state in radiative transitions to Upsilon(1S) and Upsilon(2S) at ATLAS

The chi_b(nP) quarkonium states are produced in proton-proton collisions at the Large Hadron Collider (LHC) at sqrt(s) = 7 TeV and recorded by the ATLAS detector. Using a data sample corresponding to an integrated luminosity of 4.4 fb^-1, these states are reconstructed through their radiative decays to Upsilon(1S,2S) with Upsilon->mu+mu-. In addition to the mass peaks corresponding to the decay modes chi_b(1P,2P)->Upsilon(1S)gamma, a new structure centered at a mass of 10.530+/-0.005 (stat.)+/-0.009 (syst.) GeV is also observed, in both the Upsilon(1S)gamma and Upsilon(2S)gamma decay modes. This is interpreted as the chi_b(3P) system.Comment: 5 pages plus author list (18 pages total), 2 figures, 1 table, corrected author list, matches final version in Physical Review Letter

Search for displaced vertices arising from decays of new heavy particles in 7 TeV pp collisions at ATLAS

We present the results of a search for new, heavy particles that decay at a significant distance from their production point into a final state containing charged hadrons in association with a high-momentum muon. The search is conducted in a pp-collision data sample with a center-of-mass energy of 7 TeV and an integrated luminosity of 33 pb^-1 collected in 2010 by the ATLAS detector operating at the Large Hadron Collider. Production of such particles is expected in various scenarios of physics beyond the standard model. We observe no signal and place limits on the production cross-section of supersymmetric particles in an R-parity-violating scenario as a function of the neutralino lifetime. Limits are presented for different squark and neutralino masses, enabling extension of the limits to a variety of other models.Comment: 8 pages plus author list (20 pages total), 8 figures, 1 table, final version to appear in Physics Letters

Measurement of the inclusive isolated prompt photon cross-section in pp collisions at sqrt(s)= 7 TeV using 35 pb-1 of ATLAS data

A measurement of the differential cross-section for the inclusive production of isolated prompt photons in pp collisions at a center-of-mass energy sqrt(s) = 7 TeV is presented. The measurement covers the pseudorapidity ranges |eta|<1.37 and 1.52<=|eta|<2.37 in the transverse energy range 45<=E_T<400GeV. The results are based on an integrated luminosity of 35 pb-1, collected with the ATLAS detector at the LHC. The yields of the signal photons are measured using a data-driven technique, based on the observed distribution of the hadronic energy in a narrow cone around the photon candidate and the photon selection criteria. The results are compared with next-to-leading order perturbative QCD calculations and found to be in good agreement over four orders of magnitude in cross-section.Comment: 7 pages plus author list (18 pages total), 2 figures, 4 tables, final version published in Physics Letters

Measurement of D*+/- meson production in jets from pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

This paper reports a measurement of D*+/- meson production in jets from proton-proton collisions at a center-of-mass energy of sqrt(s) = 7 TeV at the CERN Large Hadron Collider. The measurement is based on a data sample recorded with the ATLAS detector with an integrated luminosity of 0.30 pb^-1 for jets with transverse momentum between 25 and 70 GeV in the pseudorapidity range |eta| < 2.5. D*+/- mesons found in jets are fully reconstructed in the decay chain: D*+ -> D0pi+, D0 -> K-pi+, and its charge conjugate. The production rate is found to be N(D*+/-)/N(jet) = 0.025 +/- 0.001(stat.) +/- 0.004(syst.) for D*+/- mesons that carry a fraction z of the jet momentum in the range 0.3 < z < 1. Monte Carlo predictions fail to describe the data at small values of z, and this is most marked at low jet transverse momentum.Comment: 10 pages plus author list (22 pages total), 5 figures, 1 table, matches published version in Physical Review

Search for scalar top quark pair production in natural gauge mediated supersymmetry models with the ATLAS detector in pp collisions at sqrt(s) = 7 TeV

The results of a search for pair production of the lighter scalar partners of top quarks in 2.05 fb-1 of pp collisions at sqrt(s) =7 TeV using the ATLAS experiment at the LHC are reported. Scalar top quarks are searched for in events with two same flavour opposite-sign leptons (electrons or muons) with invariant mass consistent with the Z boson mass, large missing transverse momentum and jets in the final state. At least one of the jets is identified as originating from a b-quark. No excess over Standard Model expectations is found. The results are interpreted in the framework of R-parity conserving, gauge mediated Supersymmetry breaking `natural' scenarios, where the neutralino is the next-to-lightest supersymmetric particle. Scalar top quark masses up to 310 GeV are excluded for the lightest neutralino mass between 115 GeV and 230 GeV at 95% confidence level, reaching an exclusion of the scalar top quark mass of 330 GeV for the lightest neutralino mass of 190 GeV. Scalar top quark masses below 240 GeV are excluded for all values of the lightest neutralino mass above the Z boson mass.Comment: 7 pages plus author list (20 pages total), 4 figures, 1 table, matches published PLB versio