Small Matrix Research Articles

Time performance of Analog Pixel Test Structures with in-chip operational amplifier implemented in 65 nm CMOS imaging process

In the context of the CERN EP R&D on monolithic sensors and the ALICE ITS3 upgrade, the Tower Partners Semiconductor Co (TPSCo) 65nm process has been qualified for use in high energy physics, and adopted for the ALICE ITS3 upgrade. An Analog Pixel Test Structure (APTS) featuring fast per pixel operational-amplifier-based buffering for a small matrix of four by four pixels, with a sensor with a small collection electrode and a very non-uniform electric field, was designed to allow detailed characterization of the pixel performance in this technology. Several variants of this chip with different pixel designs have been characterized with a 120GeV/c positive hadron beam. This result indicates that the APTS-OA prototype variants with the best performance achieve a time resolution of 63ps with a detection efficiency exceeding 99% and a spatial resolution of 2μm, highlighting the potential of TPSCo 65 nm CMOS imaging technology for high-energy physics and other fields requiring precise time measurement, high detection efficiency, and excellent spatial resolution.

Mechanical and Geometric Characterization of a Novel 2-Ply Vacuum-Pressed Biological Scaffold Patch Design for Posterior Mitral Valve Reconstruction.

To assess the mechanical properties of small intestinal submucosal extracellular matrix (SIS-ECM) iterations and choose the optimal version for evaluating functional geometrics after posterior mitral valve reconstruction. Four SIS-ECM versions (2- and 4-ply vacuum-pressed and lyophilized) underwent uniaxial tensile testing. A posterior mitral valve reconstruction patch was developed based on MRI scans (n = 5). Posterior mitral valve reconstruction using 2-ply vacuum-pressed SIS-ECM was performed (n = 7), and geometrics were evaluated using a modified left heart simulator. The vacuum-pressed iterations displayed superior maximum stress values compared to lyophilized (2-ply: median [IQR], 15.8 [15.2-19.0] vs 7.9 [7.3-8.3] MPa, p < 0.001; 4-ply: median (IQR), 15.8 -[14.6-22.0] vs 7.9 [7.6-8.4] MPa). All reconstructed valves were competent with preserved total leaflet area, but individual leaflet segment areas were redistributed. Posterior mitral valve reconstruction with our 2-ply vacuum-pressed SIS-ECM patch design was feasible in vitro. Further in vivo evaluation is warranted.

COMPLEX WEIGHING MATRICES AND QUATERNARY CODES

Abstract Weighing matrices with entries in the complex cubic and sextic roots of unity are employed to construct Hermitian self-dual codes and Hermitian linear complementary dual codes over the finite field $\mathrm {GF}(4).$ The parameters of these codes are explored for small matrix orders and weights.

Maxwell boundary condition for discrete velocity methods: Macroscopic physical constraints and Lagrange multiplier-based implementation

In this paper, we propose an algorithm that imposes macroscopic physical constraints with Lagrange multiplier approach in implementing the Maxwell boundary condition within the framework of the discrete velocity method. For the simulation of rarefied gas flows in the presence of solid walls with complex geometry, the distribution function in the reflection region of the wall surface needs to be constructed in the discrete velocity space, to fulfill the specular reflection in the Maxwell boundary condition. The construction process should not consist of interpolation only, but include certain macroscopic physical constraints at the wall surface, so as to correctly account for gas-surface interaction on a macroscopic level. We demonstrate that for the specular reflection, keeping the symmetry of the first three moments of the distribution function between the incident and reflected region is sufficient for maintaining the conservation of mass, momentum, and energy at the wall surface. Furthermore, to strictly satisfy macroscopic physical constraints, a Lagrange multiplier method is introduced into the construction of the distribution function to correct the pure interpolation solution. In addition, the construction process requires the inversion of a large and sparse matrix (of dimension N × N, where N is the number of points in the velocity space). To improve the computational efficiency, the matrix inversion is converted into that of a much smaller matrix, i.e. (D + 2) × (D + 2) in the d-dimensional physical space. A series of numerical experiments are conducted to examine the performance of the proposed algorithm under different flow conditions. We demonstrate that the results obtained by the proposed algorithm are more accurate than the pure interpolation solution, comparing with the benchmark data. Moreover, after the validation of our results with previous studies, we find that the method significantly enhances the conservation of total mass and energy, especially for flows in an enclosed domain.

Two-Tier Modular Anharmonic Small Matrix Path Integral with Composite Spin-Boson Baths.

The anharmonic small matrix path integral (anh-SMatPI) algorithm is adapted to composite environments, where each anharmonic bath unit couples to the system and to its own harmonic bath. The propagation matrices are obtained from numerical evaluation of the composite influence functional, which is based on separate iterative path integral calculations for each bath module. The numerically exact nature of the algorithm is demonstrated on a two-level system (TLS) coupled to dissipative TLS baths of one, two, and 50 units.

Finite element analysis and experimental study on the cutting mechanism of SiCp/Al composite in laser ultrasonic elliptic vibration turning

SiCp/Al composites have excellent material properties and are increasingly being used in the aerospace, military, and electronic packaging industries. However, the inhomogeneity and non-conductivity of conventional turning (CT) results in poor material surface quality and high cutting forces, which seriously hinder the application of particlereinforced metal matrix composites. The thermal softening property of laser-assisted turning (LAT) and the intermittent cutting property of two-dimensional (2D) ultrasonic elliptical vibratory turning (UEVT) offer unique advantages in improving material surface quality. Therefore, a novel laser ultrasonic elliptical vibratory turning (LUEVT) machining technique is proposed to improve the machining of SiCp/Al composites with two different volume fractions. The effect of highfrequency intermittent machining on material processing was further investigated by adjusting the pulsed laser power. Finite element modelling was used to predict the machining state, surface roughness, and chip morphology between the workpiece and the PCD tool. The effects of varying the machining parameters during the experiment on the two composites were analysed. The results showed that the LUEVT of 25 % and 45 % SiCp/Al composites were reduced by 39.3 % and 30.7 % respectively compared to the CT cutting forces. The experiments verified the consistency of the surface roughness and chip morphology with the finite element simulation results. The stability of the cutting force during the cutting process is also improved, as the material removal process mainly takes the form of small particle fragmentation and matrix encapsulation.

Simulation of soil-structure interaction using inelasticity-separated finite element method

The soil–structure interaction (SSI) effect is nonnegligible during the nonlinear seismic response analysis of large-scale or underground structures. However, the computation of nonlinear response of a structure considering SSI effect usually is a time-consuming process because the numerical model should involve the structure itself, a wide-range soil domain and the soil–structure interface, especially for large-scale structure. This study aims to incorporate the inelasticity-separated finite element method (IS-FEM), which is an algorithm recently developed for efficient structural nonlinear analysis, to improve the efficiency of the soil-structure interaction system. To this end, an inelasticity-separated contact element model for modeling the nonlinear behavior of soil and structure interface is developed in this study. To derive the inelasticity-separated governing equation of this model, a reference elastic stiffness is defined so that the normal and shear deformation of any point pair in the presented contact element is decomposed into reference linear and reference nonlinear parts according to this reference elastic stiffness. As a result, the problem of inconsistency between the requirement of IS-FEM for constant initial linear elastic stiffness and the existence of inflection point for the stress versus relative displacement relationship of contact behavior can be solved. Then, the Goodman element formulation for depicting the nonlinear contact behavior is introduced and an interpolation scheme to establish the reference nonlinear deformation field in an element is incorporated. By combining the proposed model with existing inelasticity-separated elements, a global analytical model of the soil–structure interaction system can be established within the framework of IS-FEM. Because the nonlinearity usually occur in some local regions of the global model and the use of the IS-FEM only require performing these operations for a small scale matrix representing local nonlinearity, the computational efficiency of nonlinear structural response analysis considering SSI effect can be improved significantly. The applicability and efficiency of the proposed method is finally verified by two numerical examples.

Cranial anatomy and phylogenetic affinities of Bolosaurus major, with new information on the unique bolosaurid feeding apparatus and evolution of the impedance-matching ear.

Resolving the phylogenetic relationships of early amniotes, in particular stem reptiles, remains a difficult problem. Three-dimensional morphological analysis of well-preserved stem-reptile specimens can reveal important anatomical data and clarify regions of phylogeny. Here, we present the first thorough description of the unusual early Permian stem reptile Bolosaurus major, including the first comprehensive description of a bolosaurid braincase. We describe previously obscured details of the palate, allowing for insight into bolosaurid feeding mechanics. Aspects of the rostrum, palate, mandible, and neurocranium suggest that B. major had a particularly strong bite. We additionally found B. major has a surprisingly slender stapes, similar to that of the middle Permian stem reptile Macroleter poezicus, which may suggest enhanced hearing abilities compared to other Paleozoic amniotes (e.g., captorhinids). We incorporated our new anatomical information into a large phylogenetic matrix (150 OTUs, 590 characters) to explore the relationship of Bolosauridae among stem reptiles. Our analyses generally recovered a paraphyletic "Parareptilia," and found Bolosauridae to diverge after Captorhinidae + Araeoscelidia. We also included B. major within a smaller matrix (10 OTUs, 27 characters) designed to explore the interrelationships of Bolosauridae and found all species of Bolosaurus to be monophyletic. While reptile relationships still require further investigation, our phylogeny suggests repeated evolution of impedance-matching ears in Paleozoic stem reptiles.

Study on dynamic mechanism of granular flow erosion and entrainment based on DEM theory

BackgroundGranular flows are common on the Qinghai–Tibet Plateau and the Hengduan Mountains in China, and their dynamic process processes have obvious erosional and entrainment effects. On the one hand, the volume of the granular flow increases by a factor of several or ten, which significantly increases its ability to cause a catastrophe; on the other hand, the eroded loose material affects the granular flow dynamics process and changes its state of motion.MethodsIn this paper, the dynamic mechanism of granular flow erosion and entrainment is investigated by DEM simulation.Purpose The effects of different substrate materials and substrate boundary conditions on granular flow erosion and entrainment are analyzed, and the effects of material mixing caused by erosion and entrainment on the state of motion of granular flow are discussed. It was verified that the kinetic mechanisms of granular flow erosion and entrainment includes impact erosion, ploughing, and shear abrasion.ResultsAnd discovered that small matrix particle size, small matrix boundary friction, and small matrix thickness lead to stronger ploughing and shear abrasion; Large matrix fractal dimensions result in stronger ploughing and weaker shear abrasion, and the granular flow does not entrain large amounts of material to the accumulation zone. Meanwhile, the dynamics of erosion and entrainment of granular flow were investigated, and the results showed that: 1. The greater the erosion rate, the greater the velocity and kinetic energy of the granular flow, the greater the distance traveled, and the smaller the apparent friction angle (i.e., the greater the mobility). 2. The amount of small granules in a granular flow changes its fluidity, the more small granules there are, the more fluid it is. 3. The fit reveals that the substrate fractal dimension has the strongest effect on the velocity and kinetic energy of granular flow, followed by substrate thickness and substrate boundary friction.

A reduced order variational spectral method for efficient construction of eigenstrain‐based reduced order homogenization models

AbstractReduced order models (ROMs) are often coupled with concurrent multiscale simulations to mitigate the computational cost of nonlinear computational homogenization methods. Construction (or training) of ROMs typically requires evaluation of a series of linear or nonlinear equilibrium problems, which itself could be a computationally very expensive process. In the eigenstrain‐based reduced order homogenization method (EHM), a series of linear elastic microscale equilibrium problems are solved to compute the localization and interaction tensors that are in turn used in the evaluation of the reduced order multiscale system. These microscale equilibrium problems are typically solved using either the finite element method or semi‐analytical methods. In the present study, a reduced order variational spectral method is developed for efficient computation of the localization and interaction tensors. The proposed method leads to a small stiffness matrix that scales with the order of the reduced basis rather than the number of degrees of freedom in the finite element mesh. The reduced order variational spectral method maintains high accuracy in the computed response fields. A speedup higher than an order of magnitude can be achieved compared to the finite element method in polycrystalline microstructures. The accuracy and scalability of the method for large polycrystals and increasing phase property contrast are investigated.

One-Year Outcome of an Ongoing Pre-Clinical Growing Animal Model for a Tissue-Engineered Valved Pulmonary Conduit.

Objectives: A self-constructed valved pulmonary conduit made out of a de-cellularized porcine small intestinal submucosal extracellular matrix biological scaffold was tested in a chronic growing lamb model. Methods: The conduit was implanted in pulmonary valve position in 19 lambs. We monitored clinical, laboratory, and echocardiographic findings until 12 months after surgery. In two animals, euthanasia was planned at nine and twelve months. Pre-mortem chest computed tomography and post-mortem pathologic work up were performed. Data are presented as frequency and percentage, median and range, or mean and standard deviation. Results: Twelve (63.2%) animals survived the perioperative period. Three unexpected deaths occurred during the follow-up period: one due to aspiration pneumonia at 23 days after surgery, and two due to early and late infective endocarditis of the conduit at 18 and 256 days. In the two animals with planned scarification, the pre-mortem CT scan revealed mild or no calcification within the conduit or valve leaflets. In the echocardiographic examination at 12 months, peak and mean systolic pressure gradients across the conduit valve were 6.5 (3-21) mmHg and 3 (2-12) mmHg, while valve regurgitation was none (n = 2), trivial (n = 5), moderate (n = 1), or severe (n = 1). No clinical or laboratory signs of hemolysis were seen. After 12 months of follow-up, the animals' body weights had increased from 33 (27-38) kg to 53 (38-66) kg (p = 0.010). Conclusions: Implantation of a valved pulmonary conduit in our growing lamb model was feasible. Infective endocarditis of the implanted valved conduit remained a significant complication.

Parsing the Influence Functional: Harmonic Bath Mapping and Anharmonic Small Matrix Path Integral.

The influence functional (IF) encodes all of the information required for calculating the dynamical properties of a system in contact with its environment. A direct and simple procedure is introduced for extracting from a few numerical evaluations of the IF, without computing time correlation functions or evaluating integrals, the parameters required for path integral calculations, either within or beyond the harmonic mapping, and for assessing the accuracy of the harmonic bath approximation. In addition, the small matrix decomposition of the path integral (SMatPI) is extended to anharmonic environments and the required matrices are constructed directly from the IF.

Flexible waveband routing node architecture for spatial-division and multi-band optical networks [Invited

We propose a network architecture named flexible waveband routing (FWR) that adopts flexible path bundling and bundled path routing. This two-stage routing scheme can be implemented by a combination of small degree wavelength selective switches (WSSs) and optical matrix switches. As the optical-cross-connect (OXC) scale depends on just that of the matrix switches adopted, large-scale OXC nodes are realized by introducing recently developed high-port-count switches. FWR supports spatial division multiplexing efficiently through its introduction of spatially joint switching at WSSs. Multi-band transmission can also be supported by the use of WSSs handling different frequency bands. Numerical simulations and transmission experiments confirm its validity.

Entire Mitral Valve Reconstruction Using Porcine Extracellular Matrix: Adding a Ring Annuloplasty

PurposeThis study investigated the implications of inserting a flexible annuloplasty ring after reconstructing the entire mitral valve in a porcine model using a previously investigated tube graft design made of 2-ply small intestinal submucosa extracellular matrix (CorMatrix®).MethodsAn acute model with eight 80-kg pigs, each acting as its own control, was used. The entire mitral valve was reconstructed with a 2-ply small intestinal submucosa extracellular matrix tube graft (CorMatrix®). Subsequently, a Simulus® flexible ring was inserted. The characterization was based on mitral annular geometry and valvular dynamics with sonomicrometry and echocardiography.ResultsAfter adding the ring annuloplasty, the in-plane annular dynamics were more constant throughout the cardiac cycle compared to the reconstruction alone. However, the commissure–commissure distance was statistically significantly decreased [35.0 ± 3.4 mm vs. 27.4 ± 1.9 mm, P < 0.001, diff = − 7.6 mm, 95% CI, − 9.8 to (−5.4) mm] after ring insertion, changing the physiological annular D-shape into a circular shape which created folds at the coaptation zone resulting in a central regurgitant jet on color Doppler.ConclusionWe successfully reconstructed the entire mitral valve using 2-ply small intestinal submucosal extracellular matrix (CorMatrix®) combined with a flexible annuloplasty. The annuloplasty reduced the unphysiological systolic widening previously found with this reconstructive technique. However, the Simulus flex ring changed the physiological annular D-shape into a circular shape and hindered a correct unfolding of the leaflets. Thus, we do not recommend a flexible ring in conjunction with this reconstructive technique; further investigations are needed to discover a more suitable remodelling annuloplasty.

Band structure calculations of three-dimensional solid-fluid coupling phononic crystals using dual reciprocity boundary element method and wavelet compression method

The algorithm for calculating the band structures of two-dimensional phononic crystals (PCs) using the boundary element method (BEM) has been proposed for many years. However, it has not yet been extended to three-dimensional (3D) PCs because the fundamental solutions of 3D dynamics are complex and are related to angular frequency. In this study, the BEM is applied to calculate the band structures of 3D solid-fluid coupling PCs by combining the dual reciprocity method. The use of the dual reciprocity method can help avoid nonlinear eigenvalue problems. To address the limitation of fully populated matrices in BEM, the wavelet compression method based on B-spline wavelet on the interval is adopted. Some small matrix entries, which are generated by the vanishing moment and local support characteristics, are set to zero using the provided truncation technique. This process results in the production of sparse matrices. The constructed generalized linear eigenvalue equations are modified to tackle the ill-conditioned matrix issues arising from the distinct fundamental solutions of solids and fluids. The results show that this method is superior to the finite element method in terms of calculation efficiency. Compared to conventional BEM, the current wavelet BEM can not only generate sparse matrices but also reduce the integration calculation time when handling large-scale problems.

Kink Sum for Long-Memory Small Matrix Path Integral Dynamics.

The small matrix decomposition of the real-time path integral (SMatPI) allows for numerically exact and efficient propagation of the reduced density matrix (RDM) for system-bath Hamiltonians. Its high efficiency lies in the small size of the SMatPI matrices employed in the iterative algorithm, whose size is equal to that of the full RDM. By avoiding the storage and multiplication of large tensors, the SMatPI algorithm is applicable in multistate systems under a variety of conditions. The main computational effort is the evaluation of path sums within the entangled memory length to construct the SMatPI matrices. A number of methods are available for this task, each with its own favorable parameter regime, but calculations with strong system-bath coupling and long memory at low temperatures remain out of reach. The present paper evaluates the path sums by binning the paths (in forward time only) based on their amplitudes, which depend on the number and type of kinks they contain. The algorithm is very efficient, leading to a dramatic acceleration of path sums and significantly extending the accessible memory length in the most challenging regimes.

Comparative evaluation of two different xenogenic acellular matrices on full-thickness skin wound healing.

The purpose of the study was to compare the healing potential of bubaline small intestinal matrix (bSIM) and fish swim bladder matrix (FSBM) on full-thickness skin wounds in rabbits. Four full-thickness skin wounds (each 20×20mm) were created on the dorsum of 18 rabbits that were divided into three groups based on treatment: untreated sham control (I), implanted with double layers of bSIM (II) and implanted with double layers of FSBM (III). Macroscopic, immunologic and histologic observations were made to evaluate wound healing. Gross healing progression in the bSIM and FSBM groups showed significantly (p<0.05) less wound contraction compared with the sham group. The IgG concentration in rabbit sera was significantly (p<0.05) lower in the FSBM group compared with the bSIM group by enzyme-linked immunosorbent assay. The stimulation index of peripheral blood lymphocytes was significantly (p<0.05) lower in the FSBM group compared with the bSIM group by the MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) assay. Implantation of FSBM resulted in improved re-epithelialisation, neovascularisation and fibroplasia. The FSBM is a more effective dermal substitute when compared with the bSIM for full-thickness skin wound repair in rabbit.

Differential expression of mast cells in the small airways and alveolar septa of current smokers and patients with small airway disease and COPD.

COPD patients suffer from dysregulated and suppressed immune functionality, determined by their loss of degranulating capacity. Here we provide crucial information on the presence of degranulated mast cells (MCs) in COPD airways and demonstrate their relationship to lung physiology and airway remodelling. Small airway lung resections from non-smoking controls (NC), normal lung function smokers (NLFS), small airway disease (SAD), and mild-to-moderate COPD current smokers (COPD-CS) and ex-smokers (COPD-ES) were dual immuno-stained with MC tryptase and degranulation marker lysosome-associated membrane protein (LAMP)-1. Total MCs, degranulating MCs and non-MCs were enumerated in small airway epithelium and subepithelium, and in alveolar septa. In the small airway wall subepithelial areas, COPD-CS and COPD-ES patients had significantly lower MCs than the NC group (p<0.05), although the numbers were considerably higher in the small airway epithelium (p<0.01). Degranulating non-MCs were higher in SAD (p<0.05) than in COPD in the small airway subepithelium. In contrast, there were significant increases in total MCs (degranulated and non-degranulated) and degranulated non-MCs in the alveolar septum of COPD patients compared with the NC group (p<001). The lower numbers of MCs in the subepithelium correlated with lower forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) and forced expiratory flow at 25-75% of FVC (FEF25-75%), higher smoking rates in COPD patients, and increased small airway wall thickness and extracellular matrix. The increase in MCs in the alveolar septum negatively correlated with FEF25-75%. This study is the first to assess the differential pattern of MC, degranulating MC and non-MC populations in the small airways and alveoli of COPD patients. The spatial positioning of the MCs within the airways showed variable correlations with lung function.

GRK2-YAP signaling is implicated in pulmonary arterial hypertension development.

Pulmonary arterial hypertension (PAH) is characterized by excessive proliferation of small pulmonary arterial vascular smooth muscle cells (PASMCs), endothelial dysfunction, and extracellular matrix remodeling. G protein-coupled receptor kinase 2 (GRK2) plays an important role in the maintenance of vascular tone and blood flow. However, the role of GRK2 in the pathogenesis of PAH is unknown. GRK2 levels were detected in lung tissues from healthy people and PAH patients. C57BL/6 mice, vascular smooth muscle cell-specific Grk2 -knockout mice ( Grk2ΔSM22 ), and littermate controls ( Grk2flox/flox ) were grouped into control and hypoxia mice ( n =8). Pulmonary hypertension (PH) was induced by exposure to chronic hypoxia (10%) combined with injection of the SU5416 (cHx/SU). The expression levels of GRK2 and Yes-associated protein (YAP) in pulmonary arteries and PASMCs were detected by Western blotting and immunofluorescence staining. The mRNA expression levels of Grk2 and Yes-associated protein ( YAP ) in PASMCs were quantified with real-time polymerase chain reaction (RT-PCR). Wound-healing assay, 3-(4,5)-dimethylthiahiazo (-z-y1)-3,5-di-phenytetrazoliumromide (MTT) assay, and 5-Ethynyl-2'-deoxyuridine (EdU) staining were performed to evaluate the proliferation and migration of PASMCs. Meanwhile, the interaction among proteins was detected by immunoprecipitation assays. The expression levels of GRK2 were upregulated in the pulmonary arteries of patients with PAH and the lungs of PH mice. Moreover, cHx/SU-induced PH was attenuated in Grk2ΔSM22 mice compared with littermate controls. The amelioration of PH in Grk2ΔSM22 mice was accompanied by reduced pulmonary vascular remodeling. In vitro study further confirmed that GRK2 knock-down significantly altered hypoxia-induced PASMCs proliferation and migration, whereas this effect was severely intensified by overexpression of GRK2 . We also identified that GRK2 promoted YAP expression and nuclear translocation in PASMCs, resulting in excessive PASMCs proliferation and migration. Furthermore, GRK2 is stabilized by inhibiting phosphorylating GRK2 on Tyr86 and subsequently activating ubiquitylation under hypoxic conditions. Our findings suggest that GRK2 plays a critical role in the pathogenesis of PAH, via regulating YAP expression and nuclear translocation. Therefore, GRK2 serves as a novel therapeutic target for PAH treatment.

A high-throughput small volume matrix based calibration using isotope dilution liquid chromatography tandem mass spectrometry analysis for 42 per and polyfluoroalkyl substances in groundwater

A novel method for the determination of per- and polyfluoroalkyl substances (PFAS) in groundwater is presented using a subsample, matrix-matched calibrators, 96-well plate solid phase extraction (SPE), and ultra-performance liquid chromatography tandem mass spectrometry (UPLC-MS/MS). Accuracy, precision, measurement of uncertainty (MOU), method detection limit (MDL), method quantitation limit (MQL), analytical measurement range, interferences/ion suppression, and analyte stability were determined as part of the in-house method validation. The method quantitates 42 PFAS compounds from nine different compound classes. Accuracy for the reference material (RM) and matrix spike (MS) ranged from 52.3 to 117.8 %, and precision for the MS and matrix spike duplicate (MSD) had a coefficient of variation (CV) from 2.0 % to 23.3 %. MDLs spanned from 0.07 to 1.97 ng L−1, with MQLs ranging from 0.20 to 5.90 ng L−1. Suppression studies determined that iron and manganese have effects on analytes that do not have paired isotopically labeled standards. The results from the in-house validation indicated that this Michigan Department of Health and Human Services laboratory developed test meets the necessary accuracy, precision, MDL, MQL and reporting limits requirement established by the laboratory's quality system essentials (QSEs) and select criteria from the Department of Defense (DoD) Quality Systems Manual for Environmental Laboratories and American Industrial Hygiene Association Laboratory Accreditation Program, LLC (AIHA LAP, LLC) accrediting International Standard Organization (ISO/IEC 17025:2017) check list.