Theoretical and experimental studies on the electronic, optical, and structural properties of poly-pyrrole-2-carboxylic acid films

A theoretical approach is used here to explain experimental results obtained from the electrosynthesis of polypyrrole-2-carboxylic acid (PPY-2-COOH) films in nonaqueous medium. An analysis of the Fukui function (reactivity index) indicates that the monomer (pyrrole-2-carboxylic acid, PY-2-COOH), and dimers and trimers are oxidized in the C4 or C5 positions of the heterocyclic ring of the PY-2-COOH structure. After calculating the heat of formation using semiempirical Austin Model 1 post-Hartree-Fock parameterization for dimer species, both C4 and C5 positions adjacent to the aromatic rings of PPY-2-COOH were considered the most susceptible ones to oxidative coupling reactions. The ZINDO-S/CI semiempirical method was used to simulate the electronic transitions typically seen in the UV-VIS-NIR range in monomer and oligomers with different conjugation lengths. The use of an electrochemical quartz crystal microbalance provides sufficient information to propose a polymerization mechanism of PY-2-COOH based on molecular modeling and experimental results.

Kernel Methods for Tree Structured Data

Machine learning comprises a series of techniques for automatic extraction of meaningful information from large collections of noisy data. In many real world applications, data is naturally represented in structured form. Since traditional methods in machine learning deal with vectorial information, they require an a priori form of preprocessing. Among all the learning techniques for dealing with structured data, kernel methods are recognized to have a strong theoretical background and to be effective approaches. They do not require an explicit vectorial representation of the data in terms of features, but rely on a measure of similarity between any pair of objects of a domain, the kernel function. Designing fast and good kernel functions is a challenging problem. In the case of tree structured data two issues become relevant: kernel for trees should not be sparse and should be fast to compute. The sparsity problem arises when, given a dataset and a kernel function, most structures of the dataset are completely dissimilar to one another. In those cases the classifier has too few information for making correct predictions on unseen data. In fact, it tends to produce a discriminating function behaving as the nearest neighbour rule. Sparsity is likely to arise for some standard tree kernel functions, such as the subtree and subset tree kernel, when they are applied to datasets with node labels belonging to a large domain. A second drawback of using tree kernels is the time complexity required both in learning and classification phases. Such a complexity can sometimes prevents the kernel application in scenarios involving large amount of data. This thesis proposes three contributions for resolving the above issues of kernel for trees. A first contribution aims at creating kernel functions which adapt to the statistical properties of the dataset, thus reducing its sparsity with respect to traditional tree kernel functions. Specifically, we propose to encode the input trees by an algorithm able to project the data onto a lower dimensional space with the property that similar structures are mapped similarly. By building kernel functions on the lower dimensional representation, we are able to perform inexact matchings between different inputs in the original space. A second contribution is the proposal of a novel kernel function based on the convolution kernel framework. Convolution kernel measures the similarity of two objects in terms of the similarities of their subparts. Most convolution kernels are based on counting the number of shared substructures, partially discarding information about their position in the original structure. The kernel function we propose is, instead, especially focused on this aspect. A third contribution is devoted at reducing the computational burden related to the calculation of a kernel function between a tree and a forest of trees, which is a typical operation in the classification phase and, for some algorithms, also in the learning phase. We propose a general methodology applicable to convolution kernels. Moreover, we show an instantiation of our technique when kernels such as the subtree and subset tree kernels are employed. In those cases, Direct Acyclic Graphs can be used to compactly represent shared substructures in different trees, thus reducing the computational burden and storage requirements.

Salt stress responses in pear and quince: physiological and molecular aspects

The productivity of agricultural crops is seriously limited by salinity. This problem is rapidly increasing, particularly in irrigated lands. Like almost all the fruit tree species, Pyrus communis is generally considered a salt sensitive species, but only little information is available on its behavior under saline conditions. Previous studies, carried out in the Department of Fruit Tree and Woody Plant Science (University of Bologna), focused their attention on pear and quince salt stress responses to understand which rootstock would be the most suitable for pear in order to tolerate a salt stress condition. It has been reported that pear and quince have different ability in the uptake, translocation and accumulation of chloride (Cl-) and sodium (Na+) ions, when plants were irrigated for one season with saline water (5 dS/m). The aim of the present work was to deepen these aspects and investigate salt stress responses in pear and quince. Two different experiments have been performed: a “short-term” trial in a growth chamber and a “long-term” experiment in the open field. In the short-term experiment, three different genotypes usually adopted as pear rootstocks (MC, BA29 and Farold®40) and the pear variety Abbé Fétel own rooted have been compared under salt stress conditions. The trial was performed in a hydroponic culture system, applying a 90 mM NaCl stress to half of the plants, after five weeks of normal growth in Hoagland’s solution. During the three-weeks of salt stress treatment, physiological, mineral and molecular analyses were performed in order to monitor, for each genotype, the development of the salt stress responses in comparison with the corresponding “unstressed” plants. Farold®40 and Abbé Fétel own rooted showed the onset of leaf necrosis, due to salt toxicity, one week before quinces. Moreover, quinces displayed a significant delay in premature senescence of old leaves, while pears emerged for their ability to regenerate new leaves from apparently dead foliage with the salt stress still running. Physiological measurements, such as shoots length, chlorophyll (Chl) content, and photosynthesis, have been carried out and revealed that pears exhibited a significant reduction in water content and a wilting aspect, while for quinces a decrease in Chl content and a growth slowdown were observed. At the end of the trial, all plants were collected and organs separated for dry weight estimation and mineral analyses (Cu, Fe, Mn, Zn Mg, Ca, K, Na and Cl). Mineral contents have been affected by salinity; same macro/micro nutrients were altered in some organs or relocated within the plant. This plant response could have partially contributed to face the salt stress. Leaves and roots have been harvested for molecular analyses at four different times during stress conditions. Molecular analyses consisted of the gene expression study of three main ion transporters, well known in Arabidopsis thaliana as salt-tolerance determinants in the “SOS” pathway: NHX1 (tonoplast Na+/H+ antiporter), SOS1 (plasmalemma Na+/H+ antiporter) and HKT1 (K+ high-affinity and Na+ low-affinity transporter). These studies showed that two quince rootstocks adopted different responsive mechanisms to NaCl stress. BA29 increased its Na+ sequestration activity into leaf vacuoles, while MC enhanced temporarily the same ability, but in roots. Farold®40, instead, exhibited increases in SOS1 and HKT1 expression mainly at leaf level in the attempt to retrieve Na+ from xylem, while Abbé Fétel differently altered the expression of these genes in roots. Finally, each genotype showed a peculiar response to salt stress that was the sum of its ability in Na+ exclusion, osmotic tolerance and tissue tolerance. In the long-term experiment, potted trees of the pear variety Abbé Fétel grafted on different rootstocks (MC, BA29 and Farold®40), or own rooted and also rootstocks only were subjected to a salt stress through saline water irrigation with an electrical conductivity of 5 dS/m for two years. The purposes of this study were to evaluate salinity effects on physiological (shoot length, number of buds, photosynthesis, etc.) and yield parameters of cultivar Abbé Fétel in the different combinations and to determine the salt amount that pear is able to tolerate over the years. With this work, we confirmed the previous hypothesis that pear, despite being classified as a salt-sensitive fruit tree, can be cultivated for two years under saline water irrigation, without showing any salt toxicity symptoms or severe drawbacks on plant development and production. Among different combinations, Abbé Fétel grafted on MC resulted interesting for its peculiar behaviors under salt stress conditions. In the near future, further investigations on physiological and molecular aspects will be necessary to enrich and broaden the knowledge of salt stress responses in pear.

Electrothermal characterization and nonlinear modelling of GaN transistors for telecommunication circuit design

This thesis starts showing the main characteristics and application fields of the AlGaN/GaN HEMT technology, focusing on reliability aspects essentially due to the presence of low frequency dispersive phenomena which limit in several ways the microwave performance of this kind of devices. Based on an equivalent voltage approach, a new low frequency device model is presented where the dynamic nonlinearity of the trapping effect is taken into account for the first time allowing considerable improvements in the prediction of very important quantities for the design of power amplifier such as power added efficiency, dissipated power and internal device temperature. An innovative and low-cost measurement setup for the characterization of the device under low-frequency large-amplitude sinusoidal excitation is also presented. This setup allows the identification of the new low frequency model through suitable procedures explained in detail. In this thesis a new non-invasive empirical method for compact electrothermal modeling and thermal resistance extraction is also described. The new contribution of the proposed approach concerns the non linear dependence of the channel temperature on the dissipated power. This is very important for GaN devices since they are capable of operating at relatively high temperatures with high power densities and the dependence of the thermal resistance on the temperature is quite relevant. Finally a novel method for the device thermal simulation is investigated: based on the analytical solution of the tree-dimensional heat equation, a Visual Basic program has been developed to estimate, in real time, the temperature distribution on the hottest surface of planar multilayer structures. The developed solver is particularly useful for peak temperature estimation at the design stage when critical decisions about circuit design and packaging have to be made. It facilitates the layout optimization and reliability improvement, allowing the correct choice of the device geometry and configuration to achieve the best possible thermal performance.

Logarithmically smooth deformations of strict normal crossing logarithmically symplectic varieties

In this thesis we give a definition of the term logarithmically symplectic variety; to be precise, we distinguish even two types of such varieties. The general type is a triple $(f,nabla,omega)$ comprising a log smooth morphism $fcolon Xtomathrm{Spec}kappa$ of log schemes together with a flat log connection $nablacolon LtoOmega^1_fotimes L$ and a ($nabla$-closed) log symplectic form $omegainGamma(X,Omega^2_fotimes L)$. We define the functor of log Artin rings of log smooth deformations of such varieties $(f,nabla,omega)$ and calculate its obstruction theory, which turns out to be given by the vector spaces $H^i(X,B^bullet_{(f,nabla)}(omega))$, $i=0,1,2$. Here $B^bullet_{(f,nabla)}(omega)$ is the class of a certain complex of $mathcal{O}_X$-modules in the derived category $mathrm{D}(X/kappa)$ associated to the log symplectic form $omega$. The main results state that under certain conditions a log symplectic variety can, by a flat deformation, be smoothed to a symplectic variety in the usual sense. This may provide a new approach to the construction of new examples of irreducible symplectic manifolds.

A VEGETATION RECONSTRUCTION OF UNION COUNTY, PENNSYLVANIA DURING THE EARLY HISTORICAL PERIOD (1755-1855)

This thesis assesses relationships between vegetation and topography and the impact of human tree-cutting on the vegetation of Union County during the early historical era (1755-1855). I use early warrant maps and forestry maps from the Pennsylvania historical archives and a warrantee map from the Union County courthouse depicting the distribution of witness trees and non-tree surveyed markers (posts and stones) in early European settlement land surveys to reconstruct the vegetation and compare vegetation by broad scale (mountains and valleys) and local scale (topographic classes with mountains and valleys) topography. I calculated marker density based on 2 km x 2 km grid cells to assess tree-cutting impacts. Valleys were mostly forests dominated by white oak (Quercus alba) with abundant hickory (Carya spp.), pine (Pinus spp.), and black oak (Quercus velutina), while pine dominated what were mostly pine-oak forests in the mountains. Within the valleys, pine was strongly associated with hilltops, eastern hemlock (Tsuga canadensis) was abundant on north slopes, hickory was associated with south slopes, and riparian zones had high frequencies of ash (Fraxinus spp.) and hickory. In the mountains, white oak was infrequent on south slopes, chestnut (Castanea dentata) was more abundant on south slopes and ridgetops than north slopes and mountain coves, and white oak and maple (Acer spp.) were common in riparian zones. Marker density analysis suggests that trees were still common over most of the landscape by 1855. The findings suggest there were large differences in vegetation between valleys and mountains due in part to differences in elevation, and vegetation differed more by topographic classes in the valleys than in the mountains. Possible areas of tree-cutting were evenly distributed by topographic classes, suggesting Europeans settlers were clearing land and harvesting timber in most areas of Union County.

X-ray phase contrast imaging of metallized and non-metallized biological samples

Imaging of biological samples has been performed with a variety of techniques for example electromagnetic waves, electrons, neutrons, ultrasound and X-rays. Also conventional X-ray imaging represents the basis of medical diagnostic imaging, it remains of limited use in this application because it is based solely on the differential absorption of X-rays by tissues. Coherent and bright photon beams, such as those produced by third-generation synchrotron X-ray sources, provide further information on subtle X-ray phase changes at matter interfaces. This complements conventional X-ray absorption by edge enhancement phenomena. Thus, phase contrast imaging has the potential to improve the detection of structures on images by detecting those structures that are invisible with X-ray absorption imaging. Images of a weakly absorbing nylon fibre were recorded in in-line holography geometry using a high resolution low-noise CCD camera at the ESRF in Grenoble. The method was also applied to improve image contrast for images of biological tissues. This paper presents phase contrast microradiographs of vascular tree casts and images of a housefly. These reveal very fine structures, that remain invisible with conventional absorption contrast only.

Intra-annual radial growth and water relations of trees - implications towards a growth mechanism

There is a missing link between tree physiological and wood-anatomical knowledge which makes it impossible mechanistically to explain and predict the radial growth of individual trees from climate data. Empirical data of microclimatic factors, intra-annual growth rates, and tree-specific ratios between actual and potential transpiration (T PET−1) of trees of three species (Quercus pubescens, Pinus sylvestris, and Picea abies) at two dry sites in the central Wallis, Switzerland, were recorded from 2002 to 2004 at a 10 min resolution. This included the exceptionally hot and dry summer of 2003. These data were analysed in terms of direct (current conditions) and indirect impacts (predispositions of the past year) on growth. Rain was found to be the only factor which, to a large extent, consistently explained the radial increment for all three tree species at both sites and in the short term as well. Other factors had some explanatory power on the seasonal time-scale only. Quercus pubescens built up much of its tree ring before bud break. Pinus sylvestris and Picea abies started radial growth 1–2 weeks after Quercus pubescens and this was despite the fact that they had a high T PET−1 before budburst and radial growth started. A high T PET−1 was assumed to be related to open stomata, a very high net CO2 assimilation rate, and thus a potential carbon (C)-income for the tree. The main period of radial growth covered about 30–70% of the productive days of a year. In terms of C-allocation, these results mean that Quercus pubescens depended entirely on internal C-stores in the early phase of radial growth and that for all three species there was a long time period of C-assimilation which was not used for radial growth in above-ground wood. The results further suggest a strong dependence of radial growth on the current tree water relations and only secondarily on the C-balance. A concept is discussed which links radial growth over a feedback loop to actual tree water-relations and long-term affected C-storage to microclimate.

Assessing success of forest restoration efforts in degraded montane cloud forests in southern Mexico

Montane cloud forests are home to great biodiversity. However, non-sustainable anthropogenic activities have led to the loss of forest cover in southern Mexico. Increasing conservation, restoration and sustainable use of forest resources prevents the loss of cloud forests. In this study, success of forest restoration was evaluated in a degraded forest of Highlands Chiapas. The goal of this study was to assess the structure and composition of native tree species. We evaluated vegetation composition at three sites that had undergone enrichment plantings. Floristic composition and structure of the herbaceous, seedling, sapling, and overstory layers were measured. A total of sixty-six native tree species were recorded. Enrichment planting was found to have increased tree diversity. Moreover, 54% of the planted species were found in the understory, indicating that they were successfully recruiting. In conclusion, enrichment planting can aid in the conservation of forest cover in degraded areas.

Towards greater accuracy in individual-tree mortality regression

Background mortality is an essential component of any forest growth and yield model. Forecasts of mortality contribute largely to the variability and accuracy of model predictions at the tree, stand and forest level. In the present study, I implement and evaluate state-of-the-art techniques to increase the accuracy of individual tree mortality models, similar to those used in many of the current variants of the Forest Vegetation Simulator, using data from North Idaho and Montana. The first technique addresses methods to correct for bias induced by measurement error typically present in competition variables. The second implements survival regression and evaluates its performance against the traditional logistic regression approach. I selected the regression calibration (RC) algorithm as a good candidate for addressing the measurement error problem. Two logistic regression models for each species were fitted, one ignoring the measurement error, which is the “naïve” approach, and the other applying RC. The models fitted with RC outperformed the naïve models in terms of discrimination when the competition variable was found to be statistically significant. The effect of RC was more obvious where measurement error variance was large and for more shade-intolerant species. The process of model fitting and variable selection revealed that past emphasis on DBH as a predictor variable for mortality, while producing models with strong metrics of fit, may make models less generalizable. The evaluation of the error variance estimator developed by Stage and Wykoff (1998), and core to the implementation of RC, in different spatial patterns and diameter distributions, revealed that the Stage and Wykoff estimate notably overestimated the true variance in all simulated stands, but those that are clustered. Results show a systematic bias even when all the assumptions made by the authors are guaranteed. I argue that this is the result of the Poisson-based estimate ignoring the overlapping area of potential plots around a tree. Effects, especially in the application phase, of the variance estimate justify suggested future efforts of improving the accuracy of the variance estimate. The second technique implemented and evaluated is a survival regression model that accounts for the time dependent nature of variables, such as diameter and competition variables, and the interval-censored nature of data collected from remeasured plots. The performance of the model is compared with the traditional logistic regression model as a tool to predict individual tree mortality. Validation of both approaches shows that the survival regression approach discriminates better between dead and alive trees for all species. In conclusion, I showed that the proposed techniques do increase the accuracy of individual tree mortality models, and are a promising first step towards the next generation of background mortality models. I have also identified the next steps to undertake in order to advance mortality models further.

Automated design of self-stabilization

Self-stabilization is a property of a distributed system such that, regardless of the legitimacy of its current state, the system behavior shall eventually reach a legitimate state and shall remain legitimate thereafter. The elegance of self-stabilization stems from the fact that it distinguishes distributed systems by a strong fault tolerance property against arbitrary state perturbations. The difficulty of designing and reasoning about self-stabilization has been witnessed by many researchers; most of the existing techniques for the verification and design of self-stabilization are either brute-force, or adopt manual approaches non-amenable to automation. In this dissertation, we first investigate the possibility of automatically designing self-stabilization through global state space exploration. In particular, we develop a set of heuristics for automating the addition of recovery actions to distributed protocols on various network topologies. Our heuristics equally exploit the computational power of a single workstation and the available parallelism on computer clusters. We obtain existing and new stabilizing solutions for classical protocols like maximal matching, ring coloring, mutual exclusion, leader election and agreement. Second, we consider a foundation for local reasoning about self-stabilization; i.e., study the global behavior of the distributed system by exploring the state space of just one of its components. It turns out that local reasoning about deadlocks and livelocks is possible for an interesting class of protocols whose proof of stabilization is otherwise complex. In particular, we provide necessary and sufficient conditions – verifiable in the local state space of every process – for global deadlock- and livelock-freedom of protocols on ring topologies. Local reasoning potentially circumvents two fundamental problems that complicate the automated design and verification of distributed protocols: (1) state explosion and (2) partial state information. Moreover, local proofs of convergence are independent of the number of processes in the network, thereby enabling our assertions about deadlocks and livelocks to apply on rings of arbitrary sizes without worrying about state explosion.

Small-scale timber stand management techniques : a case study of woodlots in Isangati, Tanzania

Small-scale village woodlots of less than 0.5ha are the preferred use of land for local farmers with extra land in the village of Isangati, a small community located in the southern highlands of Tanzania. Farmers view woodlots as lucrative investments that do not involve intensive labor or time. The climate is ideal for the types of trees grown and the risks are minimal with no serious threats from insects, fires, thieves, or grazing livestock. It was hypothesized that small-scale village woodlot owners were not maximizing timber outputs with their current timber stand management and harvesting techniques. Personal interviews were conducted over a five month period and field data was collected at each farmer’s woodlots over a seven month period. Woodlot field data included woodlot size, number of trees, tree species, tree height, dbh, age, and spacing. The results indicated that the lack of proper woodlot management techniques results in failure to fully capitalize on the investment of woodlots. While farmers should continue with their current harvesting rotations, some of the reasons for not maximizing tree growth include close spacing (2m x 2m), no tree thinning, extreme pruning (60% of tree), and little to no weeding. Through education and hands-on woodlot management workshops, the farmers could increase their timber output and value of woodlots.

Development of an optimization model for biofuel facility size and location and a simulation model for design of a biofuel supply chain

To mitigate greenhouse gas (GHG) emissions and reduce U.S. dependence on imported oil, the United States (U.S.) is pursuing several options to create biofuels from renewable woody biomass (hereafter referred to as “biomass”). Because of the distributed nature of biomass feedstock, the cost and complexity of biomass recovery operations has significant challenges that hinder increased biomass utilization for energy production. To facilitate the exploration of a wide variety of conditions that promise profitable biomass utilization and tapping unused forest residues, it is proposed to develop biofuel supply chain models based on optimization and simulation approaches. The biofuel supply chain is structured around four components: biofuel facility locations and sizes, biomass harvesting/forwarding, transportation, and storage. A Geographic Information System (GIS) based approach is proposed as a first step for selecting potential facility locations for biofuel production from forest biomass based on a set of evaluation criteria, such as accessibility to biomass, railway/road transportation network, water body and workforce. The development of optimization and simulation models is also proposed. The results of the models will be used to determine (1) the number, location, and size of the biofuel facilities, and (2) the amounts of biomass to be transported between the harvesting areas and the biofuel facilities over a 20-year timeframe. The multi-criteria objective is to minimize the weighted sum of the delivered feedstock cost, energy consumption, and GHG emissions simultaneously. Finally, a series of sensitivity analyses will be conducted to identify the sensitivity of the decisions, such as the optimal site selected for the biofuel facility, to changes in influential parameters, such as biomass availability and transportation fuel price. Intellectual Merit The proposed research will facilitate the exploration of a wide variety of conditions that promise profitable biomass utilization in the renewable biofuel industry. The GIS-based facility location analysis considers a series of factors which have not been considered simultaneously in previous research. Location analysis is critical to the financial success of producing biofuel. The modeling of woody biomass supply chains using both optimization and simulation, combing with the GIS-based approach as a precursor, have not been done to date. The optimization and simulation models can help to ensure the economic and environmental viability and sustainability of the entire biofuel supply chain at both the strategic design level and the operational planning level. Broader Impacts The proposed models for biorefineries can be applied to other types of manufacturing or processing operations using biomass. This is because the biomass feedstock supply chain is similar, if not the same, for biorefineries, biomass fired or co-fired power plants, or torrefaction/pelletization operations. Additionally, the research results of this research will continue to be disseminated internationally through publications in journals, such as Biomass and Bioenergy, and Renewable Energy, and presentations at conferences, such as the 2011 Industrial Engineering Research Conference. For example, part of the research work related to biofuel facility identification has been published: Zhang, Johnson and Sutherland [2011] (see Appendix A). There will also be opportunities for the Michigan Tech campus community to learn about the research through the Sustainable Future Institute.

Life cycle assessment of biofuels produced by the new integrated hydropyrolysis-hydroconversion (IH2) process

Biofuels are alternative fuels that have the promise of reducing reliance on imported fossil fuels and decreasing emission of greenhouse gases from energy consumption. This thesis analyses the environmental impacts focusing on the greenhouse gas (GHG) emissions associated with the production and delivery of biofuel using the new Integrated Hydropyrolysis and Hydroconversion (IH2) process. The IH2 process is an innovative process for the conversion of woody biomass into hydrocarbon liquid transportation fuels in the range of gasoline and diesel. A cradle-to-grave life cycle assessment (LCA) was used to calculate the greenhouse gas emissions associated with diverse feedstocks production systems and delivery to the IH2 facility plus producing and using these new renewable liquid fuels. The biomass feedstocks analyzed include algae (microalgae), bagasse from a sugar cane-producing locations such as Brazil or extreme southern US, corn stover from Midwest US locations, and forest feedstocks from a northern Wisconsin location. The life cycle greenhouse gas (GHG) emissions savings of 58%–98% were calculated for IH2 gasoline and diesel production and combustion use in vehicles compared to fossil fuels. The range of savings is due to different biomass feedstocks and transportation modes and distances. Different scenarios were conducted to understand the uncertainties in certain input data to the LCA model, particularly in the feedstock production section, the IH2 biofuel production section, and transportation sections.

PHYLOGENETIC ANALYSES, GENE DOSAGE EFFECT AND SELECTION IN MEMBERS OF THE LATERAL BOUNDARY DOMAIN (LBD) GENES TRANSCRIPTION FACTOR FAMILY

Wood formation is an economically and environmentally important process and has played a significant role in the evolution of terrestrial plants. Despite its significance, the molecular underpinnings of the process are still poorly understood. We have previously shown that four Lateral Boundary Domain (LBD) transcription factors have important roles in the regulation of wood formation with two (LBD1 and LBD4) involved in secondary phloem and ray cell development and two (LBD15 and LBD18) in secondary xylem formation. Here, we used comparative phylogenetic analyses to test potential roles of the four LBD genes in the evolution of woodiness. We studied the copy number and variation in DNA and amino acid sequences of the four LBDs in a wide range of woody and herbaceous plant taxa with fully sequenced and annotated genomes. LBD1 showed the highest gene copy number across the studied species, and LBD1 gene copy number was strongly and significantly correlated with the level of ray seriation. The lianas, cucumber and grape, with multiseriate ray cells showed the highest gene copy number (12 and 11, respectively). Because lianas’ growth habit requires significant twisting and bending, the less lignified ray parenchyma cells likely facilitate stem flexibility and maintenance of xylem conductivity. We further demonstrate conservation of amino acids in the LBD18 protein sequences that are specific to woody taxa. Neutrality tests showed evidence for strong purifying selection on these gene regions across various orders, indicating adaptive convergent evolution of LBD18. Structural modeling demonstrates that the conserved amino acids have a significant impact on the tertiary protein structure and thus are likely of significant functional importance.