Variations of the FEAST Eigenvalue Algorithm

FEAST is a recently developed eigenvalue algorithm which computes selected interior eigenvalues of real symmetric matrices. It uses contour integral resolvent based projections. A weakness is that the existing algorithm relies on accurate reasoned estimates of the number of eigenvalues within the contour. Examining the singular values of the projections on moderately-sized, randomly-generated test problems motivates orthogonalization-based improvements to the algorithm. The singular value distributions provide experimentally robust estimates of the number of eigenvalues within the contour. The algorithm is modified to handle both Hermitian and general complex matrices. The original algorithm (based on circular contours and Gauss-Legendre quadrature) is extended to contours and quadrature schemes that are recursively subdividable. A general complex recursive algorithm is implemented on rectangular and diamond contours. The accuracy of different quadrature schemes for various contours is investigated.

Compressive strength of interbody cages in the lumbar spine: the effect of cage shape, posterior instrumentation and bone density

One goal of interbody fusion is to increase the height of the degenerated disc space. Interbody cages in particular have been promoted with the claim that they can maintain the disc space better than other methods. There are many factors that can affect the disc height maintenance, including graft or cage design, the quality of the surrounding bone and the presence of supplementary posterior fixation. The present study is an in vitro biomechanical investigation of the compressive behaviour of three different interbody cage designs in a human cadaveric model. The effect of bone density and posterior instrumentation were assessed. Thirty-six lumbar functional spinal units were instrumented with one of three interbody cages: (1) a porous titanium implant with endplate fit (Stratec), (2) a porous, rectangular carbon-fibre implant (Brantigan) and (3) a porous, cylindrical threaded implant (Ray). Posterior instrumentation (USS) was applied to half of the specimens. All specimens were subjected to axial compression displacement until failure. Correlations between both the failure load and the load at 3 mm displacement with the bone density measurements were observed. Neither the cage design nor the presence of posterior instrumentation had a significant effect on the failure load. The loads at 3 mm were slightly less for the Stratec cage, implying lower axial stiffness, but were not different with posterior instrumentation. The large range of observed failure loads overlaps the potential in vivo compressive loads, implying that failure of the bone-implant interface may occur clinically. Preoperative measurements of bone density may be an effective tool to predict settling around interbody cages.

Cell shape-dependent early responses of fibroblasts to cyclic strain

Randomly spread fibroblasts on fibronectin-coated elastomeric membranes respond to cyclic strain by a varying degree of focal adhesion assembly and actin reorganization. We speculated that the individual shape of the cells, which is linked to cytoskeletal structure and pre-stress, might tune these integrin-dependent mechanotransduction events. To this aim, fibronectin circles, squares and rectangles of identical surface area (2000μm(2)) were micro-contact printed onto elastomeric substrates. Fibroblasts plated on these patterns occupied the corresponding shapes. Cyclic 10% equibiaxial strain was applied to patterned cells for 30min, and changes in cytoskeleton and cell-matrix adhesions were quantified after fluorescence staining. After strain, megakaryocytic leukemia-1 protein translocated to the nucleus in most cells, indicating efficient RhoA activation independently of cell shape. However, circular and square cells (with radial symmetry) showed a significantly greater increase in the number of actin stress fibers and vinculin-positive focal adhesions after cyclic strain than rectangular (bipolar) cells of identical size. Conversely, cyclic strain induced larger changes in pY397-FAK positive focal complexes and zyxin relocation from focal adhesions to stress fibers in bipolar compared to symmetric cells. Thus, radially symmetric cells responded to cyclic strain with a larger increase in assembly, whereas bipolar cells reacted with more pronounced reorganization of actin stress fibers and matrix contacts. We conclude that integrin-mediated responses to external mechanical strain are differentially modulated in cells that have the same spreading area but different geometries, and do not only depend on mere cell size.

Reply to the editor

We appreciate the comments and concerns expressed by Arakawa and colleagues regarding our article, titled “Pulsatile control of rotary blood pumps: Does the modulation waveform matter?”1 Unfortunately, we have to disagree with Arakawa and colleagues. As is obvious from the title of our article, it investigates the effect of different waveforms on the heart–device interaction. In contrast to the authors' claim, this is the first article in the literature that uses basic waveforms (sine, triangle, saw tooth, and rectangular) with different phase shifts to examines their impact on left ventricular unloading. The previous publications2, 3 and 4 just varied the pump speed during systole and diastole, which was first reported by Bearnson and associates5 in 1996, and studied its effect on aortic pressure, coronary flow, and end-diastolic volume. We should mention that dp/dtmax is a load-sensitive parameter of contractility and not representative for the degree of unloading. Moreover, none of the aforementioned reports has studied mechanical unloading and in particular the stroke work of the left ventricle. Our method is unique because we do not just alternate between high and low speed but have accurate control of the waveform because of the direct drive system of Levitronix Technologies LLC (Waltham, Mass) and a custom-developed pump controller. Without referring, Arakawa and associates state “several previous studies have already reported the coronary flow diminishes as the left ventricular assist device support increases.” It should be noted that all the waveforms used in our study have 2000 rpm average value with 1000 rpm amplitude, which is not an excessive speed for the CentriMag rotary pump (Levitronix) to collapse the ventricle and diminish the coronary flow. We agree with Arakawa and coworkers that there is a need for a heart failure model to come to more relevant results with respect to clinical expectations. However, we have explored many existing models, including species and breeds that have a native proneness to cardiomyopathy, but all of them differ from the genetic presentation in humans. We certainly do not believe that the use of microembolization, in which the coronary circulation is impaired by the injection of microspheres, would form a good model from which to draw conclusions about coronary flow change under different loading conditions. A model would be needed in which either an infarct is created to mimic ischemic heart failure or the coronary circulation remains untouched to simulate, for instance, dilated cardiomyopathy. Furthermore, in discussion we clearly mention that “lack of heart failure is a major limitation of our study.” We also believe that unloading is not the only factor of the cardiac functional recovery, and an excessive unloading of the left ventricle might lead to cardiac tissue atrophy. Therefore, in our article we mention that control of the level of cardiac unloading by assist devices has been suggested as a mechanical tool to promote recovery, and more studies are required to find better strategies for the speed modulation of rotary pumps and to achieve an optimal heart load control to enhance myocardial recovery. Finally, there are many publications about pulsing rotary blood pumps and it was impossible to include them all. We preferred to reference some of the earlier basic works such as an original research by Bearnson and coworkers5 and another article published by our group,6 which is more relevant.

Continuous Flow Analysis Method for Determination of Dissolved Reactive Phosphorus in Ice Cores

Phosphorus (P) is an essential macronutrient for all living organisms. Phosphorus is often present in nature as the soluble phosphate ion PO43– and has biological, terrestrial, and marine emission sources. Thus PO43– detected in ice cores has the potential to be an important tracer for biological activity in the past. In this study a continuous and highly sensitive absorption method for detection of dissolved reactive phosphorus (DRP) in ice cores has been developed using a molybdate reagent and a 2-m liquid waveguide capillary cell (LWCC). DRP is the soluble form of the nutrient phosphorus, which reacts with molybdate. The method was optimized to meet the low concentrations of DRP in Greenland ice, with a depth resolution of approximately 2 cm and an analytical uncertainty of 1.1 nM (0.1 ppb) PO43–. The method has been applied to segments of a shallow firn core from Northeast Greenland, indicating a mean concentration level of 2.74 nM (0.26 ppb) PO43– for the period 1930–2005 with a standard deviation of 1.37 nM (0.13 ppb) PO43– and values reaching as high as 10.52 nM (1 ppb) PO43–. Similar levels were detected for the period 1771–1823. Based on impurity abundances, dust and biogenic particles were found to be the most likely sources of DRP deposited in Northeast Greenland.

The origin of oriented lakes: evidence from the Bolivian Amazon

The presence of hundreds of rectangular and oriented lakes is one of the most striking characteristics of the Llanos de Moxos (LM) landscape in the Bolivian Amazon. Oriented lakes also occur in the Arctic coastal plains of Russia, Alaska and Canada and along the Atlantic Coastal Plain from northeast Florida to southeast New Jersey and along the coast of northeast Brazil. Many different mechanisms have been proposed for their formation. In the LM, Plafker's (1964) tectonic model, in which subsidence results from the propagation of bedrock faults through the foreland sediments, is the most accepted. However, this model has not been verified. Here, we present new results from stratigraphic transects across the borders of three rectangular and oriented lakes in the LM. A paleosol buried under mid-Holocene sediments is used as a stratigraphic marker to assess the vertical displacement of sediments on both sides of the alleged faults. Our results show that there is no vertical displacement and, therefore, that Plafker's model can be ruled out. We suggest that, among all the proposed mechanisms behind lake formation, the combined action of wind and waves is the most likely. The evidence from the LM provides new hints for the formation of oriented lakes worldwide.

Radiometric Method for Emissivity Retrieval in High Reflective Materials

High reflective materials in the microwave region play a very important role in the realization of antenna reflectors for a broad range of applications, including radiometry. These reflectors have a characteristic emissivity which needs to be characterized accurately in order to perform a correct radiometric calibration of the instrument. Such a characterization can be performed by using open resonators, waveguide cavities or by radiometric measurements. The latter consists of comparative radiometric observations of absorbers, reference mirrors and the sample under test, or using the cold sky radiation as a direct reference source. While the first two mentioned techniques are suitable for the characterization of metal plates and mirrors, the latter has the advantages to be also applicable to soft materials. This paper describes how, through this radiometric techniques, it is possible to characterize the emissivity of the sample relative to a reference mirror and how to characterize the absolute emissivity of the latter by performing measurements at different incident angles. The results presented in this paper are based on our investigations on emissivity of a multilayer insulation material (MLI) for space mission, at the frequencies of 22 and 90 GHz.

The timing of significant arch dimensional changes with fixed orthodontic appliances: data from a multicenter randomised controlled trial.

OBJECTIVES To identify the timing of significant arch dimensional increases during orthodontic alignment involving round and rectangular nickel-titanium (NiTi) wires and rectangular stainless steel (SS). A secondary aim was to compare the timing of changes occurring with conventional and self-ligating fixed appliance systems. METHODS In this non-primary publication, additional data from a multicenter randomised trial initially involving 96 patients, aged 16 years and above, were analysed. The main pre-specified outcome measures were the magnitude and timing of maxillary intercanine, interpremolar, and intermolar dimensions. Each participant underwent alignment with a standard Damon (Ormco, Orange, CA) wire sequence for a minimum of 34 weeks. Blinding of clinicians and patients was not possible; however, outcome assessors and data analysts were kept blind to the appliance type during data analysis. RESULTS Complete data were obtained from 71 subjects. Significant arch dimensional changes were observed relatively early in treatment. In particular, changes in maxillary inter-first and second premolar dimensions occurred after alignment with an 0.014in. NiTi wire (P<0.05). No statistical differences in transverse dimensions were found between rectangular NiTi and working SS wires for each transverse dimension (P>0.05). Bracket type had no significant effect on the timing of the transverse dimensional changes. CONCLUSIONS Arch dimensional changes were found to occur relatively early in treatment, irrespective of the appliance type. Nickel-titanium wires may have a more profound effect on transverse dimensions than previously believed. CLINICAL SIGNIFICANCE On the basis of this research orthodontic expansion may occur relatively early in treatment. Nickel-titanium wires may have a more profound effect on transverse dimensions than previously believed.

Insight into the dynamic ligand exchange process involved in bipyridyl linked arene ruthenium metalla-rectangles

The dynamic ligand exchange behavior of cationic arene ruthenium metalla-rectangles of the type [(pcymene) 4Ru4(OOXOO)2(NXN)2]4+ (OOXOO ¼ oxalato, 2,5-dioxydo-1,4-benzoquinonato, 5,8-dioxydo-1,4-naphthoquinonato; NXN ¼ 4,40-bipyridine-H8, 4,40-bipyridine-D8) has been studied in solution. The robustness of the rectangular architecture has been evidenced by NMR and ESI mass spectrometry. Thermodynamic and kinetic aspects of the ligand exchange process have been explored using 1H/2D isotope labeling of the 4,40-bipyridine connectors. This study shows that ligand exchange does not proceed spontaneously for these metalla-assemblies, even at high temperature, unless an external stimulus is applied.

Numerical study of Anderson localization of terahertz waves in disordered waveguides

We present a numerical study of electromagnetic wave transport in disordered quasi-one-dimensional waveguides at terahertz frequencies. Finite element method calculations of terahertz wave propagation within LiNbO3 waveguides with randomly arranged air-filled circular scatterers exhibit an onset of Anderson localization at experimentally accessible length scales. Results for the average transmission as a function of waveguide length and scatterer density demonstrate a clear crossover from diffusive to localized transport regime. In addition, we find that transmission fluctuations grow dramatically when crossing into the localized regime. Our numerical results are in good quantitative agreement with theory over a wide range of experimentally accessible parameters both in the diffusive and localized regime opening the path towards experimental observation of terahertz wave localization.

Detecting single-target changes in multiple object tracking: The case of peripheral vision.

In sports games, it is often necessary to perceive a large number of moving objects (e.g., the ball and players). In this context, the role of peripheral vision for processing motion information in the periphery is often discussed especially when motor responses are required. In an attempt to test the basal functionality of peripheral vision in those sports-games situations, a Multiple Object Tracking (MOT) task that requires to track a certain number of targets amidst distractors, was chosen. Participants’ primary task was to recall four targets (out of 10 rectangular stimuli) after six seconds of quasi-random motion. As a second task, a button had to be pressed if a target change occurred (Exp 1: stop vs. form change to a diamond for 0.5 s; Exp 2: stop vs. slowdown for 0.5 s). While eccentricities of changes (5-10° vs. 15-20°) were manipulated, decision accuracy (recall and button press correct), motor response time as well as saccadic reaction time were calculated as dependent variables. Results show that participants indeed used peripheral vision to detect changes, because either no or very late saccades to the changed target were executed in correct trials. Moreover, a saccade was more often executed when eccentricities were small. Response accuracies were higher and response times were lower in the stop conditions of both experiments while larger eccentricities led to higher response times in all conditions. Summing up, it could be shown that monitoring targets and detecting changes can be processed by peripheral vision only and that a monitoring strategy on the basis of peripheral vision may be the optimal one as saccades may be afflicted with certain costs. Further research is planned to address the question whether this functionality is also evident in sports tasks.

Detecting single-target changes in multiple object tracking: The case of peripheral vision

In sports games, it is often necessary to perceive a large number of moving objects (e.g., the ball and players). In this context, the role of peripheral vision for processing motion information in the periphery is often discussed especially when motor responses are required. In an attempt to test the capability of using peripheral vision in those sports-games situations, a Multiple-Object-Tracking task that requires to track a certain number of targets amidst distractors, was chosen to determine the sensitivity of detecting target changes with peripheral vision only. Participants’ primary task was to recall four targets (out of 10 rectangular stimuli) after six seconds of quasi-random motion. As a second task, a button had to be pressed if a target change occurred (Exp 1: stop vs. form change to a diamond for 0.5 s; Exp 2: stop vs. slowdown for 0.5 s). Eccentricities of changes (5-10° vs. 15-20°) were manipulated, decision accuracy (recall and button press correct), motor response time and saccadic reaction time (change onset to saccade onset) were calculated and eye-movements were recorded. Results show that participants indeed used peripheral vision to detect changes, because either no or very late saccades to the changed target were executed in correct trials. Moreover, a saccade was more often executed when eccentricities were small. Response accuracies were higher and response times were lower in the stop conditions of both experiments while larger eccentricities led to higher response times in all conditions. Summing up, it could be shown that monitoring targets and detecting changes can be processed by peripheral vision only and that a monitoring strategy on the basis of peripheral vision may be the optimal one as saccades may be afflicted with certain costs. Further research is planned to address the question whether this functionality is also evident in sports tasks.

Effects of Strand Geometry on Selected Properties of Long-Stand Structural Composite Lumber Made from Northeastern Hardwoods

Structural composite lumber (SCL) products often possess significantly higher design values than the top grades of solid lumber, making it a popular choice for both residential and commercial applications. The enhanced mechanical properties of SCL are mainly due to defect randomization and densification of the wood fiber, both largely functions of the size, shape and composition (species) of the wood element. Traditionally, SCL manufacturers have used thin, rectangular elements produced from either moderate density softwoods or low density hardwoods. Higher density hardwood species have been avoided, as they require higher pressures to adequately densify and consolidate the wood furnish. These higher pressures can lead to increased manufacturing costs, damage to the wood fiber and/or a product that is too dense, making it heavy and unreceptive to common mechanical fastening techniques. In the northeastern United States high density, diffuse-porous hardwoods (such as maple, beech and birch) are abundant. Use of these species as primary furnish for a SCL product may allow for a competitive advantage in terms of resource cost against products that rely on veneer grade logs. Proximity to this abundant and relatively inexpensive resource may facilitate entry of SCL production facilities in the northeastern United States, where currently none exist. However, modifications to current strand sizes, geometries or production techniques will likely be required to allow for use of these species. A new SCL product concept has been invented allowing for use of these high density hardwoods. The product, referred to as long-strand structural composite lumber (LSSCL), uses strands of significantly larger cross sectional areas and volumes than existing SCL products. In spite of the large strand size, satisfactory consolidation is achieved without excessive densification of the wood fiber through use of a symmetrical strand geometric cross-section. LSSCL density is similar to that of existing SCL products, but is due mainly to the inherent density of the species, rather than through densification. An experiment was designed and conducted producing LSSCL from both large (7/16”) and small (1/4”) strands, of both square and triangular geometric cross sections. Testing results indicate that the large, triangular strands produce LSSCL beams with projected design values of: Modulus of elasticity (MOEapp) – 1,750,000 psi; Allowable bending stress (Fb) – 2750 psi; Allowable shear stress (Fv) – 260 psi. Several modifications are recommended which may lead to improvement of these values, likely allowing for competition against existing SCL products.

Equity Markets, the Money Market, and Long-Run Monetary Neutrality

This paper outlines a process for teaching long-run neutrality of money, drawing an analogy between equity markets and the money market. The key points in the discussion include the following: (1) What is the price of money? (2) Why does the long-run demand for money trace out a rectangular hyperbola? (3) Why does the slow adjustment of goods and service prices to changes in the stock of money lead to a different short-run demand for money? and (4) Why does a successful currency reform generate similar short-run movements in the price of money as movements in equity share prices after a change in the supply of shares? I have used this approach successfully for over 30 years at all levels, wherever I need to discuss the money market in a macroeconomic model.