Abstract

A tracking detector system has been constructed with an innovative approach to the classical multiwire proportional chamber concept, using contemporary technologies. The detectors, covering an area of 0.58 square meters each, are optimized for the application of muon radiography. The main features are high (>99.5%) and uniform detection efficiency, 9 mm FWHM position resolution, and filling gas consumption below 2 liters per hour for the nontoxic, nonflammable argon and carbon dioxide mixture. These parameters, along with the simplicity of the construction and the tolerance for mechanical effects, make the detectors a viable option for a large area muography observation system.

1. Introduction

Cosmic muons are highly penetrating particles, which are created by interaction of high energy cosmic particles impinging on the upper atmosphere of Earth. Muons are decay products of pions, which make up the majority of the charged particles reaching the Earth surface. Production mechanisms and interaction properties are extensively discussed in the literature [1]. The mean energy of muons is around 2 GeV, and such an average muon would be able to cross 3-4 meters of rock. This considerable penetrating power gives rise to the practical applicability of cosmic muon radiography, or “muography” in short.

There is a broad range of applications of muography. Very large objects can be imaged by measuring muon absorption, such as in the classical measurement by Alvarez et al. [2]. Notable, high impact examples are volcano imaging [35] and mapping the core of nuclear reactors [6]. Underground measurements have smaller background due to the natural shielding and can be exploited during ore mining [7] or natural cave search [811]. If the object of interest fits into a detector system, measuring the muon scattering [1215] can provide information about the internal structure.

When it comes to particle detectors for muon radiography, one of the key considerations is the flux of the signal particles. On the surface, the rate is in the order of one hundred particles per square meter per second, a relatively low value compared to high energy physics (HEP) experimentation. On the other hand, if a large object is to be imaged, the flux of the muons can be drastically reduced further [3]: this implies that large size detectors are mandatory for most meaningful measurements. An additional important point is the suppression of background against the relevant (transmitted) muon signal: this background turns out to be mostly low energy particles from various sources [16, 17].

From the above it follows that the requirements for muography detectors differ from those generally expected from high energy physics detectors (for the latter, see textbooks such as [18]). The key points may be summarized as follows:(i)Large Size at Sensible Cost. Due to the low and fixed cosmic rate, the statistics becomes the ultimate limiting factor of the measurements on large objects.(ii)Long Term Operation under Ambient Outdoor Environmental Conditions. Considerable temperature, pressure, and humidity variations need to be tolerated.(iii)Tracking Capability. In typical HEP conditions, the detector size is usually constrained by available space. On the contrary, most applications of muography require a specific angular resolution, which by increasing the measurable track length allows for weaker position resolution. However, clear identification of the straight muon tracks is mandatory to suppress low energy background.(iv)Rate Capability. The cosmic particle rate is low on the Earth surface; therefore none of the issues apply which characterize “high rate” detectors in HEP. Eliminating those design elements which are meant to improve rate capability potentially results in considerable simplification of the structure.(v)Sustainable and Safe Operation. Given the long continuous measurement times, low maintenance and low running costs are desirable, whereas outdoor installation requires safe and environmentally friendly operation.

Various detector technologies have been considered for muon radiography. The three notable groups are scintillators, emulsion, and gaseous detectors. So far scintillator-based trackers proved to be the most successful, applied for volcanology [1921] and mining [7]. Emulsions present very high position resolution [22], and hence the detector layers require small space [16, 23, 24]. The main drawback of this technology is that the information can be made available only after the film development, making real time measurements impossible. Gaseous detectors have the possibility to reach relatively large surfaces at reasonable cost, and therefore these tend to be the choice for the outer layers of most HEP experiments. Various groups proposed viable solutions with Resistive Plate Chambers [25], drift chambers [13, 14], and GEMs [15].

Our objective was to conceive a tracking detector system which is optimized for muography from the very beginning of the design. Based on our own experience with various gaseous detectors, an innovative approach to the classical Multiwire Proportional Chambers (MWPCs) [18] was chosen. In smaller sizes, MWPCs with simplified geometry may be used for educational purposes [26].

In order to address the “muography-requirements” above, design features include consistent material choice (Section 2) for improved thermal behavior, an internal support structure which causes little efficiency loss (Section 7) and improves mechanical stability (Section 8). As working gas, the mixture of argon and carbon dioxide was studied in detail (Section 10), being the most obvious choice from cost, environmental, and safety considerations.

The paper presents the detector construction, followed by relevant performance studies, for a sensitive area of about 0.58 m2. The key guiding principle was the simplification of the design, without compromise on detection efficiency and tracking capability: our studies below attempt to give a comprehensive performance report.

2. Construction

The detector uses a relatively large, 22 mm gas gap between two grounded cathode planes. The gas amplification takes place on an anode (sense, SW) wire plane nearly half way between the cathodes. In the same plane, field shaping wires (FW) are fixed as well parallel to the anode wires, as shown in Figure 1. There is another wire plane, for the so-called “pick-up wires” (PW), which run perpendicularly to the anode wires. The concept is similar to the CCC case [27, 28]: the anode wires are used as trigger and amplitude measurement, whereas the field shaping wires and the pick-up wires give the position information. The pick-up wire plane is 2 mm above one of the cathodes, as shown in Figure 1. The advantage of using such pick-up wires is cost consideration and simplification of the structure: an alternative would be using printed circuit boards (PCBs) with strips, however PCBs at large sizes tend to be expensive and their surfaces are not sufficiently flat. The anode wires for the specific detectors presented in this paper are 22 and 25 micrometers thick Au-plated tungsten, whereas both FW and PW are brass wires of 100-micrometer diameter. The two different anode wire diameters presented no practical difference except for a slight change in operational voltage.

In one detector panel, there are 64 field shaping wires (hence 63 sense wires), with wire spacing of 12 mm, covering a 768 mm long region. The pick-up wires cover the same length with 4 mm spacing (hence there are 192 of them), which are connected in groups of 3 (i.e., each of the 64 pick-up wire groups is 12 mm wide). All the wires are approximately 770 mm long. The front view outline is shown in Figure 2.

From the electronic connection point of view, the 64 field wires and the 64 pick-up wire groups are directly connected to the readout (see in Section 4), whereas the 63 sense wires are all interconnected. From the sense wire signal, two output signals are formed after preamplification and shaping: one is a TTL-compatible discriminated trigger signal, and the other is a direct analog output for amplitude measurement.

The detector is constructed out of fiberglass (glass fiber reinforced epoxy). Specifically the sidewalls are G-10, whereas the cathode planes and other printed circuit components (soldering points, electronics support) are standard FR-4. This way, all parts of the detector feature a consistent thermal behavior. The cathode planes are 1.5 mm thick FR-4 layers, fully covered with copper on both sides. The complete detector is glued with epoxy: the units are thus not repairable in case of broken wires (however, no such case was observed so far after a year of operation). The total weight of a detector panel is about 6 kg including electronics.

Given the relatively large area, the distance between the cathode planes needs to be kept at the nominal 22 mm with reasonable precision. In order to achieve this, internal support pillars were installed in a particular way between the wires, as shown in Figure 1. The pillars have a rectangular 3 mm by 4 mm cross section around the anode wire plane. These pillars, made out of insulating material (PLA) using 3D printing technology, do not touch any of the field and sense wires. The pillars were placed closer to the field wires, with around 0.5 mm clearance, which leaves approximately 2.5 mm distance from the sense wire. One can expect that, near a support pillar, the field lines still terminate on the sense wires; therefore the detector remains sensitive to the traversing particles, a feature which proves to be correct according to the studies presented in Section 7. For our particular detectors, 6 such pillars have been installed, according to Figure 2. The possibility of installing many of such pillars without introducing a large dead zone in the detector is of key importance for mechanical stability.

The detector layers were operated with mixture of argon and carbon dioxide, in proportion of 82 : 18, at various flow rates ranging from 1 to 10 liters per hour. Detailed studies on gas behavior are presented in Section 10.

3. Tracking Detector Design

A tracking system has been set up, using six individual detector layers, arranged each with 20 cm spacing (i.e., 1 m distance between the first and the last layer). The detector layers were mounted vertically, standing on the side without connectors (bottom side in Figure 2); that is, the anode wires were horizontal. The individual detectors are counted from 1 to 6 for the rest of the paper and will be referred to as “detector layers” for clarity, as opposed to the “tracking system” denoting the complete stack. For most of the measurements, the detector layers shared the same anode high voltage. The gas line was connecting the layers in series: the gas entered to 1 and was vented to the atmosphere from 6.

The setting with vertical detector layers resembles most to the application in volcanology, where the object is imaged from the “side,” and where, because of the reduced muon flux close to the horizon, the system is most sensitive to background. Such dedicated issues are however outside of the scope of the present paper: for the tracking system we did not use absorbers, which are otherwise mandatory for such arrangements [5] and will be subject to later detailed studies. The image of the actual system with vertical detector layers is shown in Figure 3, including two additional nonfunctional layers.

For various performance studies, we applied a similar arrangement, but with horizontal layers inside the laboratory: that is, capturing the most abundant nearly vertical muons, high statistics mappings of the detector parameters became possible.

The six detector layers are sufficient to identify clear particle tracks (with at least 5 measured points) and allow for the actual detector performance studies. We have used tracking algorithms from our own earlier systems [10]. From the tracking information one can constrain the detector geometry (by requiring straight tracks on average), which is needed for a more precise track finding. The tracking system has shown a consistently excellent performance, enabling a reliable characterization of the individual detector layers using cosmic muons.

4. Readout and Data Acquisition

The detector readout was adapted from our own earlier system [11], with a single Raspberry Pi microcomputer as the core of the data acquisition [29]. An event sequence starts with the arrival of the trigger signals: if three coincident triggers appear (from any of the 6 layers), the readout will wait for a predefined delay time and then issue a “master” trigger pulse. At the instance of this master trigger, the position information is recorded from discriminated signals of all the channels of the field wires and the pick-up wires. Also the master trigger time defines when the sense wire amplitude measurement, via 10 bit ADCs, takes place.

The total data from the tracking system is moderate: there are 6 × 2 × 64 = 768 bits from the position information, and including the converted ADC values as well as recording the trigger pattern, the raw data is still around 100 bytes per event. This explains why a single Raspberry Pi can handle the DAQ task.

The input signals from the wires are protected from external electronics noise as much as possible, and since the grounded cathodes act as a Faraday cage, a reasonable noise performance has been observed consistently, with input equivalent noise below 3 × 104 electrons (5 fC). The trigger level from each of the detector layers is set approximately to 2 × 105 electrons (32 fC) on input; for the field shaping and pick-up wires the discrimination threshold was also set to 2 × 105 electrons (measured as charge on the sense wires).

5. Efficiency of Cosmic Muon Detection

A key figure of merit of a particle detector system is the tracking detection efficiency, that is, referring to a single layer, the probability that a valid particle hit (position) is actually detected. The tracking efficiency is defined in the following way for the rest of the paper: considering a specific detector layer, a particle trajectory is formed from the 5 detectors excluding the one under study. If this 5-point pattern is a unique, straight track from both directions (see Section 6), one checks if there is a hit at the expected position (within 5 sigma of the position resolution) in the detector layer under study. The probability of finding the hit is defined as the “tracking efficiency.” In a similar way, one can check if the specific detector layer has provided a trigger signal for the event: this is defined as the “trigger efficiency.” Trigger and tracking efficiencies are found to be well above 90% for the nominal operational conditions.

First, the tracking and trigger efficiency has been determined as a function of the sense wire voltage. A clear plateau is observed, shown in Figure 4. The efficiency reaches around 95% (>99%) at 1650 V (1700 V) for the detector layers constructed with 25 micron anode wires. The plateau continues to 1800 V, where stable operation has still been observed.

The efficiency values shown in Figure 4 are consistent with the expectation from the signal amplitude spectrum, shown in Figure 5 for a sense wire voltage of 1700 V for the first tracking layer. The spectrum of all tracks with at least 5 points is overlayed with the spectrum of 6-point tracks. The difference between these two sets, in case of this nominal high efficiency operating point, includes those tracks which miss the first layer and hence result in a 5-point track. The approximate detection level is 120 ADC units which has also been indicated with an arrow. The horizontal scale of Figure 5 is expressed in ADC units, which have been calibrated by direct charge injection: 1 ADC unit corresponds to 1.8 × 103 e, or 0.3 fC.

The true amplification gain, that is, the multiplication factor of the avalanches, has been estimated measuring the mean amplitude on the sense wires. Taking into account that the mean number of electrons over the 2 cm of the drift gap is around 200 in the applied gas mixture [18], one finds that the detectors with 25 micron anode wires at 1700 V (full efficiency) have a gain of about 7 × 103, a conveniently low figure.

The trigger signal is created at the point when the input signal reaches the trigger threshold. There is a certain time delay needed after this to reach the highest amplitude on the field wire and pick-up wire discriminators and hence to reach optimal tracking efficiency. The tracking efficiency as a function of this trigger delay is shown in Figure 6, with a maximum around 2-3 μs consistent with the peaking of the signal. The data has been taken keeping detectors 1, 2, 5, and 6 at 1700 V and reducing the voltage on the middle detectors, 3 and 4, to the values indicated (1600 V and 1550 V). This procedure ensures reasonably high trigger efficiency of the tracking system (from four detector layers) and a measurable efficiency (i.e., well below 100%) for the inner layers. The setting for all other measurements uses a fixed 2 μs delay.

6. Position Resolution

With the 6-layer tracking system setup, the position resolution of the layers can be quantified. Figure 7 shows the distance between the measured hit in one of the middle detector layers and the extrapolation from the 5-point track (i.e., the fitted track excluding the specific layer). One finds that the position resolution (RMS) is = 3.73 mm in the field wire and = 3.92 mm in the pick-up wire direction, that is, around 9 mm FWHM. This is roughly consistent with the 12 mm spacing expectation (). The tracking algorithm incorporates this position resolution to find a set of points on 5 or 6 of the chambers which form a straight line from both directions.

7. Dead Zones due to Support Pillars

The individual detector layers, as seen above, are highly efficient except for the support pillars, which deserve a closer look. In the tracking system configuration with horizontal detectors, the local efficiency has been determined, as shown in Figure 8. For this measurement, nearly vertical tracks were selected, with less than 200 mrad relative to the zenith. The pixel size in this map is one half of the readout unit (6 mm), and one immediately notes the surprisingly small reduction (the measured efficiency does not go to zero due to the fact that the pillars are smaller than the tracking detector position resolution). Taking a 3 cm wide strip over the pillars (indicated in the figure by two dashed lines), the mean efficiency is indicated in Figure 9. To first order, one can conclude that the detector is sensitive for any muons which enter the gas volume around the pillars, or in other words, the presence of the support pillars does not distort the field structure to the extent to considerably reduce efficiency. There seems to be some fluctuation on the apparent loss caused by the pillars, which may have to do with the details of the actual positioning.

8. Mechanical Stability

The mechanical stability and required tolerances of the detector are a key question for practical outdoor applications. In the studies below, some of the mechanical effects on the detector layers are quantified in relation to their performance. Despite the simplicity of the structure, it turned out that most of the relevant mechanical variations (e.g., difference between vertical or horizontal position) give insignificant changes; therefore we decided to expose the detectors to considerable mechanical effects, outside of the normally expected conditions.

The gain map, that is, the mean signal amplitude over the sensitive surface under reference conditions, is shown in Figure 10(a). In the first study, the detectors were ran with an increased inner gas pressure of 1.5 mbar above the ambient atmosphere. This results in a slight bulging; the inner support pillars in this case act against the forces on the cathode layers. The corresponding gain map is shown in Figure 10(b), demonstrating that the gain drop under this condition is below 10% (and hence, under normal operating conditions, one expects practically no efficiency loss).

The second study involved placing a weight directly on the detector. In this case, the detector was supported at the four corners, and two weights, 1 kg each, were placed at the positions indicated in Figure 10. One position was between the pillars; the other was approximately on one of the pillars. The gain increase (due to the reduction of the distance between the cathode layers) is apparent in the first case (by about 15%), whereas at the position of the pillar, the change is smaller. There is no effect observable outside of the weight positions (e.g., no global detector bending).

One can conclude that the mechanical structure, despite the simplicity, is very robust: the internal support pillars ensure the constant distance between cathode planes and hence allow for acceptable low gain variation.

9. Rate Capability

In the standard setting, the tracking system was operated in outdoor conditions on ground level, which is the maximum available cosmic muon rate. In order to quantify if a higher local rate is tolerable without efficiency loss, a beta source was placed close to one of the layers. The data was taken with a trigger condition requiring a coincidence of three detector layers, just as earlier; that is, master triggers were issued independently from the source: beta particles from cannot fully penetrate a single detector layer. The native trigger count rate (including natural radiation sources and detector noise) of the specific detector layer increased by 3.5 kHz (from around 1 kHz without source in horizontal position) from this concentrated irradiation covering about 20 cm2. Correspondingly, the number of hits appearing in triggered events (events in the sense of the tracking system) increases locally over the homogeneous distribution dominated by cosmic events, as shown in Figure 11. The reason why one can observe a modest local increase of hits is due to the finite time resolution of the detector.

The tracking efficiency has been determined along strips of 36 mm width (indicated by the dashed lines in Figure 11) in both directions. The data was taken at an anode voltage of 1700 V, that is, at efficiency of 99.5%. According to Figure 12, there is no efficiency loss observed for the source region (shown with the peaked continuous line).

The total observed count rate for a single layer corresponds to 0.2 Hz/cm2 (1 kHz over 0.5 m2 area), whereas in our experiment the induced local rate is above 170 Hz/cm2 (3.5 kHz over 20 cm2 area), corresponding to a three orders of magnitude increase. From the result, one can conclude that rate capability limitations, if any, are not relevant for the presented detector design in cosmic radiography applications.

10. Required Gas Flow and Tolerance on Gas Purity

The detector runs with industrial grade argon and carbon dioxide mixture in 82 : 18 proportion (Ferroline C18). The total volume of a single detector panel is around 13 liters, which implies long gas equilibration times for typical operational gas flow rates of 1–3 liters per hour.

In order to understand the behavior of the gas system, we performed two studies. One of these involved direct injection of ambient air into the gas line, whereas the other was part of the realistic outdoor studies presented in Section 12. In the first case, the gas flow was set to 1 liter per hour, and the detectors were connected in series to the gas line (sequentially from 1 to 6). Using a gas bypass between detector layers 1 and 2, 166 cm3 ambient air has been injected into the gas line, which then entered the second layer and later gradually moved along the subsequent panels. This way, the signal amplitude change relative to the first detector layer quantifies the effect of air (oxygen) contamination.

Figure 13 shows the mean anode wire signal amplitude in the layers as a function of time , normalized to the mean signal in detector layer 1 as reference before the gas injection (at in the figure). The pattern is very clear, with a sharp decrease in the second layer, followed with an exponential purging.

Assuming complete mixing of the air bubble in the gas volume, indeed an exponential relaxation is expected for the mean signal amplitude , described bywhere the time constant is fixed by the ratio of the detector volume and the flow rate and found to be hours. The only free parameter , the relative gain drop, is determined to be %. The function is shown in Figure 13 with the continuous line. In fact, the assumption of “complete mixing” predicts the gain drop in the subsequent detector layers, without any free parameters:for , where is the number of the detector layers, and with we get back the pure exponential behavior. Figure 13 confirms this behavior, indicated with the lines: the air bubble walks along the detectors.

The gain drop of 34% corresponds to 1.3% (13000 ppm) air contamination or, assuming that only the oxygen content of air counts, 2500 ppm of contamination level. One can conclude that the gas distribution inside the detector system is fully predictable, with the fixed time constant determined as the detector volume filling time. Oxygen contamination of 1000 ppm, leading to less than 15% gain drop (and hence no efficiency loss), thus seems to be tolerable, which is indeed a relaxed requirement for gas quality.

11. Background and Noise Performance

The instrumental background induced by the detector system is of vital importance from the point of view of practical muon radiography. Vaguely, we can call “background” anything which looks like a muon that traversed the object under study, while it actually did not. The most relevant component of such backgrounds is given by (low energy) particles which actually reach the detector [17]. This component can only be eliminated by a careful detector design and will be subject of future detailed studies; however, at this point we aim at quantifying instrumental artifacts (background which may not even be connected to any particle depositing charge in the detector). One possibility is to apply a “random” (uncorrelated oscillator) trigger and observe patterns resembling trajectories. In the first place, it has been verified that the pure “electronics” noise is exceedingly small: with the high voltage off, simply there were no triggers, and there were no observable trajectories. The probability of any of the single channels (on FW or PW) firing is in the order of ; that is, the chance to have 5 or 6 of them in a line to fake a track is negligible.

In order to have a relevant background estimate, the tracking system was operated in a fully efficient configuration () with vertical detector layers and was running for a day with oscillator trigger at 110 Hz. For the background estimation two data sets are compared: this “random” (uncorrelated oscillator) trigger data, which contains reconstructed tracks and can be used to estimate the “background” flux , and the normal trigger data, resulting in the measured true particle flux . Standard procedures to calculate are described in [9]. is proportional to the number of observed tracks per event:where the number of tracks is measured over time of , whereas the number of triggered events divided by is the trigger rate . The coefficient of proportionality includes the detector geometry. is the mean track number per event. Most triggered events do not contain acceptable (5 or 6 point) tracks due to the weak trigger condition; therefore the background tracks having 5 or 6 points add up to the “true” track set without the possibility to be suppressed as double track. In the data taken with oscillator trigger, the number of identified tracks is denoted by . The probability of random tracks per event can be used to estimate the background flux :where is the event number taken with the oscillator (random) trigger. The background flux has been evaluated as a function of elevation angle and is shown in Figure 14. There is a clear minimum at the horizon (elevation angle 0), reaching a background flux value of around . The background flux happens to be around 4 orders of magnitude below the true measured flux (also indicated in the Figure). One has to note that this measured flux has nothing to do with the flux of cosmic muons (which should reduce to lower values near the horizon): the measurement includes a considerable and undetermined low energy background.

The key conclusion of the observations above is that there is no apparent uncorrelated track background beyond that expected from the truly random appearance of otherwise physical tracks; such tracks are expected to be suppressed along with the low energy background component in true muography tracking systems.

12. Experience under Ambient Outdoor Conditions

For testing under ambient outdoor conditions, the tracking system with vertical detectors was installed on a support frame and placed on the flat roof of a building at the Wigner Research Centre campus. To allow for considerable environmental influence, the system was only lightly covered with a box, thus receiving strong sunlight in the early afternoon. In Figure 15, the outside temperature is shown on (a), featuring more than 25°C daily temperature variations. During the 6 days of continuous data taking in our example, the system sustained a light raining on the 5th day; otherwise the weather was clear.

The change of the temperature by 30°C corresponds to nearly 10% change in absolute temperature and hence in gas density. The system thus cannot run at constant anode voltage but needs to be compensated to first order: in our system, the anode HV was reduced by 1 V for a 1°C temperature increase (from the nominal 1700 V at 25°C). Though the compensation is far from perfect (possibly containing nonlinear terms), it has been verified that this ensures only a slight, 10% mean amplitude variation (demonstrated in Figure 15(b)).

The key figure of merit for the detector is the tracking efficiency, shown in Figure 15(c) for two specific detector layers. The measured tracking efficiency is reliably high (>99.5%) and stable for the whole period (shown only for two detector layers but found consistently for all panels).

The detector system was running with a constant, 2 liters per hour gas flow rate. This presents an interesting issue during decreasing temperatures. The total tracking system gas volume was about 100 liters. Assuming that the temperature decreases by 6°C in one hour, the 2% relative change in absolute temperature causes the total gas volume to reduce accordingly, that is, by 2 liters. In this case, actually air may be sucked in from the outside atmosphere to the last detector. In order to avoid this, a long buffer tube has been added after the last detector, with sufficiently large (in our case 1 liter) total volume.

13. Summary and Conclusions

A new approach to the classical MWPCs has led to a highly efficient, low weight, mechanically and environmentally tolerant tracking detector, which by its simple and cost efficient construction can be an ideal part of large size muon radiography detectors. An internal support pillar structure has been shown to present a negligible dead zone. The tracking system has been demonstrated to tolerate broad temperature variations maintaining very high detection efficiency (>99.5%) of the detector layers. The position resolution of better than 4 mm RMS (9 mm FWHM) is compatible with the segmentation. Gas flow rate using nonflammable, nontoxic Ar and CO2 mixture of 1-2 liters per hour is sufficient for the operation, with a well predictable gas flow pattern. Future muography tracking systems will be completed with absorber (scatterer) layers and are intended to be applied for large size object imaging, such as the internal structure of volcanos.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

Support was provided by members of the REGARD group of the Wigner Research Centre for Physics of the Hungarian Academy of Sciences. The authors are grateful to Professor H. K. M. Tanaka for the useful discussions. This work was supported by the Momentum (Lendület) grant of the Hungarian Academy of Sciences under Grant no. LP2013-60.