A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Social value orientation as a moral intuition: Decision-making in the dictator game

We studied the decision making process in the Dictator Game and showed that decisions are the result of a two-step process. In a first step, decision makers generate an automatic, intuitive proposal. Given sufficient motivation and cognitive resources, they adjust this in a second, more deliberated phase. In line with the social intuitionist model, we show that one s Social Value Orientation determines intuitive choice tendencies in the first step, and that this effect is mediated by the dictator s perceived interpersonal closeness with the receiver. Self-interested concerns subsequently leadto a reduction of donation size in step 2. Finally, we show that increasing interpersonal closeness can promote pro-social decision-making.

Simultaneous identification of multiple mammalian species from mixed forensic samples based on mtDNA control region length polymorphism

Molecular species identification in mixed or contaminated biological material has always been problematic. We developed a simple and accurate method for mammal DNA identification in mixtures, based on interspecific mitochondrial DNA control region length polymorphism. Contrary to other published methods dealing with species mixtures, our protocol requires a single universal primer pair and amplification step, and is not based on a pre-defined panel of species. This protocol has been routinely employed by our laboratory for species identification in dozens of human and animal forensic caseworks. Six representative forensic caseworks involving the specific identification of mixed animal samples are reported in this paper, in order to demonstrate the applicability and usefulness of the method.

Induced aggregation operators in decision making with the Dempster-Shafer belief structure

We study the induced aggregation operators. The analysis begins with a revision of some basic concepts such as the induced ordered weighted averaging (IOWA) operator and the induced ordered weighted geometric (IOWG) operator. We then analyze the problem of decision making with Dempster-Shafer theory of evidence. We suggest the use of induced aggregation operators in decision making with Dempster-Shafer theory. We focus on the aggregation step and examine some of its main properties, including the distinction between descending and ascending orders and different families of induced operators. Finally, we present an illustrative example in which the results obtained using different types of aggregation operators can be seen.

Induced aggregation operators in decision making with the Dempster-Shafer belief structure

We study the induced aggregation operators. The analysis begins with a revision of some basic concepts such as the induced ordered weighted averaging (IOWA) operator and the induced ordered weighted geometric (IOWG) operator. We then analyze the problem of decision making with Dempster-Shafer theory of evidence. We suggest the use of induced aggregation operators in decision making with Dempster-Shafer theory. We focus on the aggregation step and examine some of its main properties, including the distinction between descending and ascending orders and different families of induced operators. Finally, we present an illustrative example in which the results obtained using different types of aggregation operators can be seen.

Impacts on the deep-sea ecosystem by a severe coastal storm

Major coastal storms, associated with strong winds, high waves and intensified currents, and occasionally with heavy rains and flash floods, are mostly known because of the serious damage they can cause along the shoreline and the threats they pose to navigation. However, there is a profound lack of knowledge on the deep-sea impacts of severe coastal storms. Concurrent measurements of key parameters along the coast and in the deep-sea are extremely rare. Here we present a unique data set showing how one of the most extreme coastal storms of the last decades lashing the Western Mediterranean Sea rapidly impacted the deep-sea ecosystem. The storm peaked the 26th of December 2008 leading to the remobilization of a shallow-water reservoir of marine organic carbon associated with fine particles and resulting in its redistribution across the deep basin. The storm also initiated the movement of large amounts of coarse shelf sediment, which abraded and buried benthic communities. Our findings demonstrate, first, that severe coastal storms are highly efficient in transporting organic carbon from shallow water to deep water, thus contributing to its sequestration and, second, that natural, intermittent atmospheric drivers sensitive to global climate change have the potential to tremendously impact the largest and least known ecosystem on Earth, the deep-sea ecosystem.

A Monte Carlo-based procedure for independent monitor unit calculation in IMRT treatment plans.

Intensity-modulated radiotherapy (IMRT) treatment plan verification by comparison with measured data requires having access to the linear accelerator and is time consuming. In this paper, we propose a method for monitor unit (MU) calculation and plan comparison for step and shoot IMRT based on the Monte Carlo code EGSnrc/BEAMnrc. The beamlets of an IMRT treatment plan are individually simulated using Monte Carlo and converted into absorbed dose to water per MU. The dose of the whole treatment can be expressed through a linear matrix equation of the MU and dose per MU of every beamlet. Due to the positivity of the absorbed dose and MU values, this equation is solved for the MU values using a non-negative least-squares fit optimization algorithm (NNLS). The Monte Carlo plan is formed by multiplying the Monte Carlo absorbed dose to water per MU with the Monte Carlo/NNLS MU. Several treatment plan localizations calculated with a commercial treatment planning system (TPS) are compared with the proposed method for validation. The Monte Carlo/NNLS MUs are close to the ones calculated by the TPS and lead to a treatment dose distribution which is clinically equivalent to the one calculated by the TPS. This procedure can be used as an IMRT QA and further development could allow this technique to be used for other radiotherapy techniques like tomotherapy or volumetric modulated arc therapy.

Impacts on the deep-sea ecosystem by a severe coastal storm

Major coastal storms, associated with strong winds, high waves and intensified currents, and occasionally with heavy rains and flash floods, are mostly known because of the serious damage they can cause along the shoreline and the threats they pose to navigation. However, there is a profound lack of knowledge on the deep-sea impacts of severe coastal storms. Concurrent measurements of key parameters along the coast and in the deep-sea are extremely rare. Here we present a unique data set showing how one of the most extreme coastal storms of the last decades lashing the Western Mediterranean Sea rapidly impacted the deep-sea ecosystem. The storm peaked the 26th of December 2008 leading to the remobilization of a shallow-water reservoir of marine organic carbon associated with fine particles and resulting in its redistribution across the deep basin. The storm also initiated the movement of large amounts of coarse shelf sediment, which abraded and buried benthic communities. Our findings demonstrate, first, that severe coastal storms are highly efficient in transporting organic carbon from shallow water to deep water, thus contributing to its sequestration and, second, that natural, intermittent atmospheric drivers sensitive to global climate change have the potential to tremendously impact the largest and least known ecosystem on Earth, the deep-sea ecosystem.

Prediction of free imatinib concentrations based on total plasma concentrations in patients with gastrointestinal stromal tumours.

AIM: Total imatinib concentrations are currently measured for the therapeutic drug monitoring of imatinib, whereas only free drug equilibrates with cells for pharmacological action. Due to technical and cost limitations, routine measurement of free concentrations is generally not performed. In this study, free and total imatinib concentrations were measured to establish a model allowing the confident prediction of imatinib free concentrations based on total concentrations and plasma proteins measurements. METHODS: One hundred and fifty total and free plasma concentrations of imatinib were measured in 49 patients with gastrointestinal stromal tumours. A population pharmacokinetic model was built up to characterize mean total and free concentrations with inter-patient and intrapatient variability, while taking into account α1 -acid glycoprotein (AGP) and human serum albumin (HSA) concentrations, in addition to other demographic and environmental covariates. RESULTS: A one compartment model with first order absorption was used to characterize total and free imatinib concentrations. Only AGP influenced imatinib total clearance. Imatinib free concentrations were best predicted using a non-linear binding model to AGP, with a dissociation constant Kd of 319 ng ml(-1) , assuming a 1:1 molar binding ratio. The addition of HSA in the equation did not improve the prediction of imatinib unbound concentrations. CONCLUSION: Although free concentration monitoring is probably more appropriate than total concentrations, it requires an additional ultrafiltration step and sensitive analytical technology, not always available in clinical laboratories. The model proposed might represent a convenient approach to estimate imatinib free concentrations. However, therapeutic ranges for free imatinib concentrations remain to be established before it enters into routine practice.

Experience-dependent plasticity in the mouse barrel cortex

Abstract : Host-Cell Factor 1 (HCF-1) was first discovered in the study of the herpes simplex virus (HSV) infection. HCF-1 is one of the two cellular proteins that compose the VP16-induced complex, a key activator of HSV lytic infection. lncleed, when HSV infects human cells, it is able to enter two modes of infection: lytic or latent. The V`P16-induced complex promotes the lytic mode and in so doing the virus targets important cellular regulatory proteins, such as HCF-1, to manipulate the status of the infected cell. Indeed, HCF-1 regulates human cell proliferation and the cell cycle at different steps. In human, HCF-1 is unusual in that it undergoes a process of proteolytic maturation that results from cleavages at six centrally located 26 amino acid repeats called HCF-1pro repeats. This generates a heterodimeric complex of stably associated amino- (HCF-1n) and carboxy- (HCF-1c) terminal subunits. The absence of the HCF-1 N or HCF-1; subunit leads predominantly to either G1 or M phase defects, respectively. We have hypothesized that HCF-1 forms a heterodimeric complex to permit communication between the two subunits of HCF-1 involved in regulating different phases of the cell cycle. Indeed, there is evidence for such inter-subunit communication because a point mutation called P134S in the HCF-1N subunit in the temperature-sensitive hamster cell line tsBN67 causes, addition to G1- phase defects associated with the HCF-1n subunit, M-phase defects similar to the defects seen upon loss of HCF-1 function. Furthermore, inhibition of the proteolytic maturation of HCF-1 by deletion of the six HCF-1pro repeats (HCF-1Aimo) also leads to M-phase defects, specifically cytokinesis defects leading to binucleation, indicating that there is loss of HCF-15 function in the absence of HCF-1 maturation. I demonstrate that individual point mutations in each of the six HCF-1pro repeats that prevent HCF-1 proteolytic maturation also lead to binucleation; however, this defect can be latgely rescued by the presence of just one HCF-1pRO sequence in I-ICF»1. These results argue that processing itself is important for the HCF-1g function. In fact, until now, the hypothesis was that the proteolytic processing per se is more important for HCF-1C function than the proteolytic processing region. But I show that processing per se is not sufticient to rescue multinucleation, but that the HCF-lpm sequence itself is crucial. This discovery leads to the conclusion that the I-ICF-1pRO repeats have an additional function important for HCF-le function. From the studies of others, one potential function of the HCF-lrxo tepeats is as a binding site for O-link NAcetyl glycosamine tansferase (OGT) to glycosylate an HCF-1n-sunbunit region called the Basic region. This new function suggests the Basic region of HCF-1n is also implicated in the communication between the two subunits. This inter-subunit communication was analyzed in more detail with the studies of the Pl34S mutation and the residues 382-450 region of HCF-l that when removed prevents HCF-l subunit association. I demonstrate that the point mutation also leads to a binucleation defect in Hela cells as well as in the tsBN67 cells. In addition, the effect of this mutation on the regulation of HCF-1c activity seems to interfere with that of the HCF-lpgg repeats because the sum of the deletion of the proteolytic processing region and the point mutation surprisingly leads to re-establishment of correct cytokinesis. The study of the 382-450 HCF-lN region also yielded surprising results. This region important for the association of the two subunits is also important for both HCF-1c function in M phase and G1 phase progression. Thus, I have discovered two main functions of this region: its role in the regulation of HCF-lc function in M phase and its involvement in the regulation of G1/S phase ?- an HCF-1n function. These results support the importance of inter-subunit communication in HCF-1 functions. My research illuminates the understanding of the interaction of the two subunits by showing that the whole HCF-1n subunit is involved in the inter-subunit communication in order to regulate HCF-1c function. For this work, I was concentrated on the study of cytokinesis; the first phenotype showing the role of HCF-1c in the M phase. Then, I extended the study of the M phase with analysis of steps earlier to cytokinesis. Because some defects in the chromosome segregation was already described in the absence of HCF-1, I decided to continue the study of M phase by checking effects on the chromosome segregation. I showed that the HCF-1n subunit and HCF-1pro repeats are both important for this key step of M phase. I show that the binucleation phenotype resulting from deletion or mutation in HCF-1pro repeats, Pl34S point mutation or the lack of the region 382-450 are correlated with micronuclei, and chromosome segregation and alignment defects. This suggests that HCF«lç already regulates M phase during an early step and could be involved in the complex regulation of chromosome segregation. Because one of the major roles of HCF-1 is to be a transcription regulator, I also checked the capacity of HCF-1 to bind to the chromatin in my different cell lines. All my recombinant proteins can bind the chromatin, except for, as previously described, the HCF-1 with the P134S point mutation, This suggests that the binding of HCF-1 to the chromatin is not dependant to the Basic and proteolytic regions but more to the Kelch domain. Thus, if the function of HCF-ig in M phase is dependant to its chromatin association, the intercommunication and the proteolytic region are not involved in the ability to bind to the chromatin but more to bind to the right place of the chromatin or to be associated with the co-factors. Résumé : L'étude de l'infection par le virus Herpes Simplex (HSV) a permis la découverte de la protéine HCF-1 (Host-Cell Factor). HCF-1 est une des protéines cellulaires qui font partie du complexe induit par VP16 ; ce complexe est la clef pour l'activation de la phase lytique de HSV. Afin de manipuler les cellules infectées, le complexe induit pas le VPIG devrait donc cibler les protéines importantes pour la régulation cellulaire, telles que la protéine HCF-1. Cette dernière s'avère donc être un senseur pour la cellule et devrait également jouer un rôle de régulation lors des différentes phases du cycle cellulaire. Chez l'humain, HCF-1 a la particularité de devoir passer par une phase de maturation pour devenir active. Lors de cette maturation, la protéine subit une coupure protéolytique au niveau de six répétitions composées de 26 acides aminés, appelé HCF-1pro repeats. Cette coupure engendre la formation d'un complexe formé de deux sous-unités, HCF-1n et HCF-1c, associées l'une à l'autre de façon stable. Enlever la sous-unité HCF-IN ou C entraîne respectivement des défauts dans la phase G1 et M. Nous pensons donc que HCF-1 forme un complexe hétérodimérique afin de permettre la communication entre les molécules impliquées dans la régulation des différentes phases du cycle cellulaire. Cette hypothèse est déduite suite à deux études: l'une réalisée sur la lignée cellulaire tsBN67 et l'autre portant sur l'inhibition de la maturation protéolytique. La lignée cellulaire tsBN67, sensible à la température, porte la mutation Pl 345 dans la sous-unité HCF-1n. Cette mutation, en plus d'occasionner des défauts dans la phase G1 (défauts liés à la sous-unité HCF-1N), a aussi pour conséquence d'entrainer des défauts dans la phase M, défauts similaires à ceux dus a la perte de la sous-unité HCF-1c. Quant à la maturation protéolytique, l'absence de la région de la protéolyse provoque la binucléation, défaut lié à la cytokinèse, indiquant la perte de la fonction de la sous-unité HCF-1c. Au cours de ma thèse, j'ai démontré que des mutations dans les HCF-1=no repeats, qui bloquent la protéolyse, engendrent la binucléation ; cependant ce défaut peut être corrigé pas l'ajout d'un HCF-1pro repeat dans un HCF-1 ne contenant pas la région protéolytique. Ces résultats soutiennent l'idée que la région protéolytique est importante pour le bon fonctionnement de HCF-1c. En réalité jusqu'a maintenant on supposait que le mécanisme de coupure était plus important que la région impliquée pour la régulation de la fonction de HCF-1;. Mais mon étude montre que la protéolyse n'est pas suffisante pour éviter la binucléation ; en effet, les HCF-1pro repeats semblent jouer le rôle essentiel dans le cycle cellulaire. Cette découverte conduit à la conclusion que les HCF-1pro repeats ont sûrement une fonction autre qui serait cruciale pour la foncton de HCF-1c. Une des fonctions possibles est d'être le site de liaison de l'O-linked N-acetylglucosamine transférase (OGT) qui glycosylerait la région Basique de HCF-1n. Cette nouvelle fonction suggère que la région Basique est aussi impliquée dans la communication entre les deux sous- unités. L'intercommunication entre les deux sous-unités ai été d'ailleurs analysée plus en détail dans mon travail à travers l'étude de la mutation Pl34S et de la région 382-450, essentielle pour l'association des deux sous»unités. J'ai ainsi démontré que la mutation P134S entraînait aussi des défauts dans la cytokinése dans la lignée cellulaire Hela, de plus, son influence sur HCF-1c semble interférer avec celle de la région protéolytique. En effet, la superposition de ces deux modifications dans HCF-1 conduit au rétablissement d'une cytokinése correcte. Concernant la région 382 à 450, les résultats ont été assez surprenants, la perte de cette région provoque l'arrêt du cycle en G1 et la binucléation, ce qui tend à prouver son importance pour le bon fonctionnement de HCF-1n et de HCF-1c. Cette découverte appuie par conséquent l'hypotl1èse d'une intercommunicatzion entre les deux sous-unités mettant en jeu les différentes régions de HCF-1n. Grâce à mes recherches, j'ai pu améliorer la compréhension de l'interaction des deux sous-unités de HCF-1 en montrant que toutes les régions de HCF-1n sont engagées dans un processus d'intercommunication, dont le but est de réguler l'action de HCF-1c. J'ai également mis en évidence une nouvelle étape de la maturation de HCF-1 qui représente une phase importante pour l'activation de la fonction de HCF-1c. Afin de mettre à jour cette découverte, je me suis concentrée sur l'étude de l'impact de ces régions au niveau de la cytokinése qui fut le premier phénotype démontrant le rôle de HCF-1c dans la phase M. A ce jour, nous savons que HCF-1c joue un rôle dans la cytokinèse, nous ne connaissons pas encore sa fonction précise. Dans le but de cerner plus précisément cette fonction, j'ai investigué des étapes ultérieures ai la cytokinèse. Des défauts dans la ségrégation des chromosomes avaient déjà été observés, ai donc continué l'étude en prouvant que HCF-1n et les HCF-1pro repeats sont aussi importants pour le bon fonctionnement de cette étape clef également régulée par HCF-1c. J' ai aussi montré que la région 382-450 et la mutation P134S sont associées à un taux élevé de micronoyaux, de défauts dans la ségrégation des chromosomes. L'une des fonctions principales de HCF-1 étant la régulation de la transcription, j'ai aussi contrôlé la capacité de HCF-1 à se lier à la chromatine après insertion de mutations ou délétions dans HCF-1n et dans la région protéolytique. Or, à l'exception des HCF-1 contenant la mutation P134S, la sous-unité HCF-1c des HCF-1 tronquées se lie correctement à la chromatine. Cette constatation suggère que la liaison entre HCF-1c et chromatine n'est pas dépendante de la région Basique ou Protéolytique mais peut-être vraisemblablement de la région Kelch. Donc si le rôle de HCF-1c est dépendant de sa capacité â activer la transcription, l'intercommunication entre les deux sous-unités et la région protéolytique joueraient un rôle important non pas dans son habileté à se lier à la chromatine, mais dans la capacité de HCF-1 à s'associer aux co-facteurs ou à se placer sur les bonnes régions du génome.

Helical tomotherapy in the treatment of anal canal cancer: 3-year clinical experience

Objective: To report a single-center experience treating patients with squamous- cell carcinoma of the anal canal using helical Tomotherapy (HT) and concurrent chemotherapy (CT).Materials/Methods: From October 2007 to February 2011, 55 patients were treated with HT and concurrent CT (5-fluorouracil/capecitabin and mitomycin) for anal squamous-cell carcinoma. All patients underwent computed- tomography-based treatment planning, with pelvic and inguinal nodes receiving 36 Gy in 1.8 Gy/fraction. Following a planned 1-week break, primary tumor site and involved nodes were boosted to a total dose 59.4 Gy in 1.8 Gy/fraction. Dose-volume histograms of several organs at risk (OAR; bladder, small intestine, rectum, femoral heads, penile bulb, external genitalia) were assessed in terms of conformal avoidance. All toxicity was scored according to the CTCAE, v.3.0. HT plans and treatment were implemented using the Tomotherapy, Inc. software and hardware. For dosimetric comparisons, 3D RT and/or IMRT plans were also computed for some of the patients using the CMS planning system, for treatment with 6-18 MV photons and/or electrons with suitable energies from a Siemens Primus linear accelerator equipped with a multileaf collimator.Locoregional control and survival curves were compared with the log-rank test, and multivariate analysis by the Cox model.Results: With 360-degree-of-freedom beam projection, HT has an advantage over other RT techniques (3D or 5-field step-and-shot IMRT). There is significant improvement over 3D or 5-field IMRT plans in terms of dose conformity around the PTV, and dose gradients are steeper outside the target volume, resulting in reduced doses to OARs. Using HT, acute toxicity was acceptable, and seemed to be better than historical standards.Conclusions: Our results suggest that HT combined with concurrent CT for anal cancer is effective and tolerable. Compared to 3D RT or 5-field step-andshot IMRT, there is better conformity around the PTV, and better OAR sparing.

Psychomotor retardation in depression: bradykinesia or paucity of movement?

Objectives: Psychomotor retardation is part of Major Depression (MD) diagnosis criteria and has been assimilated to bradykinesia, even though there is a clear lack of objective measurement of motor activity in MD. We conducted a study to evaluate bradykinesia, posture and gait parameters in MD patients with an ambulatory system, allowing continuous motor measurements. Methods: Patients with DSM-IV MD and healthy controls matched for age and sex were asked to carry on with their usual activities while being recorded for 6 hours by a wireless autonomous ambulatory system, containing miniature gyroscopes, data-logger, battery and flash memory. allowing continuous recording of upper limbs movements (speed, amplitude and activity (% of time with movement)), posture (% of time standing, walking, lying or sitting) and gait parameters (speed, cadence, stance, double support, stride). Results: Hands activity was significantly lower in depressed patients, as compared to controls (MD: 40%, controls: 60%; p<0.05). Speed of hand movements (p= 0.13) and their amplitude (p=0.71) were similar to controls. MD patients had a trend to spend more time lying or sitting than controls (p=0.06) but did not differ in terms of any gait parameters. Conclusion: Patients with MD displayed less hand movements than controls and tended to spend more time lying or sitting over 6 hours, but did not differ in terms of speed and amplitude of movement, nor in gait parameters. These results suggest that psychomotor retardation classically described in MD might be the expression of a paucity of movement rather than a bradykinesia as observed in parkinsonism and might involved different (nondopaminergic) mechanisms.

Development of a simulation tool for torsional vibration of branched power transmission system

In this master's thesis a mechanical model that is driven with variable speed synchronous machine was developed. The developed mechanical model simulates the mechanics of power transmission and its torsional vibrations. The mechanical model was developed for the need of the branched mechanics of a rolling mill and the propulsion system of a tanker. First, the scope of the thesis was to clarify the concepts connected to the mechanical model. The clarified concepts are the variable speed drive, the mechanics of power transmission and the vibrationsin the power transmission. Next, the mechanical model with straight shaft line and twelve moments of inertia that existed in the beginning was developed to be branched considering the case of parallel machines and the case of parallel rolls. Additionally, the model was expanded for the need of moreaccurate simulation to up to thirty moments of inertia. The model was also enhanced to enable three phase short circuit situation of the simulated machine. After that the mechanical model was validated by comparing the results of the developed simulation tool to results of other simulation tools. The compared results are the natural frequencies and mode shapes of torsional vibration, the response of the load torque step and the stress in the mechanical system occurred by the permutation of the magnetic field that is arisen from the three phase short circuit situation. The comparisons were accomplished well and the mechanical model was validated for the compared cases. Further development to be made is to develop the load torque to be time-dependent and to install two frequency converters and two FEM modeled machines to be simulated parallel.

Lifecycle Information Management of a Mass Customized Product

Diplomityön tavoitteena on tutkia mitä uusia tiedonhallinnallisia ongelmia ilmenee, kun massaräätälöidyn tuotteen tuotetieto hallitaan läpi tuotteen elinkaaren, sekä miten nämä ongelmat voitaisiin ratkaista. Ongelmat ja haasteet kerätään kirjallisuuslähteistä ja massaräätälöintiprosessi yhdistetään PLM-vaiheisiin. Ratkaisua tutkitaan testaamalla kuinka standardit STEP ja PLCS sekä standardeja tukeva PLM järjestelmä voisivat tukea massaräätälöidyn tuotteen elinkaaren tiedonhallintaa. MC tuotteiden ongelmia ovat tuoterakenteen monimutkaisuus, jäljitettävyys ja muutosten hallinta läpi elinkaaren. STEP ja PLCS pystyvät kummatkin tahollaan tukemaan tiedonhallintaa. MC-tuotteen geneerinen tuoterakenne on kuitenkin manuaalisesti liittettävä elinkaaritiedon tukemiseen. PLM-järjestelmä pystyy tukemaan MC-tuotteiden elinkaarta, mutta koska toiminto ei ole järjestelmään sisäänrakennettuna, MC-tuotteiden tukemisen parantamisessa on edelleen haasteita.

Développement et validation d'une méthode de calcul indépendant des unités moniteur basée sur les simulations Monte Carlo pour la vérification des traitements IMRT

Résumé : La radiothérapie par modulation d'intensité (IMRT) est une technique de traitement qui utilise des faisceaux dont la fluence de rayonnement est modulée. L'IMRT, largement utilisée dans les pays industrialisés, permet d'atteindre une meilleure homogénéité de la dose à l'intérieur du volume cible et de réduire la dose aux organes à risque. Une méthode usuelle pour réaliser pratiquement la modulation des faisceaux est de sommer de petits faisceaux (segments) qui ont la même incidence. Cette technique est appelée IMRT step-and-shoot. Dans le contexte clinique, il est nécessaire de vérifier les plans de traitement des patients avant la première irradiation. Cette question n'est toujours pas résolue de manière satisfaisante. En effet, un calcul indépendant des unités moniteur (représentatif de la pondération des chaque segment) ne peut pas être réalisé pour les traitements IMRT step-and-shoot, car les poids des segments ne sont pas connus à priori, mais calculés au moment de la planification inverse. Par ailleurs, la vérification des plans de traitement par comparaison avec des mesures prend du temps et ne restitue pas la géométrie exacte du traitement. Dans ce travail, une méthode indépendante de calcul des plans de traitement IMRT step-and-shoot est décrite. Cette méthode est basée sur le code Monte Carlo EGSnrc/BEAMnrc, dont la modélisation de la tête de l'accélérateur linéaire a été validée dans une large gamme de situations. Les segments d'un plan de traitement IMRT sont simulés individuellement dans la géométrie exacte du traitement. Ensuite, les distributions de dose sont converties en dose absorbée dans l'eau par unité moniteur. La dose totale du traitement dans chaque élément de volume du patient (voxel) peut être exprimée comme une équation matricielle linéaire des unités moniteur et de la dose par unité moniteur de chacun des faisceaux. La résolution de cette équation est effectuée par l'inversion d'une matrice à l'aide de l'algorithme dit Non-Negative Least Square fit (NNLS). L'ensemble des voxels contenus dans le volume patient ne pouvant être utilisés dans le calcul pour des raisons de limitations informatiques, plusieurs possibilités de sélection ont été testées. Le meilleur choix consiste à utiliser les voxels contenus dans le Volume Cible de Planification (PTV). La méthode proposée dans ce travail a été testée avec huit cas cliniques représentatifs des traitements habituels de radiothérapie. Les unités moniteur obtenues conduisent à des distributions de dose globale cliniquement équivalentes à celles issues du logiciel de planification des traitements. Ainsi, cette méthode indépendante de calcul des unités moniteur pour l'IMRT step-andshootest validée pour une utilisation clinique. Par analogie, il serait possible d'envisager d'appliquer une méthode similaire pour d'autres modalités de traitement comme par exemple la tomothérapie. Abstract : Intensity Modulated RadioTherapy (IMRT) is a treatment technique that uses modulated beam fluence. IMRT is now widespread in more advanced countries, due to its improvement of dose conformation around target volume, and its ability to lower doses to organs at risk in complex clinical cases. One way to carry out beam modulation is to sum smaller beams (beamlets) with the same incidence. This technique is called step-and-shoot IMRT. In a clinical context, it is necessary to verify treatment plans before the first irradiation. IMRT Plan verification is still an issue for this technique. Independent monitor unit calculation (representative of the weight of each beamlet) can indeed not be performed for IMRT step-and-shoot, because beamlet weights are not known a priori, but calculated by inverse planning. Besides, treatment plan verification by comparison with measured data is time consuming and performed in a simple geometry, usually in a cubic water phantom with all machine angles set to zero. In this work, an independent method for monitor unit calculation for step-and-shoot IMRT is described. This method is based on the Monte Carlo code EGSnrc/BEAMnrc. The Monte Carlo model of the head of the linear accelerator is validated by comparison of simulated and measured dose distributions in a large range of situations. The beamlets of an IMRT treatment plan are calculated individually by Monte Carlo, in the exact geometry of the treatment. Then, the dose distributions of the beamlets are converted in absorbed dose to water per monitor unit. The dose of the whole treatment in each volume element (voxel) can be expressed through a linear matrix equation of the monitor units and dose per monitor unit of every beamlets. This equation is solved by a Non-Negative Least Sqvare fif algorithm (NNLS). However, not every voxels inside the patient volume can be used in order to solve this equation, because of computer limitations. Several ways of voxel selection have been tested and the best choice consists in using voxels inside the Planning Target Volume (PTV). The method presented in this work was tested with eight clinical cases, which were representative of usual radiotherapy treatments. The monitor units obtained lead to clinically equivalent global dose distributions. Thus, this independent monitor unit calculation method for step-and-shoot IMRT is validated and can therefore be used in a clinical routine. It would be possible to consider applying a similar method for other treatment modalities, such as for instance tomotherapy or volumetric modulated arc therapy.